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Piezoelectric Crystals in Oscillator Circuits 
By I. E. FAIR 


12.00 INTRODUCTION 

STUDY or an explanation of the performance of a piezoelectric crystal 

in an oscillator circuit involves a study or explanation of oscillator 
circuits in general and a study of the crystal as a circuit element. Nicolson’ 
appears to have been the first to discover that a piezoelectric crystal had 
sufficient coupling between electrical electrodes and mechanical vibratory 
movement so that when the electrodes were suitably connected to a vacuum 
tube circuit, sustained oscillations were produced. In such an oscillator 
the mechanical oscillatory movement of the crystal functions as does the 
electrical oscillatory circuit of the usual vacuum tube oscillator. His 
circuit is shown in Fig. 12.1. Cady’ independently though later made the 
same discovery, but he utilized it somewhat differently and expressed it 
differently. He found that when the electrodes of a quartz crystal are 
connected in certain ways to an electric oscillator circuit, the frequency is 
held very constant at a value which coincides with the period of the vibrat- 
ing crystal. He made the further discovery that due to the very sharp 
resonance properties of the quartz crystal, the constancy in frequency to be 
secured was far greater than could be obtained by any purely electric 
oscillator. 

The development of analytical explanations of the crystal controlled 
oscillator came along rather slowly. Cady explained the control in terms 
of operation upon the electrical oscillator to which the crystal was attached. 
He said that the ‘‘capacity”’ of the crystal changes rapidly with frequency 
in the neighborhood of mechanical resonance, even becoming negative. 
This “‘capacity”’ connected across the oscillator tuned circuit or in other 
places prevented the frequency from changing to any extent, as any fre- 
quency change caused such a “‘capacity”’ change in the crystal as to tend to 
tune the circuit in the other direction. Cady, however, devised one circuit, 
Fig. 12.2, in which no tuned electrical circuit was used, but he confined his 
explanation to ‘‘a mechanically tuned feedback path from the plate to the 
grid of the amplifier”. Pierce*® came along later with a two-electrode crystal 
connected between plate and grid, and no tuned circuit, and also with a 
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Fig. 12.2—Cady’s oscillator circuit using a crystal as a “mechanically tuned 
feedback path” 
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Fig. 12.3—Equivalent electrical circuit of a piezoelectric crystal near its 
resonant frequency 


two-electrode crystal connected between grid and cathode and no tuned cir- 
cuit, where the operation would not be satisfactorily explained by Cady’s 
method. His circuits would require the crystal to exhibit inductive react- 
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ance, rather than the capacitance Cady spoke of. Miller* also produced a 
circuit with a two-electrode crystal connected between grid and cathode but 
with a tuned circuit in the plate lead, which circuit required the crystal to 
provide inductive reactance. 

It was not until after Van Dyke’ showed that the crystal could be repre- 
sented by the circuit network of Fig. 12.3 that it was possible to explain 
these various phenomena. With this view of the crystal, and using the 
differential equation method of circuit analysis, Terry® pointed out that, 
as with electrical oscillators, the frequency is not completely governed by 
the resonant element, in this case the crystal, but is influenced somewhat 
by the circuit elements. The circuit as a whole is quite complex and the 
equations are difficult to use. Wright’ and Vigoureux” also made analyses 
of the Pierce type oscillator. Because of the complexity of the equations, 
the frequency, amplitude, or activity are not computed directly, but the 
effects of the circuit variables are analyzed in a qualitative manner and 
the results compared with experimental data. 

Oscillators employing crystals may be classified in a number of ways. 
One classification is based upon whether or not the circuit without the 
crystal is in itself an oscillator. If it is, the oscillator is called a “crystal 
controlled” oscillator. If it is not, it is called a ‘‘crystal” oscillator. All 
of Cady’s oscillator circuits, except the one shown in Fig. 12.2, are of the 
first named class. This type of circuit will oscillate at a frequency deter- 
mined by the tuned circuit if the crystal becomes broken or disconnected, 
or if high resistance develops in the crystal, or if the electric tuned circuit 
should become tuned too far from the resonant frequency of the crystal. 
This property at times is an advantage and at other times a disadvantage. 
This type of circuit will oscillate under control of the crystal with much 
less active crystals than most of the other types. 

Nicolson’s, Pierce’s, Cady’s of Fig. 12.2 and Miller’s oscillators belong to 
the second named class. They will cease oscillating if the crystal breaks, 
develops high resistance or is disconnected. Failure of the oscillator to 
function at all then serves as a warning that something has happened to 
the crystal. 

This second named class of crystal oscillators has been used much more 
than the first named. The crystal is the principal frequency determining 
element in the circuit. Often there are required only resistances, or re- 
sistances and an inductance, as the other elements to embody along with 
the vacuum tube and crystal. The simplicity, low costs, and usually no 
tuning, have made this class attractive. Most analytical studies of oscilla- 
tor circuits have been made upon this class. For that reason the discussion 
in this chapter will be limited to this class. 

An analytic study of the crystal oscillator can readily start by looking 
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upon the oscillator as consisting only of inductances, capacitances, and 
resistances, along with the vacuum tube. The crystal is replaced by the 
proper circuital elements arranged as in Fig. 12.3. ‘This circuit or equivalent 
of the crystal is that of a series resonant circuit having capacitance parallel- 
ing it. The circuit will show both phenomena of series resonance and 


parallel resonance, the two frequencies being very close together. By 


making suitable measurements on a crystal, the magnitudes of the in- 
ductance, resistance, and the two capacitances can be determined. It is 
usually found that the series inductance is computed as hundreds or thou- 
sands of henries, and the series capacitance is a small fraction of a micro- 
microfarad. The magnitudes of the inductance and capacitance are beyond 
what it is possible to construct in the usual forms of building inductances 
and capacitances. This accounts for its superior frequency contro! 
properties, 

Although reducing the crystal to an equivalent electrical circuit provides 
one notable step in understanding the performance of the crystal oscillator, 
it does not readily lead to a full understanding. The electric oscillator in 
itself is not fully and completely analyzed in all its ramifications, although 
it has been under study for over 25 years. These studies have been mathe- 
matical and experimental in character, but in all cases it appears there have 
been approximations of some kind, made because the variable impedance 
characteristics both of the plate circuit and the grid circuit of the tubes 
did not lend themselves readily to a rigorous analysis. The earlier investi- 
gations assumed a linear relation between grid voltage and plate current 
and assumed constant plate impedance. Later investigations brought in 
further elements and further variables, the different investigators attacking 
the problem in different ways and attempting to prove different points. 
By this means a large number of factors in oscillators have been ascertained 
to a first degree of approximation so that a qualitative review of the per- 
formance of the electric oscillator is very well known. It is the quantitative 
view upon the first order magnitude which is still difficult or uncertain. 
This is particularly true of the crystal oscillator because of the slightly 
different circuit. 

It is proposed, therefore, in this paper to cover briefly a number of the 
studies on crystal oscillators so as to point out the different modes of attack 
and the different behavior points in the oscillators which the various investi- 
gators have studied. After covering these points, there will be discussed 
the frequency control properties of the crystal and the frequency stability of 
crystal oscillators. The performance of the crystal in the oscillator with 
respect to activity is then treated. There will be introduced two new yard- 
sticks for measuring or indicating crystal quality, one called “figure of 
merit” and the other called “performance index.” These are related to 
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(he crystal constants and paralleling capacitances which are usually involved. 


rhey will be defined and their method of use and application in oscillators 
will be pointed out. 


12.10 SoLuTion BY DIFFERENTIAL EQUATIONS 


The most direct method of determining the oscillating conditions in a 
circuit is to analyze the differential equation for the current in some particu- 
lar branch of the circuit. The relations existing between the coefficients 
determine whether the current builds up, dies out, or is maintained at a 
constant value and frequency. Unfortunately the equations resulting from 
the application of this method to the crystal oscillator circuit are quite 
complicated. However, lower order differential equations result from the 





























Fig. 12.4—Equivalent circuit of oscillator with crystal connected between 
grid and plate 


application of this method to similar electric oscillator circuits, and certain 
qualitative information obtained from the latter is applicable to crystal 
oscillators. Thus Heising’s’ analysis of the Colpitts and Hartley circuits 
gives much information directly applicable to the Pierce and Miller types 
of crystal oscillators. From this the circuit conditions necessary for oscilla- 
tions to exist and the effect of certain circuit variables upon the frequency 
are ascertained. The more complex qualitative view is given by Terry” 
who shows the relations of the coefficients of linear differential equations 
of the 2nd, 3rd, and 4th orders, and applies them to the analysis of three 
common types of crystal oscillator circuits. The resulting equations, 
together with certain qualitative information regarding their interpretation, 
are repeated here. In making this analysis the grid current is disregarded 
and the static tube characteristic is considered linear. 
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The equation is the same for the three types of circuits considered and is 
derived for the current 7; , in Figs. 12.4 and 12.5, although it may be set up 
in terms of any of the currents or voltages existing in the circuit. It is 
of the form 

4+ ° 
d (al d’ i; 


2 
dn tt ap 


2. 
» ah 


ili 
+ Ps ap 


+ P3— + Pui, = 0 (12.1) 
dt 
The P coefficients are functions of the circuit elements and are defined for 
each type of circuit in the following sections. 
The solution of (12.1) normally represents a doubly periodic function 
arising from the two coupled antiresonant meshes (a) and (b). The norma! 
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Fig. 12.5—Equivalent circuit of oscillator with crystal connected between 
grid and cathode 

















modes of oscillation consist of two currents in each mesh with frequency 
and damping factors 8; and a1, 82 and a: respectively. 

The conditions for undamped oscillation as derived from the general 
equation (12.1) are expressed in terms of the coefficients by 


Py _ Ps VP} — 4Py (12.2) 
P; 2 


and the angular frequencies are 


ene 
Ps. wa (12.3) 


where the plus sign gives the condition for one damping factor to be zero 
and the minus sign that for the other to be zero. 
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The frequency at which oscillations are maintained is determined by the 
required phase relation of voltages applied to the tube. With crystal from 
grid to plate, as in Fig. 12.4, the phase difference of grid and plate voltages 
is such that the circuit oscillates at only one of the normal modes, and with 
crystal connected between grid and cathode, as in Fig. 12.5, it oscillates at 
the other only. 


12.11 
With the crystal connected between the grid and plate of the tube, as in 
Fig. 12.4, the coefficients of the general equation (12.1) are 
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The uncoupled damping factors, a, and a, the uncoupled undamped 
angular frequencies, 8, and % , and the coupling coefficient + may be intro- 


duced as follows 
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Note that C, is the total capacitance across L; and R,, and C; is the total 
capacitance across Ly and Re. 
The coefficients of (12.4) become 


1 
>» or ’ 
P, 2(Qa a ap) + FA 


P r 1 2 
P, Ba + 4acayn + Bs + > (aa + a) Z 
R, Cy 


i 1 * 1 1 
P3 2(an Ba + aaBs) + R. (B, + 4a, as) | 
p 


Go LOG 
P= eet}i - 2+ 2(5 - 
a x, Cy 





The coefficients as given by (12.5) satisfy (12.2) and (12.3) only when the 
plus sign is used. 
The equations are simplified by dividing through by 8; thus 


Pe > \2 ? 
at a/ (By - 
ee Ba Ba Ba 


2 
: Pr» 2 4P, 
2 + V/ (3) ~ gt (12.7) 


2 


which gives the ratio of driven frequency of the crystal to its undriven value. 
The common variable R, must satisfy both (12.6) and (12.7). The method 
of computing the frequency would be to solve for R, in (12.6) and substitute 
in (12.7). However, the equations are too complicated a function of R, 
for this to be practical. Terry solved them graphically by plotting (12.6) 
and (12.7) as functions of R, for assigned values of the circuit, and the inter- 
section of these curves gave the frequency for the different circuit conditions. 
The results are shown in Fig. 12.6. The G-P curves show the frequency 
change as a function of plate circuit tuning for the grid to plate connection 
of the crystal. 


12.12 Crystat BETWEEN GRID AND CATHODE 


With the crystal connected between the grid and cathode of the tube, 
the circuit is as shown in Fig. 12.5. The coefficients of equation (12.1) 
are as follows: 
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With the substitution of uncoupled frequencies, damping factors and 


coupling coefficient as described in the previous section, they become 
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These equations of conditions for oscillation in this case satisfy (12.4) 
and (12.5) only when the minus sign is used. That is 
Po, P, a P: — 4P; 
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Again dividing by 83 to obtain the frequency as a ratio of driven to undriven 
crystal frequency, we have 


Pr 2) 4Py 
P3 a 3, a2 (12.12) 
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Fig. 12.6—The oscillating frequency as a function of the plate circuit frequency for the 
crystal connected grid to plate (G-P) and grid to cathode (G-C) 





The frequency change as a function of plate circuit tuning was determined 
graphically in the manner described in section (12.11) and the curves are 
shown in Fig. 12.6 as the G-C curves. 


12.13. RersIstANCE LOAD Circuit 


This is a special case of Plate-Grid connection of the crystal described in 
section (12.11) in which the plate circuit consists of a capacitance and re- 
sistance in parallel. This is a very common Pierce type of oscillator circuit 
and has the advantage that no tuning adjustment is necessary when using 
crystals of different frequencies. 
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Since this circuit is singly periodic, the differential equation for 7 is 
F of the third order and is derived from (12.1) by setting the plate inductance 
L, of the P coefficients equal to zero. The general equation then becomes 
da 2 di; 


Ie 
du? ah 


= RG. + BLO + BUC _— 





1 1 1 1 1 
aly a R, — GC rors | 


| With the substitution of the uncoupled damping factors and frequencies, 
> (12.15) becomes 


Py = 
P, = , = (12.16) 


P3; = 





RG Ric * 
The frequency as obtained from (12.14) is 


8° (12.17) 


with the conditions for oscillation 
P, — (12.18) 
1 


obtained by setting the damping factor a equal to zero. The ratio of driven 
to undriven frequency is obtained by dividing (12.17) and (12.18) by Bs. 
That is 

e Pr P; 


st Bt RP, — 


12.14 INTERPRETATION OF THE EQUATIONS 


It is learned from this analysis that the frequency of oscillation while 
governed principally by the frequency of the crystal also depends upon all 
the constants of the circuit. The effect of the plate circuit impedance is 
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12.7—Calculated increase in mean plate resistance against capacitance of the 
oscillatory circuit 
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Fig. 12.8—Experimental curves, showing the influence of interelectrode capacitances on 
the frequency 

shown in Fig. 12.6. It is pointed out that the effect of the crystal resistance 

R, is to decrease the frequency for the G-C connection and increase the 

frequency for the G-P connection. The discrepancy between the measured 
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and experimental values shown on the curves is attributed to the difference 
between chosen and actual value of R,;. The effect of the input loss of the 
tube is not shown because the grid current was disregarded; however, this 
loss may be reduced to an equivalent R;. The resistance of the plate cir- 
cuit Re affects the frequency in a similar manner. The effects of these 
resistances on frequency are less for low values of plate circuit impedances. 

The required value of R, gives a measure of amplitude of oscillation 
because it is necessary for oscillations to build up until the internal plate 
resistance is equal to the calculated value. It is found that R, increases 
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Fig. 12.9—Experimental curves, showing the relation between the frequency and the 
resistance of the oscillatory circuit 
gradually to a maximum as the common frequency for the two types of 
circuits is approached then abruptly drops. 

Vigoureux” analyzes the crystal oscillator in a manner similar to 
Terry and correlates his interpretations of the equations with considerable 
experimental data, some of which are shown in Figs. 12.7, 12.8, 12.9 and 
12.10. He points out that there is an optimum value of grid capacitance 
with the crystal connected between grid and plate and a certain amount 
of grid-plate capacitance is required when the crystal is connected between 
grid and cathode. 

Wheeler" does not assume a linear static tube characteristic but 
represents it by a three-term nonlinear expression. The results are more 
complex and it is necessary in the end to disregard certain resistance terms. 
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Fig. 12.10—Experimental curves, showing the relation between the frequency of a quartz 
oscillator and the capacitance of the oscillatory circuit for various values of 
the grid leak 
12.20 SoLUTION BY COMPLEX FUNCTIONS 
The analysis of oscillator circuits may be simplified when only steady 
state conditions are of interest, all circuit elements are considered linear, 
and certain requirements which define the conditions necessary for oscilla- 
tions are known. Under these conditions the common circuit equations 
of complex numbers give the information desired. In this method the 
voltage induced in the plate circuit is considered the driving voltage which 
produces a current in the grid circuit (see Fig. 12.11). The network be- 
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Fig. 12.11—Equivalent circuit of Pierce and Miller types of oscillators shown 
in Fig. 12.12 
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tween plate and grid may be of any type and oscillations are maintained 
when the total gain through the circuit is unity (gain of tubes = attenuation 
through circuit) and the phase relation between the induced plate voltage 
(uV,) and the grid voltage (V7) is 180° (the phase shift is zero when 4 is 
considered negative). The expression u8 = 1 defines these requirements. 
Llewellyn” applies this method to oscillator circuits in general and Koga” 
uses it to study the crystal oscillator in particular. 

The equations are developed on the assumption that the grid-voltage vs. 
plate-current characteristic of the tube is linear. The fundamental equa- 
tion of w8 is given by the ratio of the voltage developed across the grid 
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Fig. 12.12—Circuit diagrams of crystal oscillators with crystal connected from grid 
to cathode (A) and grid to plate (B) 


circuit by the fictitious driving voltage uV, to the voltage V,. For the 
general circuit, Fig. 12.11, it is 
I2Z2 —yuZ,Z2 
B= — = —_ 12.20 
‘ Ve -RpZ, + 2i(Z2 + Zs) a 
where 
Z, - Z; + Z2 + Z3 


It is more convenient to write this in the reciprocal form 


Lo RZ + 21 (Zo + 2s) _ 
ups —plZ, Ze 

In applying this to the crystal oscillator, the additional assumptions made 

are that the grid current is negligible and the resistance in the plate im- 

pedance Z, is zero. 


1 (12.21) 


12.21 CrysTat GRID TO CATHODE 


With the assumptions made above and the crystal connected from grid 
to cathode of the tube according to Fig. 12.12A, the impedances are 


Z4,= 5X1 22> Reg tjXeq 23 = jXs 
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where R,, is the effective resistance and NX,, the effective reactance of the 
crystal, the grid resistance R, and the circuit capacitance C, in paralle) at 
the oscillating frequency. Upon substitution of these in (12.21) 

i - [Reg R, “2 X1(Xeg + X;)| + VX Rey + R,X:;) 

ug uX, Xeg a) puxXy Reg 


where 
Xe = Xi + Xeg + Xz 


1 
Thus — is of the form 


up 


which means that P = 1 and Q = 0. 
This results in the following two equations obtained from the real and 
imaginary parts of (12.22) both of which must be satisfied for oscillations 
to be maintained. 
The real part of (12.22) gives 
Xi(u + 1) (Reg + Xeg) + Xeg Xs 
Reg(X1 + Xs) 
and from the imaginary part is obtained 
X, X3 =, Rey R,, 
Ry deo 


—R, = 


X, = 


cg 


where $y = (This ratio of reactance to resistance of the « rystal circuit 


\eg 
will appear in various equations later.) 
Equation (12.24) may be said to detine the oscillating frequency and is in 


a convenient form to examine the effect of the various circuit variables upon 


the frequency. The impedances Y,, R.,, X,, and X3 may be thought of as 
forming an oscillating loop (See Fig. 12.11). For oxcillations to be main- 
tained in such a loop the sum of the reactances must equal zero and the 
sum of the resistances must equal zero. But the sum of the resistances 
cannot equal zero since R,, is the only resistance in the loop and it is posi- 
tive. It is therefore necessary for the driving voltage uV, to act upon the 
circuit and supply the energy dissipated by the resistance R., (and also R, 
through which the energy is supplied). This alters the frequency some- 
what and it is no longer determined by setting the three reactances equa! 
to zero as may be seen by equation (12.24). Nevertheless, the right side 
of this equation is small and approaches zero when R., approaches zero. 
It also becomes very small when the reactance VY; becomes small and R,., 
is not too great. This is the same condition as found by the differentia! 
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equation method and illustrated in Fig. 12.6 by the G-C curves. As the 
plate reactance XY, is made small the frequency increases and approaches a 
limiting value but does not quite reach it. This limiting value is the fre- 
quency at which XY, = 0. The dotted G-C curve shows that R., tends to 
lower the frequency and determines how close the limiting frequency is 
approached. The plate circuit resistance R. (component of 7Z;), if con- 
sidered, would have a similar effect as shown by the experimental curves 
12.9. The grid resistance R, (component of Z2) has an opposite effect as 
shown in Figure 12.10 because increasing R, is equivalent to decreasing 
the effective resistance R,,. 

The effect of the various constants of the crystal and circuit upon the 
oscillating frequency may be obtained from (12.24) upon substitution of 
these constants for the reactances and resistance R.,. The equation is put 
in a more convenient form for this purpose by Koga.” Equation (12.21) 
is written, 

1 1 m 


re, aes = 0 
, i? - R, 12. 
2 3 Za(1 + R,/Zi + Z et ( 
42 43 


It is assumed that the current in the grid branch is small compared to the 
plate current. This reduces the equation to 
1 1 Mm 


‘19 
Z2 7 Z3 7 Z3(1 + Rp Z)) (12. 


The admittance expression for the crystal is 


1 1 
R, - jel be wl’ oe w(Co rm = 


R I 1 l 2K r Gy as + C 
i+tjol; wa, ale +a 
(12.27) 


Note that Koga considers the air gap capacitance Cy as a separate factor 
but it may be included in the other constants of the crystal in which case 
the equivalent circuit is as shown in Fig. 12.3. With the crystal con- 
nected between grid and cathode the various circuit admittances are: 


Be, ete ss 
“7 Jules 


3 Joly 
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After substitution of these values of the admittances in (12.26) and setting 
the real and imaginary parts equal to zero, the following two equations are 
obtained: 


Ry 


lene (<“a) +t 
‘ 1 1 "\Co + Cs R 
slid Ez ~ wr w(Cyo + 55 | 


(12.29) 


CoC uC 


aa ett : 
spe 1+ oe - WC) 


Equation (12.28) gives the conditions necessary for oscillations and (12.29) 
gives the oscillating frequency as explained below: 
12.22 FREQUENCY OF OSCILLATIONS FOR G-C CONNECTION OF CRYSTAL 


Equation (12.29) for frequency is simplified by the fact that over the 
narrow frequency range considered, the reactances of Lz and C2 do not 
change appreciably. Also at the oscillating frequency, 


. 1 1 3 
Ri < [ots = — w(Co + at 
With these approximations (12.29) may be written 
1 4 1 1 a 
~ LC, © L1Cy Co+C: , Co (12.30) 
—— is 
Cy C4 
where 


’ ; uC: 
Cc, =C 9 + C3 + . ‘ 2 
1+ R(- - wvCs) 

wo Le 


and wy, is a constant approximating the oscillating frequency. 


FREQUENCY INCREASE IN CYCLES PER SECOND 
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Since the frequency is a function of the internal plate resistance of the 
tube (R,) and this is in turn a function of the other circuit variables, the 
frequency equation (12.30) is not sufficient to calculate the frequency. 
However, qualitative effects of the various circuit components upon fre- 
quency are obtained by assuming R, an independent variable. It is readily 
seen that an increase in R, increases the frequency. The effect of the air 
gap between crystal and electrodes, which is represented by the capacitance 
Cy, and the effect of the capacitance across the crystal C, are illustrated in 
Fig. (12.13).* To determine the frequency change caused by tuning of 
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FREQUENCY INCREASE IN CYCLES PER SECOND 


Fig. 12.13—Experimental curves, showing the effect of crystal air gap and grid 
capacitance on the frequency of oscillations 
the plate circuit (variations of C2) requires the calculation of the change of 
the variable part of C,. This quantity is 


ee uC; 
on si 3 F 12.31 
1 -t K; ( = woC:) \ 

wo Le 


The plot of C, is shown in Fig. (12.14A). The frequency decrease is pro- 
portional to the increase in C,. This is indicated in Fig. (12.14B).  Oscil- 
1 
lations stop before the point weC2 = 1. reached. The frequency thus 
Wyle 

varies in the same manner as shown in Fig. (12.6) but the curve is reversed 
because of the fact that the independent variable is taken as C2 instead of 
the frequency function of C2. 

The frequency change resulting from variations in the grid-plate capaci- 
tance C; depends also upon the value of C, as seen from (12.31). It is also 


* See also: “The Piezoelectric Resonator and the Effect of Electrode Spacing upon Fre- 
quency,” Walter G. Cady, Physics, Vol. 7, July 1936. 
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seen that the smaller the value of Cs (lower the plate reactance) the les: 
effect will the tube constants wu , R, and C; have upon the frequency. The 
circuit is therefore more stable. For this reason it has become customary 
to measure the frequency of crystals with the capacitance C2 reduced to ; 
value below that which gives maximum amplitude of oscillations. 


12.23 AMPLITUDE OF OSCILLATIONS FOR G-C CONNECTION OF CRYSTAI 


A measure of the amplitude of oscillations is obtained from (12.28 
which expresses the necessary conditions for oscillations to be maintained. 
In order for oscillations to start the expression must be negative, and, as 
the amplitude builds up, R, increases which reduces the negative terms 


Cy 


| 
| 
| 
Cs 
B 


Fig. 12.14—The variation of grid to cathode capacitance (A) and oscillator frequency 
(B) with change in plate circuit capacitance. Crystal connected grid to cathode 
until the equality is satisfied. The difference between the positive and 
negative terms is therefore a measure of the amplitude of oscillations. 

Equation (12.28) may be written 


1 
y-[w+Q|-4 (12.32) 


where A is a measure of the amplitude, 


1 
R, & * mast ao C :) 


¥ = wC3a» 1 2 (12.33) 
i+ Ki ( “ ws) 
Wo Ls 


and 
fCot AVP CoC. a ia uC; 
®) = Rw + . +C, + Cs + : 


y ; : ‘ 1 z 
C 4 ( 0 + € 4 1 + R’, (. F — os Gs) 
042 


(12.34) 


FUNCTION 





PIEZOELECTRIC CRYSTALS IN OSCILLATOR CIRCUITS 
where again Rj is assumed small compared to 
I 1 1 ‘ 
tc aaa Tica ~ n 
woC,  awo(Co + Cs) 
and wo is considered a constant. 


Equation (12.32) shows that in order to obtain a large amplitude y 
should be large and ®» should be small. With this in mind equations (12.33) 











FUNCTION 





* “Ca sails, 
Fig. 12.15—Functions from which the activity variations (A) are determined as the 
plate circuit capacitance is varied. Crystal connected grid to cathode 


and (12.34) may be analyzed to determine the relation between the circuit 
components and amplitude. It is found that for maximum amplitude 

C, and R; should be small, 

C, should be large, 

C’; has an optimum value, and 

R, should be large. 


1 1 
As to the plate circuit, the amplitude is maximum when =~ = ] 
Pp Wol.e2 


1 
aCe. A plot of Y and Bp + R. is shown in Fig. 12.15. The difference 
g 


between these two curves is a measure of the amplitude and is sho 


1 

curve A. Oscillations can exist only where y lies over ®y + R.” The shar 
gy 

ness of y varies considerably with the value of R, and the resistance of th 


L, — C2 circuit. The latter is disregarded for simplicity. Here again the 
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results can only be considered a first approximation, but agree with actua! 
conditions sufficiently to be of considerable interest. 


12.24 Crystat Grip TO PLATE 


The equation (12.20) is general and for the condition of crystal connected 
between grid and plate of the tube (See Figure 12.12B) Z; represents the 
crystal impedance which will be called Z, = R. + jX., also: Z; = jX,, 
Z. = jX2,and X, = X,+ X2+ X,. 

Note that R, and C; are disregarded in this case because their effects are 
similar to those determined for the foregoing case of crystal connected grid 
to cathode. 

After substitution of these values in (12.20) the real part is found to be 


= (u + 1)X1X2 7 X; Xe 


R, R 


(12.35) 


and the imaginary part is 


_ RX 


ae = 
Ry 


(12.36) 
which shows the effect of the various variables on the frequency. The right 
hand side of equation (12.36) is comparatively small and the frequency is 
therefore close to a value fo which makes X, = 0. In this case the frequency 
is above the limiting frequency fy because the right hand side is positive 
since X, is negative, whereas it was found that the frequency was below fp 
for the crystal connected between grid and cathode. As R, and X, are 
increased the frequency will increase and as Ry, increases the frequency 
decreases. These interpretations are verified by the G-P curves of Figures 
12.6, 12.9 and 12.10. 

The effects of the various circuit and crystal constants are determined 
by Koga” by writing the general uf equation as 


uZe 


Zs+ 22+ 7 + R,/Z,~ 


0 (12.37) 


After substitution for the Z’s, the real and imaginary parts are respectively, 


ol- = - 


wl} wl 


+ Ri . Be 
i) R 4 | 1 1 | 2 wl, 


(12.38) 
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Fig. 12.16—Frequency and activity change for variations in the plate circuit 
capacitance. Crystal connected grid to plate 
There are two values of w which satisfy (12.39) but only one of these w, 
will satisfy (12.38). At this value of w, 


. 1 i? 
Ri < | ots sage Sige 21 (12.40) 


By introduction of this and the assumption that wo is essentially constant, 
(12.38) may be written 
4 4 + 
1Lo], ~ . : 
, C C 2 
Se ee 
wo 


(12.41) 
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This is an approximation for the conditions for oscillation and relative 
amplitude. 
The frequency equation (12.39) becomes 


‘7% 4 4 
nn” as 


be 1 


G 2/_! .\ 
: 1 + R; 7~— oa Ce 
(On) ly 
and wy is a fixed value written in place of w,. Figure 12.16 shows the fre- 
quency and amplitude changes as a function of C2 for the crystal connected 
between grid and plate. 


where 
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Fig. 12.17—Generalized oscillator circuit in the form of a filter network 
12.25 CONDITION p8 = 1 FOR CIRCUITS IN GENERAL 


It is convenient to apply the rule 48 = 1 as the condition for sustained 
oscillations to more complex oscillator circuits. The circuits may be drawn 
as shown in Figure 12.17 and the characteristics of the filter network between 
transmitting and receiving end may be analyzed by conventional filter 
theory to determine the conditions which fulfill the oscillation requirements. 
An example of this is the oscillator shown in Figure 12.18A. The equiva- 
lent configuration, Figure 12.18B, indicates that the crystal is part of a low 
pass filter and the frequency of operation is that at which the total phase 
shift is 180°, 

Oscillators involving more than one tube may also be inspected in this 
manner. Figure 12.19 is a two tube oscillator designed to operate at a 
frequency close to the resonant frequency of the crystal. The proper phase 
shift is obtained by a two-stage amplifier and, therefore, no phase shift is 
required through the crystal network. The crystal thus must operate as a 
resistance which it can only do at its resonant or antiresonant frequency. 
Since the transmission through the crystal branch is very low at the anti- 
resonant frequency of the crystal, it will oscillate only at the resonant 
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° 13 . e . ° ° 
frequency. Heegner” explains a number of crystal oscillator circuits by 
the method briefly outlined above. 


{{]+ 
1 IT 








a 7o9 3 
: 








>Rp 
to, Tt 








| 
eee eee 











B 


Fig. 12.18—The oscillator circuit (A) is equivalent to the filter circuit (B 





























Fig. 12.19—Oscillator circuit in which the crystal operates at its series 
resonant frequency 


12.30 Vercror MetHop oF OSCILLATOR ANALYSIS 


A convenient method of examining the effect of certain circuit variables 
on frequency and the necessary conditions for oscillation is by the vector 
representation of the voltages and currents in the circuit. Much of Heis- 
ing’s early work on the analysis of electric oscillators by vector methods is 
directly applicable to crystal oscillators. Boella’’ analyzed the crystal 
oscillator circuit by this method and treated in detail the effect of the 
decrement of the crystal on the oscillating frequency. Since some engineers 
prefer this method of qualitative analysis to approximate equations it will 
be briefly explained. 

The vector diagrams for the two conditions, crystal between grid and 
plate and between grid and cathode, are shown in Figure 12.20A and B as 
applied to the circuit diagraras, Figure 12.12A and B, respectively when in 
the simplified form of Figure 12.11. The necessary conditions for oscilla- 
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tions are that V, is in phase with and equal to uV, (note that » is considered 
negative). Like Koga, Boella assumes the current J, small compared to /, , 
hence the voltage drop across Z, is approximately 7,/,. The angle this 
makes with V, is determined by the value of Z; and the internal plate 
impedance R,. Any change in either of these requires a change in the 
angles y and y’ in order that V, shall be in phase with u.V,. This means 
that the frequency must vary to produce this change in andy’. Because 
of the rapid change in the reactance and resistance of the crystal with 
frequency, these requirements are met with very little change in frequency, 
which accounts for the high degree of frequency stability obtained with 
crystals. This is described more in detail in a later section. 





B 


Fig. 12.20—Vector diagrams of currents and voltages in the oscillator circuit Figure 
12.11 with crystal connected grid to cathode (A) and grid to plate (B) 


12.31 CHANGE IN FREQUENCY WITH DECREMENT OF CRYSTAL 


It has been found that for the crystal connected from grid to cathode 
there is a maximum theoretical frequency at which the circuit can be made 
to oscillate by reducing the plate circuit impedance. This also corresponds 
to the minimum frequency which can be obtained with the crystal connected 
between grid and plate. This was called the limiting frequency fo. It is 
interesting to note that fp is determined by the intersection of the reactance 
curve of the crystal plotted as a function of frequency and the reactance 
curve of the capacitance in series with the crystal. This series capacitance 
is the grid-plate capacitance for one case and the grid-cathode capacitance 
for the other. As illustrated in the curves Figure 12.21, the limiting fre- 
quency fo increases as the decrement of the crystal increases. 

The difference between the true frequency of oscillations and fo increases 
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as the plate impedance is increased and as the losses in any of the circuit 
elements increase. This is necessary for the proper angle of Y + y’ in the 
vector diagram. With the G-P connections, the departure from fo and 
change in fy as the decrement of the quartz varies are in the same direction, 
while for the grid-cathode connection they vary in opposite directions, 
and the net result will depend upon the value of the internal plate resistance 
and plate circuit impedance. The curves of Figure 12.21 show that the 
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x 10° 0.6 
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f —- 
Fig. 12.21—The change in reactance characteristic of a crystal resulting from a change 
in decrement 
change in fy for a given change in decrement is less for smaller values of 
aC. (larger values of series capacitance C3). That is, the effect of the de- 
WL 3 
crement of the crystal upon the oscillating frequency is small when the crys- 


tal is operated near its frequency of resonance. 


12.40 NEGATIVE RESISTANCE METHOD OF ANALYSIS 


The methods of analyzing oscillator circuits described in the previous 
sections define the operation in terms of the individual circuit elements and 
the crystal is treated as one of the circuit elements. Certain advantages 
result, however, by grouping all the circuit elements, except the crystal, 
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Here Z, represents the 
impedance looking into the oscillator from the crystal terminals. 

The requirements for sustained oscillations are that the sum of the 
reactances around the loop equal zero and the sum of the resistances equa! 
zero as previously stated in section 12.21. These conditions are obtained 
when Z; is a negative resistance p in parallel with (or in series with) a 
capacitance C; as shown in Fig. 12.22C. The crystal is considered to be 
operating as an inductance L, and resistance R, as determined in the pre- 
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Fig. 12.22—-Equivalent representations of crystal and oscillator circuit 





























vious sections. The frequency equation has been derived by Reich” 
from the differential equation for the current in the loop. It is 


p+R 1 
p L.C, 
and the condition for oscillation is shown to be 
(12.44) 


We shall consider the crystal connected between the grid and cathode of 
the tube, in which case Z, is the input impedance of the vacuum tube. 


rr . . 17 . . . . 
rhe expression for g, Was developed by Chaffee from which it is possible 
+t 


to determine the circuit conditions necessary for the input resistance and 
reactance to be negative. The effect of the circuit variables upon the abso- 
lute values of p and C, determines their effect upon the frequency and activity 
according to equations (12.43) and (12.44). 
12.41 InpuT ADMITTANCE OF THE VACUUM TUBE 
With the assumption that the grid current is negligible and the static 
tube capacitances C, and Cy are part of the external circuit, Chaffee’s 
equation for input conductance becomes 
_ Cyw(K + Gi) + CzauK(Cyw — By) 


(K + Gi) + Ge — BP (12.45) 
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and for the input susceptance 
CywuK(K + Gi) — Cpa (C3x — By) 
b = —C3w — — \ Bed. . = ss = 2 (12.46) 
(K + G))- + (C38 — B,)? 


where A and yu are detined as follows: 


Ot, 
(e, constant) 


ey 


de, 3 
u=—-|- (?, constant) 


and G; and B, are the conductance and susceptance of the plate circuit. 
If we let 

Gy 

Bi 


wl "3 
= : and B 


K 


(12.45) becomes 
(1 + A#B) + u(t 


(1 + AB)? + a(t - s) 


g = C3wA 


and (12.46) becomes 


w(l + AB) —. ( 


(12.48) 
(1 + AB) + A’ 


r 
| 
b= — +. 


When the resistance of the plate circuit is neglected (i.e. = 0), and n> 1 
we may write 


u(A — B) 
(12.49 
1+ (A — By ) 


b —- BA -B : 
-= —4l® \ u (12.50) 
K i+ (4 — BP 
These equations are in a convenient form to determine the effect of the 
plate tuning f(B) and grid-plate capacitance f(.1) upon the resistance p 
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and capacitance C, with the assumptions of no grid current, low plate 
circuit resistance, and u > 1. 

From (12.49) it is seen that in order for g to be negative, B must be positive 
and greater than A, since A is normally positive. That is, the plate circuit 
reactance must be positive and less than the grid-plate reactance when the 
latter is a capacitance. Under these conditions 6/K and hence the input 
reactance will be negative according to (12.50). 


Curves of 6K are shown in Fig. 12.23 with B as independent variable 
and A as parameter. These curves indicate frequency change. On the 











Fig. 12.23—Variations in the input impedance functions of an oscillator circuit for 
changes in plate circuit tuning 


same figure is plotted 6, g called ¢,. This may be considered the sensitivity 
of the oscillator or, for a given value of w/,/R, of the crystal, it represents 
the activity. The similarity between these curves and the actual change in 
frequency and activity normally experienced is apparent. 

It should be pointed out here that the presence of harmonics is effective 
in changing the input impedance of the vacuum tube and hence the fre- 
quency and activity of the oscillator. The presence of harmonics results 
from the non-linear characteristics of the vacuum tube. Llewellyn” 
explains that a non-linear resistance may be represented by a linear re- 
sistance plus a linear reactance. From what has been said concerning the 
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frequency of the oscillating loop, it is apparent that this effective reactance 
will alter the frequency. However, this reactance is small when the im- 
pedance of the circuit is low at the harmonic frequencies and is zero when 
the external circuit is a pure resistance. 


12.50 EFFICIENCY AND POWER OUTPUT OF OSCILLATORS 


In many applications of crystal oscillators the efficiency and power output 
are important factors. These are not treated here but reference is made to 
the work of Heising™ which covers this aspect for various electric oscillator 
circuits. Much of the analysis is directly applicable to crystal oscillators. 

12.60 FREQUENCY STABILITY OF CRYSTAL OSCILLATORS 

The equations for frequency show that the frequency is governed some- 
what by the amplification factor, the grid resistance and internal plate 
resistance of the vacuum tube. Since these factors are functions of voltages 
applied to the tube and amplitude of oscillation, they cannot be considered 
fixed. If the frequency change resulting from these variables is great, the 
frequency stability is said to be low, and if very little frequency change 
takes place the frequency is determined principally by the circuit constants 
and the frequency stability is said to be high. 


re le : . 
Llewellyn” shows how it is possible to compensate for the change in plate 


resistance by the proper value of circuit elements. This was done by deter- 
mining the relations necessary for R, to be eliminated from the frequency 
equation. It is sometimes helpful in designing very stable oscillators for 
frequency standards to select circuit elements which will reduce the effect 
of plate voltage changes on the frequency. It is more the purpose of this 
section, however, to show Llewellyn’s derivation of the equations for fre- 
quency stability which have not heretofore been published and from them 
point out the characteristic of crystals which enable them to stabilize 
oscillators. 


12.61 THe FREQUENCY STABILITY EQUATION 
The steady state oscillating condition is 
up = 1 (12.51) 
In general 8 is a function of the frequency, the amplitude of oscillations, 
and of some independent variable V. This independent variable is the one 
for which it is desired to stabilize the frequency. It may be the potential 


applied to the tube, or it may be a capacitance located somewhere in the 
circuit. § depends upon these three variables thus: 


u8 = f(p,a, V) (12.52) 
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Instead of the frequency, a more general symbol p is used and may be 
thought of as the differential operator d/d¢ which occurs in the fundamen ta| 
linear differential equations taken as describing the oscillatory system. 
That is 

d 


lla Nar Ha 


The function uf may have the form 


up = Ae” 


The result of taking a general variation 6 of (12.54) is then 


= + i560 = 0 


Since (12.54) is a function of the three variables p, a, and V the variational 
equation (12.55) may be expressed in terms of partial derivatives with 
respect to these three variables. That is 
1] 0A OA On. .| 00 00 OF =o, 
— 6 —- da + ~~ 6V |] +7] — 66 + — ba + — 6V | =0 
AL Op e+ da al ap P da ay 
(12.56) 


The solution of (12.56) for the variation in p is 


.! ya ie 
1 ( = 6] aa ca) 
~: = & {6557} 


1 0A 4 ; 90 
A Op Op 
It is a property of functions of complex variables that, provided they 
possess derivatives at all, then the value of the derivative is the same regard- 
less of the direction in which the limiting point is approached. This fact is 
expressed by 
04 O41 — , dA] 
Op da . Ow 
00 00 j 00 


Op da Ow 





and 6p = da + tdw 


and provides means by which the real and imaginary parts of (12.57) may 
be separated to yield the two equations 


104/09... . 30 df1 eA... 9. 
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10A/10A.,, 1 0A 00/00... 00 
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The variable p in general may be written as the sum of @ and iw. With 


and 


bw 


the remembrance that # is the differential operator d/dt and that a set of 


linear equations expresses the transient condition, it is evident that the 


a+iw)t 


current will have the form Je” which is equivalent to Je Inspection 
of this shows that the real part of p, namely a, determines whether the cur- 
rents in the system are going to build up with time, or die away with time, 
or remain constant, depending respectively upon whether a is greater than 
zero, is less than zero, or is actually equal to zero. With this in mind we 
see that (12.59) and (12.60) state the change in @ and w respectively which 
would result from some change in the circuit condition. Initially the 
circuit was oscillating in a steady manner so that a was zero and w had some 
particular value. A change in V then occurred. This produced a change in 
the amplitude accompanied by a change in the frequency as expressed by 
(12.60) and a change in the transient terma. Suppose now that the change 
in V were very small. Then in order for oscillations again to assume a 
steady value it is necessary for the amplitude ‘‘a’’ to change a sufficient 
amount to cause a to become zero. Thus in (12.59) we put da equal to 
zero and solve for the required amplitude change. This may then be elimi- 
nated from (12.60) resulting in the final expression 


104 00 190A 06 
by <1 OV Ga A 0a Vg, (1261) 
104 00 1 AA 08 


A dw da A Oa Aw 
which gives the frequency change éw in terms of the change of the inde- 
pendent variable 6V. 
12.62 FREQUENCY STABILITY OF CONVENTIONAL OSCILLATOR 
In applying this equation to the oscillator circuit, Fig. 12.24, we must 
first set up the conditions for oscillations. The u@ equation is 
uX,X2R, 
WO ™ 3X,RpR, — Xi Xa Xa — (Rp Xa(X + Xs) + R, Xe + Xd 
(12.62) 
The oscillating conditions u8 = 1 requires 
X,R,R, = Xi:X2X3 
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and 

MXiX2R, + RpX2(X1 + Xs) + RgXi(X2+ X3) =O (12.63 
It will be assumed that the following relations exist: 


p= fil V), Re = fo(a), Ks = NX, + Xe ot Ns; => fs(w), R, = a constant 
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Fig. 12.24—Equivalent oscillator circuit analyzed for frequency stability 
Then we obtain from (12.63) 
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By substitution of these values in (12.61) and disregard of Y, in comparison 
with all other X’s the equation for frequency stability is obtained as 
ee 
— >, Xi X2Xs3 
dw pw Ot 


iV = . aX aX aX (12.65) 
t OXe 0 

a Aj As A3 | 

( +) ( Ow ” Ow + Ow ) Ry R, 


From this we learn that the values of the reactances X,, No, and NV, 
should be small and the values of R, and R, large to give small changes in 
w when V is varied. These variables are more or less limited, however, by 
the conditions necessary for sustained oscillations according to equation 
(12.63). It is important to notice that the denominator of (12.65) contains 
functions which do not appear in equation (12.63) and hence may be of any 
value. These factors are the rates of change of the various reactances with 
frequency. For given values of circuit constants, the equation shows that 
the frequency stability increases as these rates of change increase. 

12.63 FREQUENCY STABILITY COEFFICIENT OF CRYSTALS 

The rate of change of the reactance of an element is referred to as the 
“frequency stability coefficient”* of the element. Expressed in per cent, 
we have for the frequency stability coefficient of a reactance 


dX w 


a eas 


(12.66) 
Let us now examine the frequency stability coefficient of a crystal which 
is used as the reactance Ys when connected between grid and cathode of the 
tube and as X3 when connected between grid and plate (See Fig. 12.24). 
The resistance of the crystal will be assumed to equal zero due to the negli- 
gible effect of the resistance variations upon the reactance for crystals with 
average Q and operated at a frequency not too near the anti-resonant 
frequency. (This may be observed in Fig. 12.21.) 
The reactance of the crystal! then is 
2 


x 2% 
7, eo Oy 


X, = (12.67) 


= 9 ” 
wow — w 


where 


= 2x X frequency 
= 2m X resonant frequency 
= 27 X anti-resonant frequency 


* First suggested by N. E. Sowers. 
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By substitution of the relations 


9 
2 


be 


Cy ws —- w _ 
_ = — and — = Xp 
Co wi wl 


into (12.67) we obtained 


X.= Xo a | (12.68 
Cows — w 
and by differentiation 


* , . 2 ’ 24 
dX. _ -*/ 1  . | - xe] os | (12.69) 
dw w Cow, — w CoL(w — w)” 


w . “ye - . 
Multiply by <> to obtain the stability coefficient 
«le 


x 

1X Xo CG w | 
F(X,) = > = —]1 - 2. 

(X-) X,. dw S| Co ws — w" 


12.70) 
~ Xo GQ o wt ( 


X. Co wt (ak — wv)? 
and eliminate w by substituting in (12.70) the relations obtained from 
equation (12.68). These are 





d o(-He 7 
“o-e XoJ Ci 
d 
ws : (12.71) 
w Ci Ci 1 
ae . > 
wi Co + Co Xe 
1 a 
Ao_] 
Thus 
aa Sy wein4 OS . ote 
F(X.) = -1l ~is i= X Cy 1 ra X. (12.72) 


which may be written 


: , a 
F(X.) = -¥ E + (1 _ *) + a (1 -**) (12.73) 


c 





The stability coefficient of a coil and condenser F(X)’ may be obtained 
from (12.73) by letting C; = «. Then 


F(X)’ = + E + (1 —_ «| (12.74) 
<10 <he 








ee ee ee ea Sati ocla’ abakb 





The 
rati 


Th 


fun 


det 


10° 
ant 
coe 
eq 
fre 








PIEZOFLECTRIC CRYSTALS IN OSCILLATOR CIRCUITS 





197 


The comparative stability of the crystal and tuned circuit is given by the 














ratio 
“Spe X\ 
F(X, C X we 
( ewe te ee sf (12.75) 
7 F(X) C 1 Xo 
4 ous = Y. 
2.68) FF = 
: 10,000 
2.69) & 
& 
by 
4 
y 1000 
Si 
a 
Som 
- A Vin 
2 70) = x 100 
ufo 
; 
from & 
10 
71) 
71 ) i : . F r Al 
f 1 10 100 1000 10,000 
i Co 
/ Cy 
i Fig. 12.25—The stability coefficient of a crystal as compared to a coil and condenser 
i for variations of the ratio of capacitances 
rm. * . . . + - Xo . Co . 
) 79) rhis ratio is plotted in Fig. 12.25 for — = —1 and with C the independent 
fae <le 1 
. , le, , R 
function. It is apparent that the value of C of a crystal is the factor which 
1 
determines its frequency stability for given operating conditions. For an 
73) AT crystal in an air gap holder, the ratio of capacitances is of the order of 
3 ° eye 7s : 3 : 
10° and its stability coefficient is therefore 2.6 X 10° greater than for a simple 
: y g | 
ined anti-resonant circuit. Since this is so much greater than the stability 
coefficients for the other reactances which appear in the denominator of 
equation (12.65) it represents the order of magnitude of improvement of the 
a 
74) 





frequency stability of an oscillator obtained by the use of a crystal. 
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The fact that the frequency stability of a crystal oscillator is a function 
Co ‘ a rg 
of C explains why a BT cut crystal is in general more stable than an {7 

1 
cut. The two may be made equal, however, by adding capacitance across 
the AT cut crystal. 

Actually we have compared the frequency stability obtained by the use 
of one type of circuit (the equivalent crystal circuit) with one of a different 
configuration obtained by making C; = «. In practice this is usually the 
case since C; must be large to obtain oscillations when using coils and con- 


0 : 1) 

~ at which oscilla- 
1 

tions stop and this is determined by the Q of the circuit elements as shown 

in the next section which deals with activity. It will be shown that the () 


densers. The limiting factor is therefore the value of 


‘ . : Co ‘ 
required is proportional to C and therefore the maximum frequency sta- 
1 


bility that can be obtained is directly related to Q. 


12.70 RELATION BETWEEN CRYSTAL QUALITY AND AMPLITUDE OF 
OSCILLATIONS 


The activity of a crystal is usually thought of as the relative amount oi 
grid current produced in an oscillator circuit. This method of defining 
activity affords a means of comparing the quality of one crystal with another 
for a particular set of conditions. The disadvantages are first; it is only a 
relative measure, and second; it is not possible to compute the activity as 
thus defined by any method of oscillator analysis so far presented. Curves 
have been shown of amplitude of oscillations as a function of certain circuit 
variables, but these represent only qualitative changes associated with 
plate resistance variations. The first objection has been somewhat recti- 
fied by the use of reference oscillators* in which all the circuit elements 
including the tubes have been carefully matched. There is still the diffi- 
culty, however, of comparing crystals of different frequencies for it cannot 
be assumed that the measurements are independent of this variable. It 
would be more desirable to have some absolute measure of activity and 
particularly one which would lend itself to convenient computation from 
readily measurable constants of the crystal. 


12.71 DEFINITION OF CRYSTAL QUALITY FOR OSCILLATOR PURPOSES 


In deriving an expression for the quality of a crystal, it is convenient to 
use the negative resistance concept of the oscillator as described in section 
12.40. The equations are general and in a form which admit of separating 


* Developed by G. M. Thurston. 
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he crystal from the oscillator circuit. Equation (12.44) which gives the 
conditions necessary for oscillations to exist, may be written in the form 
- wl puta 
Le = —— (12.76) 
R. 
In order for oscillations to start, the right side of this equation must be 
equal to or greater than the left side. If it is greater, oscillations build up 
causing p to increase until the equality is satisfied. The difference between 
these two terms before oscillations start is therefore a relative measure of 
the final amplitude for a particular oscillator. The absolute value of am- 
plitude cannot be obtained from equation (12.76) since we do not know the 
relation between p and amplitude. However the greater the magnitude of 
the greater will be the amplitude of oscillations for a given set of oscil- 
c 
lator conditions. This term may therefore be considered a measure of 
crystal quality. It is the effective Q of the crystal unit as measured at its 
two terminals and at the operating frequency. To distinguish this from 
the Q of the crystal as usually spoken of, it will be called ¢. . 
In the same respect the left side of equation (12.76) may be thought of asa 
1 
measure of quality of the oscillator circuit, that is, pwC, = — , then (12.76) 
Yo 
becomes 
SL og = 1 
12.72 Ficure or Merit or Crystats —M 
It is very inconvenient to use ¢, as a figure of merit of the crystal because 
it is a complex function of the constants of the crystal circuits, Figure 12.3, 
and also the frequency. The computation of ¢, from such measurable 
characteristics as frequency of resonance f/,, frequency of anti-resonance 


fs, resonant resistance R; , and static capacity Co, requires considerable 


time and effort. 
It is highly desirable that a simple, easily determined expression for a 
figure of merit be found. The steps to indicate a suitable one are as follows: 
The equation for ¢, in terms of measurable quantities for computing it is 
derived from equation (12.27) and given by the formula 


whi (w — wi)(w — ws) 


wl, Rj w! 


—1] 
R. rs why ws — wi 
R, w 
By letting 
why ws =" 


R, aw" 
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2 
— ow 


| 
Ee 


n 


tote tore 


2 (12.80) 
~ os 

and with the assumption that over the narrow frequency range that the 
crystal will operate 


(12.81) 


w 
equation (12.78) is reduced to 


ee nM’ (n — 1) (1 
¥e M 


2.82) 


In this equation it will be observed are only two variables, namely, » 
which varies widely as the frequency is varied between f; and fy and M 
which is substantially constant over the same range. 

A set of curves is plotted in Fig. 12.26 for a hypothetical set of crystals 
having values of M of 1, 2, 5, and 10 with varied over a range that falls 
between measured frequencies f,; and fz. Studies will show that whenever 
M increases ¢, will increase. M is readily calculated from measured 
constants as seen from the following: From (12.79) 

M = wl Ww — wi 4 W1 Ly ws 5 wi (12.83) 


R, w?* R, WW 
With the assumption in (12.81) 


M a wl Cy =a 1 
Ri Co w1Co R, 
which is a simple expression containing three of the four measured quanti- 
ties mentioned above, and which bears a direct relation to activity for a 
given value of the frequency variable ». M is the new figure of merit of 
the crystal. 


(12.84) 


Figure 12.26 contains a further indication which is useful on occasions. 
Here ¢, is not positive at any frequency unless M is greater than 2. But ¢, 
must be positive for the crystal to oscillate in the two general types of cir- 
cuits considered here.* It provides a measurable index to separate com- 
pletely non-useful crystals from those that can oscillate in a given circuit. 

Equation (12.84) may be written 


M = 2 (12.85) 


* For a description of oscillator circuits which do not require the crystal to exhibit a 
yositive reactance see: “A New Direct Crystal-Controlled Oscillator for Ultra-Short- 
Nave Frequencies” by W. P. Mason and I. E. Fair, Proc. [.R.E., Vol. 30, p. 464, Oct 
1943. 
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where Q is the Q of the crystal and r its ratio of capacitances. Thus the 
figure of merit involves the dissipation in the crystal determined by Q 
and the piezo-electric effect determined by r.”° 

It was pointed out in the preceding section that the frequency stability 
increases as r is increased. The above equation shows that Q must increase 
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Fig. 12.26—The dependence of the quality function ¢. of a crystal upon frequency 
and figure of merit M 
proportionately if the same figure of merit is to be maintained. The fre- 
quency stability obtainable in a particular oscillator is therefore limited 
by the Q of the crystal and the frequency stability coefficients should be 
compared on this basis. 


12.80 Activiry OF CRYSTALS 


In deriving a figure of merit for crystals as oscillators, it was found that 
the amplitude of oscillations in a given circuit not only depends upon M 
but also it is a function of frequency relative to the resonant frequency 
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of the crystal. This may be explained by referring to Fig. 12.27 which 
shows curves of the reactance Y, of the crystal plotted as a function of fre- 
quency. The frequency at which oscillations occur depends principally 
upon the value of circuit capacitance C;. Equation (12.43) shows that 
the frequency must adjust itself to a value at which C; resonates with 
the reactance of the crystal. This value of reactance is represented on the 
curve as X,. and the corresponding frequency of oscillations as f,. ‘The 


























f2 
Fig. 12.27—Two crystals having the same figure of merit but with different reactance 
characteristics Y. and XY; will operate with different amplitudes according to the 
relative values of ¢eo and gj respectively 


circuit capacitance (or Y,.) has been so chosen in the illustration that /, 
lies equidistant between the resonant frequency /; and the anti-resonant 
frequency fe. 


Then at this value of f, , g. is a maximum, as shown by the 
¢e curve, and for a given value of ¢, this will result in the greatest activity. 
Now let the capacitance Cy of the crystal be increased. The frequencies /; 
and f2 will then become closer and at the same time the height of the ¢. 
curve is reduced. Assume also that the Q is increased in order to maintain 
the same value of M and hence the same maximum value for g,. The 
reactance-frequency curve for the modified crystal and the corresponding 
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curve for g, are represented by the dotted curves. Note that the oscillating 
frequency fo for the new curve is closer to f; than it is for f; , therefore gi, is 
less than ¢- and the amplitude of oscillations will be less. Thus two 
crystals may have the same value of M but will not give the same output 
unless operated at the same relative frequency with respect to their resonant 
and anti-resonant frequencies. It would not have been possible to increase 
the oscillator output by increasing C; so as to lower the frequency f, because 
by doing so ¢y is decreased more rapidly than g. is increased, and the result 
would be a further decrease of activity. It is therefore necessary in deriving 
an expression for activity to include the variable of relative frequency or 
which we have shown to be a function of the reactance of the circuit and 
the crystal constants. 


12.81 DERIVATION OF PERFORMANCE INDEX OF CRySTALS—PI 


It will be assumed in the first derivation that the negative resistance p of 
the circuit is much greater than the effective resistance of the crystal R. 
under stable oscillating conditions. Equation (12.43) which expresses the 
frequency then becomes 


= | 
= 12.86 
w ‘7. G ( ) 


This leads to a very simple solution from which a more exact expression is 
later obtained. 
Equation (12.44), which expresses conditions necessary for oscillations, 
may be written 
| wh, | * 
ip| Ss aC, R, (12.87) 
As before, the numerical difference between p and the right side of the 
equation is a measure of activity. In fact, the right side of the equation 
may be considered to be an expression of the activity performance provided 
the terms are themselves independent of p. This is not quite true, since 
previous sections show that the capacitance C; is not entirely independent of 
the activity. (See equations (12.30) and (12.46).) However, this effect 
may be considered negligible for most practical purposes and the value of 
the right side of (12.87) called the Performance Index (PJ) of the 
crystal. From equations (12.86) and (12.87) the performance index is found 
to be 


1 
R. w C t 
This equation may be greatly in error under operating condition which makes 
R. large compared to p. Also R, varies rapidly with frequency and is 


PI = (12.88) 
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difficult to evaluate. R, is most readily eliminated from the equation by 
revising the picture slightly. With reference to the simplified oscillator 
circuit, Fig. 12.22B, it is apparent that the static crystal capacitance Cy and 
the circuit capacitance C, may be combined. This leaves for the crystal 
branch the inductance L; (different from L.) which is a function of frequency 
and the resonant resistance of the crystal R; which is not a function of 
frequency. Now R, may be assumed small compared to p with considerable 
accuracy. It is only necessary, then, to replace C; in equation (12.88) with 
(Cy + C,) and R, by R;. This equation then becomes 
1 


” Re Cs + OF 


An exact equation for PI is derived in section 12.83 and it is shown that the 
error in the simple expression above will in most cases be very small. 

An approximate equation for the relation between R; and R, is obtained 
by dividing (12.88) by (12.89). We thus find 


Ri (Co + Ci)? 
R.Ci 
or the effective resistance of the crystal at the operating frequency is 


2 
R. = R; (@ + 1) (12.91) 
C; 


Because of the approximation in equation (12.88) the equation for R, above 


PI (12.89 


(12.90) 


‘ Co $ 2 ; , ; 
is accurate only when C. +1) <« M*as will be shown in section 12.83. 
t 


The expression for PJ as given by (12.89) may be written 
1 


ror ae 29 
ri ie (1 + =) (12.92 
‘ Co 


which is the most convenient form for calculating PJ from the constants 
of the crystal and the oscillator circuit. 


12.82 RELATION BETWEEN M anp PI 


It was found that 


1 
= 2.93 
M ath (1 ) 


over the narrow frequency range 


. . . 1 
and this is essentially equal to oC Ri 


considered. Therefore, 


(12.94) 
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This gives a relation between the performance index and the figure of merit 
of the crystal. 
Another useful relation between M and PI is obtained from (12.89) and 
(12.93). Equation (12.93) may be written 
Xo 
M = — 12.95) 
Ri ( 
This is of the same form as the Q of a coil when Xp is considered to be the 
reactance of the coil and R, its resistance. Like the Q of a coil M is essen- 
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Fig. 12.28—The change in P/ and oscillating frequency of a crystal as the shunt 
capacitance is increased 
tially constant over a wide frequency range. Now if we let C; approach 
zero in (12.89) it becomes 


 . 
Rw C R; 


PI = (12.96) 
This equation for P/ is of the same form as the anti-resonant impedance of a 
coil and condenser in parallel, and like this impedance it changes rapidiy 
with frequency. The maximum value of PJ is therefore Xo times M and is 
obtained when C, = 0. Figure 12.28 shows a curve of ©] PI plotted as a 


, ae as ; 
function of —. This curve represents the change in activity as capacitance 
0 
is added across the crystals (increase in C;). 
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12.83 Exact EXPRESSIONS FOR PI AND R, 


The error in PJ caused by the assumption that the frequency is independ- 
ent of the crystal resistance Rj, that is, by use of approximate equation 
(12.86) for the frequency, may be investigated as follows: 

The impedance of the crystal and C; in parallel is given by 


1 


a. es ra i ale 
Z= wlCe re C)a i m? P*) [P 7 +. mP (m 1))] (12.97 
where 
_ _MCo _%-@ 
Cot C: a ~ a 


w3 = 2m times frequency of anti-resonance of the crystal and C; combina- 
tion when R; = 0 

w, = 2m times frequency of resonance of the crystal and C; combination 
when R; = 0 

w = 2m times operating frequency 


(Note that P is the figure of merit of the crystal and C; in parallel.) The 
imaginary part of Z is 


1 + mP*(m — 1) 


X= ~ Co + CQ PP) 


(12.98 


The condition for stable oscillations requires Y = 0. For this condition 


eg 
=i a 2.99) 
i, Vv; P? ( 


which defines the exact frequency of oscillation. The negative sign before 
the radical is used since the effective resistance is greater at this frequency, 
thus requiring less negative conductance for oscillation. 

With P large (m — 0) the frequency of oscillations coincides with w; 
which is the same as given by the approximate frequency equation (12.86). 
The real part of (12.97) is 


P 
R= (Q+Col +m P) (12.100) 
and when m = 0 
Rn ne a 2 (12.101) 
w(Co + C2) wo (1 + oi) 
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This is identical to the expression for PJ of equation (12.94). PI is therefore 
the anti-resonant resistance of the crystal and capacitance C; in parallel. 
Substitution of the value of m as given by (12.99) into (12.100) gives the 
anti-resonant resistance at the oscillating frequency. Thus 


_ 
a. l / ~ Pp (12.102) 


oe = | 1 - — - 
“ w(Co + Cr) 2 
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Fig. 12.29—The error in PJ resulting from the use of approximate equation (12.94) 


which is the exact expression for PJ, The differential error resulting from 
the use of approximate equation (12.94) is then 


i 
nok te 
R, < 4 
r+ 4/1-3 


The per cent error as a function of P is shown in Fig. (12.29). The error 
diminishes rapidly with increase in P and is negligible for crystals that are of 
sufficient quality for most oscillator purposes. 

Equation (12.91) for R, is also approximate because of the assumption 
that the frequency is independent of R.. A more exact expression will be 
derived. 

The impedance of the crystal alone is 

1 
~ @Cy(1 + 2? M’) 


(12.103) 


Z. 


[M —j(l + M?(n — 1))] (12.104) 
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the figure of merit of the crystal 


2 


2 
We —- Ww 


w — wi 
w, = 2m times frequency of resonance of the crystal 
we = 2m times frequency of anti-resonance of the crystal 
w the independent variable, 27 times frequency 
Co = the static capacitance of the crystal 


The resistive component of Z, is 

; M R 

wCo(1 + n?M*)~ 1 A 
WE +n 


~~ (12.105) 


In order to express R, in terms of Cy and C;, , the quantity 7 must be expressed 
in terms of these variables which define the oscillating frequency. This is 
accomplished as follows: 

The equation of ratio of capacitances of a crystal is 
Ci va Ws oie ar 


; 2 (12.106) 
Co w} 
Similarly, when the capacitance C; is placed across the crystal 


9 


C) w3 — WI ~ 

= ;*#—, (12.107) 
Cot CG W] 

where ws; is 27 times the anti-resonant frequency of the crystal and C; ii 

parallel. The ratio of (12.107) and (12.106) is 

Co ws = ‘wl 


, ’ = ” ” 
Co + Ci Wo. — Ww] 


(12.108) 


The oscillating frequency is given by (12.99) in which m is as defined 
under (12.97). The oscillating frequency w is therefore given by 


4 
adh y/ ~ Pp (12.109) 


2? 


r 
wi w | shied y/' - Pp (12.110) 
= 2 


2 2 
3 ~~ WI 
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Che angular frequency ws is eliminated by multiplying this by (12.108). 


Thus 
w} oe a —Co shi V'-> (12.111) 
w> —_ w} ~ & + C: 2 
or 
4 
a wat gal _ oily 1+ / — 44 (12.112) 
w—w Cot, . 


This is the value for 7 at the oscillating frequency and may be reduced to 


the form, 
i en 1+ 4/1-3 4 
(12.113) 
2 Co 
| (C41) | 


When this value of 7 is substituted in (12.105), the value of R.is found to 


be 
R. = 
1 E a ~ Pr: /\- | (12.114) 
TE + : 
- « Co | 
| (Os) 


4 : 
For crystals of usable quality Pp? << land by this assumption the equation 


n= 


reduces to 


Ri 
1 1 


+ ; (12.115) 
uM: (4 4 1) 


i : Co 7 
his again reduces to (12.91) when M° >> (¢ + 1). 


t 


R, = 


12.84 FREQUENCY CHANGE RESULTING FROM PARALLELING CAPACITANCE 

It is often desirable to know how much the frequency of an oscillator may 
be changed by varying the capacitance C, across the crystal. This is de- 
termined from (12.112) which gives the oscillating frequency as a function 
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4 
of C,. For practical considerations we may assume Pp 1 which reduces 
the equation to 
.. Soa * (12.11 
oy & 2 116) 
Ws — @] Co “bk C; 4 
From this and (12.106) we obtain 
2 2 
— CiC e 
——— “= aa bat (12.117) 
wij Co (Co + C:) 
Since 
w—w — (uw — w) (we +6) _ 2(w — w) 
Ww? wi isin wy 
then 
wW — we —1 
~~ , \ (12.118) 
Ww} 2(' 0 + i) 
C; 


where r is the ratio of the capacitances of the crystal. 
A curve of per cent frequency change multiplied by r as a function of 


’ 


Ce. a3 ea ‘ : ; 
cs shown on Fig. 12.28 for comparison with the associated PJ change. 
0 


12.85 RELATION BETWEEN PI AND OSCILLATOR ACTIVITY 

The relation between PJ and activity obtained in a particular oscillator 
will now be examined. Let the curves of Fig. 12.30 represent the variations 
of p with amplitude for two oscillators A and B, or they might be for the 
same oscillator at widely different frequencies. These are characteristics of 
the. oscillator circuits and may be of any shape. However, for oscillators 
with grid leak bias, the curves normally have no negative slopes. The rate 
of change of p depends upon the rate of change of u and plate resistance of 
the vacuum tube as shown by (12.45) for input conductance.* Since p 
builds up to a value equal to PJ we may plot PJ for p. The grid current 
T, is usually taken as a measure of amplitude. Therefore, Fig. 12.30 may 
be plotted as shown in Fig. 12.31 where PJ is the independent variable. 
These curves are the characteristics of the oscillator circuits A and B with 7] 
detining the quality of the crystal when used with a particular value of C,. 
It is characteristic of oscillators to ‘saturate’ as shown by the curves. 


* It is also a function of grid resistance but this does not appear in the approximate 
equation (12.45) because of the assumption of no grid current. See Chaffee’s complete 
equation for input admittance. 
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Fig. 12.30—Hypothetical curves illustrating the normal relation between the negative 
resistance of oscillator circuits and the amplitude of oscillations 
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Fig. 12.31—By interchanging the coordinates of Figure 12.30 the curves will represent 
the relation between PJ and oscillator grid current 
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Some oscillators saturate very rapidly and completely according to curve B 
and no further output is obtained regardless of the improvement in the 
crystal quality. For this reason it has not been possible in the past to 
separate the performance of the oscillator and the crystal since both were 
based upon the grid current as a measure of quality. By defining crystal 
activities and oscillator sensitivity in the manner outlined, the crystal and 
circuit can be studied separately. The per cent of crystals obtainable with 
PT above a certain value will be known and the design and improvement of 
oscillator circuits will be facilitated. 
12.86 Usk or PI ix Crystat DESIGN 

The expression of P/ in terms of the crystal constants and C; as given by 
equations (12.89) or (12.92) assists in the design of crystals. As an ex- 
ample, the effect of changing the area of the crystal electrodes will be com- 
puted. The Q of a crystal is defined as 


1 
) = : 12.119) 
( w C 1Ry ( I 
By introduction of the ratio of capacitances of the crystal r = a equation 
1 
(12.119) becomes 
) f 12.120 
e= wi Co Ri iceman 
or 
QO 1 
= : = M 12.121) 
r wi Co Ri ( 


Assuming Q and r do not vary, that is, disregarding effects such as secondary 
modes, change in damping produced by the mounting etc., and substituting 
(12.121) in (12.94) we obtain 

Q Co 


OF . > 
ail r w(Co + C1)” 


(12.122) 


- ; ; a 
where v is considered constant. Differentiating (12.122) with respect to 
r 


Co we find that PJ is a maximum when Cy = C,;. Since Co is proportional 
to the area of the electrodes this establishes the optimum area for a par- 
ticular value of circuit capacitance. 

The capacitance of BT-cut plates is 1.68 mmf. per square centimeter pcr 
megacycle.* Substitution of this for Cy in (12.122) gives 


6 
py = 268 X 10°MA 
(1.68 Af + C,? 


* All frequencies are referred to the time interval of one second throughout this paper: 
i.e. megacycles per second is called simply megacycles as is customary in the radio field. 


(12.123) 
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where A = the area of the crystal in square centimeters. 
f = the frequency in megacycles 
C, = circuit capacitance in mmf. 
M = figure of merit of the crystal (assumed constant) 


Thus for crystals of a given area, the performance index should decrease as 
the frequency increases. Figure (12.32) shows the theoretical variations of 
PI as the function of the diameter of the electrodes of three frequencies and 


10 
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Fig. 12.32—Theoretical curves showing the relations of PJ, electrode diameter, and 
crystal frequency for BT crystals and a circuit capacitance of 50 wuf 


for a circuit capacitance of 50 mmf. The activity of a 4-megacycle crystal 
with 11-mm. diameter electrodes is about the same as a 10-megacycle 
crystal with 18 mm. electrodes. It must be remembered in making this 
comparison that it is assumed that the damping introduced by the mounting 
is the same in both cases. Actually the damping is much greater for low- 
frequency crystals of this type than for high-frequency ones and maximum 
PI occurs at some intermediate frequency as shown by the curves of Fig. 
12.33. These curves show that the damping caused by the particular 
mounting used was small for frequencies above 6 megacycles but increases 
rapidly below this value. 
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Fig. 12.33——Calculated and measured values of PI for BT crystals. The discrepancy 
is a measure of mounting loss 


12.87 MEASUREMENT OF PI AND M 


In all the discussions so far regarding the performance of crystals in 
oscillator circuits, the crystal has been represented by the equivalent circuit 
of Fig. 12.3 in which all the elements were considered constant. It is 
possible to obtain crystals in which this is essentially the case, but in genera! 
there are three secondary effects which complicate the picture. These are, 
first, the effect of other modes of vibration of the crystal, second, variations 
in the crystal constants resulting from variations in the amplitude of vibra- 
tion, and third, the leakage or dielectric loss in the crystal holder. These 
factors will be considered in the order named. 

Secondary modes of vibration affect the crystal for oscillator purposes 
only when the frequencies of these modes are sufficiently close to the princi- 
pal one to alter its impedance characteristic in the frequency range of 
oscillation; that is, to alter the reactance as shown in Fig. 12.27 between the 
frequency f; and fy and the corresponding effective resistance between these 
two frequencies. With interfering’ modes present, the equivalent crystal 
circuit is so complicated as to make it impractical to compute PJ or M from 
such measurable quantities as resonant resistance R;, series resonant 
frequency /;, anti-resonant frequency f2, etc. For this reason it is 
necessary to measure the reactance and effective resistance of the crystal at 
the operating frequency in order to obtain a measure of crystal quality 
which will correlate with the crystal performance. For the same reason 
it is important when comparing oscillator circuits that the crystal should be 
operated at the same frequency in each case. 

It is believed that the non-linear effect noticed in crystals when used as 
oscillators is produced by the changes in the mounting as the amplitude of 
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vibration varies. The PI of some clamped or pressure-mounted crystals 
has been found to vary as much as 50% with change in drive. Noticeable 
frequency change also occurs. A change in the nature of secondary modes 
as the amplitude is varied has also been observed. Some secondary modes 
which interfere with large amplitude of vibrations practically disappear 
when the amplitude is reduced. This may be explained by the fact that 
certain modes are damped out by the pressure of the mounting and with 
large amplitude of vibration the effect of the pressure is reduced. 

The dielectric loss of the holder was considered negligible in the theory 
but it is found that certain phenolic holders have equivalent high-fre- 
quency leakage resistances less than 100,000 ohms. This resistance is in 
parallel with the crystal and will therefore reduce the PJ according to the 
equation 


PI,Ri 
PI = 12.12 
PI. + Ri 7 
where 
PI = resulting P/ 
PI, = calculated PT 
R, = equivalent high-frequency leakage resistance 


Because of these secondary effects which are not negligible it is essential 
in measuring crystal activity that ;he frequency and voltage across the 
crystal be known. Standard test circuits should simulate operating condi- 
tions in this respect. With these considerations, a crystal PJ meter has 
been developed in which the frequency and amplitude may be adjusted to 
correlate with various oscillators. The principle of operation and perform- 
ance of this meter is described by C. W. Harrison.* 
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The Measurement of the Performance Index of Quartz Plates 


By C. W. HARRISON 
15.00 INTRODUCTION 


HE theory of the general behavior of crystals in oscillator circuits has 

been described by I. E. Fair’. In Fair’s paper as well as in others’, it 
has been pointed out that in the neighborhood of the operating frequency a 
crystal is equivalent to the circuit shown in Fig. 15.1A. The crystal 
possesses two resonant frequencies, a series resonant frequency determined 
by the effective inductance, Z; , and effective capacitance, C; , and an anti- 
resonant frequency determined by these same elements plus the paralleling 
capacitance, Cy. This paralleling capacitance is the static capacitance 
between electrodes of the crystal and any capacitance connected thereto by 
the crystal holder and lead wires within the holder. The dotted resistor, 
R,, shunting the equivalent crystal circuit represents the effective shunt 
loss of the holder. In the ideal case and in many practica! instances this 
loss is negligible. 

It is rather difficult to express the circuital merit of a crystal quantitatively 
ina single term such as has been fou ad useful for inductancesand capacitances. 
It is customary to express the circuital merit of these two e.ements in the 
form of the ratio of reactance to resistance. That is, for an inductance 


Q = “ (15.1) 
and for a capacitance 
i ie, (15.2) 
‘ wR 


For filter purposes, the Q of a crystal involving only the inductance, 1, and 
resistance, R, , of Fig. 15.1A is adequate to express its usefulness in certain 
parts of a filter network, but for oscillator purposes it is insufficient. At 
frequencies other than the series resonant frequency the paralleling capacitor 
Co together with the associated shunt loss of the holder enters into the 
determination of a crystal’s performance. The term ( therefore is not com- 
pletely indicative of the crystal performance. There has been devised, as 
pointed out in Fair’s paper’, a term called “figure of merit” for a crystal 
1T. E. Fair, “Piezoelectric Crystals in Oscillator Circuits,” this issue of the B.S.T.J. 


2K. S. Van Dyke, “The Electrical Network of a Piezo-Electric Resonator”, Physical 
Review, Vol. 25, pp. 895, 1925. 
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which involves all the elements in the effective crystal circuit, and this term 
is much more expressive of the quality of a crystal. The figure of merit is: 


why C; ne Q 


M = 


” 


where ‘‘r 


capacitance, Figure of merit is useful for expressing the quality of a crysi 
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Vig. 15.1-—-Electrical equivalent circuits of a piezoelectric crystal—(A) At any fre- 
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Fig. 15.2—Generalized oscillator circuit of the Pierce or Miller type. 


in its holder or mount; however by definition it is independent of the value of 
RK, , and does not permit a ready evaluation of the performance of the crystal 
in an oscillator of the type that may be represented by Fig. 15.2. Any 
oscillator that operates the crystal in the positive region of the reactance vs. 
frequency characteristic exhibits capacitive reactance and negative 
resistance paralleled across the terminals to which the crystal is connected. 
The operation of the crystal when connected to an oscillator will be influenced 
by the magnitude of these two terms, and the combination must operate at 
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such a frequency that the total reactance is zero and at such an amplitude 
that the total resistance is zero. The performance of the crystal will 
therefore not depend solely upon its figure of merit, but will involve the 
impedance of the remainder of the oscillator. Up to the present time circuit 
design engineers have not devised standards or units to express the quality 
of their oscillator circuits without the crystal, so there are no corresponding 
circuital units of quality with which to correlate figures of merit of crystals 
to ascertain the suitability of one for the other. 

It was a practice for many years for manufacturers to test crystals in a 
model of the oscillator in which the crystal was to be used. This required 
manufacturers to keep on hand models of all oscillators for which they 
expected to make crystals. To avoid the mounting number of such test 
oscillators a special test set was developed which could be adjusted to simu- 
late any oscillator. By correlating various oscillator circuits to a set of 
adjustments on the test set, the actual model of the oscillator can be dis- 
pensed with. This special test set usually referred to asthe ‘‘‘D” spec. test 
set’, eliminated the “‘file” of oscillators, and substituted a file of adjustment 
readings that would be theirequivalent. However, the ‘“D” spec. test set’ is 
still inadequate to the development engineer since it defines ‘‘activity”’ in 
terms of oscillator grid current rather than in terms of the electrical equiva- 
lent circuit of the crystal. The activity as expressed by grid current is a 
purely arbitrary standard and serves only as a means of determining the 
relative activity as against other crystals of the same frequency operated 
under the same circuit conditions. 

The need for a system of measurement using units that are fundamental 
and not empirical has led to the proposal of ‘‘Performance Index”. An 
instrument to make such measurements is to be described in this paper. 

Specifically the Performance Index is 


_ ol 
~ wR 
where C; is the paralleling capacitance that is found in the oscillator circuit 
to which the crystal is attached, and L and R represent the effective induc- 
tance and resistance of the crystal as measured at the operating frequency 
indicated in Fig. 15.1C which is its equivalent at that frequency. If the loss 
in the holder is so low that the resistance, R,, may be neglected, then P/ 
may be expressed in other relations that are more useful such as, 
M ) 


PI = ee 
Coli + °t) | 
& o( “ i 2) (1 


or PI = PR, | 


PI (15.4) 


on 
on 
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where the symbols R; and C» are as shown in Figs. 15.1A and 15.2 and /’ js 
expressed as 


M 
eo 
— 1 4. © (15.6) 


With the effective capacitance, C;, of the remainder of the oscillator added 
to the paralleling capacitance, Co, in Fig. 15.2, the operating frequency will be 
that frequency at which the combination will exhibit a pure resistance at 
the terminals 4B (excluding the generator ““X” which is involved in the 
measuring technique). This leads to the definition: 

The Performance Index is the anti-resonant resistance of the crystal and 


holder having in parallel with it the capacitance introduced by the remainder of 


the oscillator. 

The Performance Index is therefore a term to express performance not in 
terms of the grid current of some particular oscillator, but in fundamental 
circuital units—impedance. The Performance Index is a term that may be 
used to compare performance of crystals at different frequencies. Its 
value is independent of plate voltage, grid leak resistance, or of plate 
impedance. It provides a measuring stick that should replace the “activity” 
figures of grid current in so far as the crystal is concerned. It paves the way 
for the oscillator circuit designers to come forth with standards of measure- 
ment for the oscillator circuit without the crystal in the hope that the 
two may be quantitatively associated and lend themselves to theoretical 
calculation of full oscillator performance. 


15.10 THrorY oF MEASUREMENT 


The problems of measurement are most readily explained by reference to 
Fig. 15.2. The crystal provides elements 1;,C;, Riand Cy. The circuit of 
the oscillator provides an effective capacitance, C,, which is composed of grid 
and lead wire capacitances plus capacitance introduced from the plate cir- 
cuit. The frequency at which this combination exhibits anti-resonance as 
measured at AB is the oscillating frequency. The resistance when added to 
negative resistance, p, will be zero. Oscillations will start with p numeri- 
cally smaller than the anti-resonant resistance measured at AB, but the 
amplitude of oscillations will increase causing p to increase until p and Z,, 
are equal numerically. The primary problem is to measure the anti-resonant 
resistance at AB at the anti-resonant frequency with p disconnected. 

Measurement of anti-resonant resistance directly is very difficult. The 
current flowing into an anti-resonant circuit is too small to measure with 
the usual meters. Other devices for measuring the current are likely to 
introduce paralleling capacitance that will vitiate the readings. The sug- 
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vested method of measurement utilizes a suitable driving voltage at “X”’ 
(Fig. 15.2) and a means to indicate the voltage at ““X”’ as well as at points 
1B. From these and other measured constants, the anti-resonant im- 
pedance can be computed. 

This method of measurement has its own difficulties, but it is believed 
corrections can be made to allow for errors introduced. Fundamentally, 
the series resonant frequency and the anti-resonant frequency are the same 
only when the resistances in the inductive and capacitive branches are equal. 
When the resistance is practically all in the inductive branch, which is true 
in this case, the impedance between terminals 14, at the series resonant 
frequency will exhibit capacitive reactance, though the total impedance will 
scarcely be different from that at the anti-resonant frequency. In the 
Performance Index meter, although the voltage is introduced in series with 
the circuit, the frequency is adjusted to the point of maximum voltage across 
AB, which further minimizes this frequency difference. A second error is 
inherently introduced by the loss in the crystal holder. This means 
that the series resonant frequency is also altered by the presence of this loss. 
Errors of any seriousness will result from the assumption that the series 
resonant and anti-resonant frequencies are identical only when the resistance 
in the inductive branch and loss in the crystal holder approach the effective 
crystal reactance in magnitude. These errors will be discussed in greater 
detail in a succeeding section. 

The development and operation of a satisfactory meter to measure P/ 
(Performance Index) depends upon a number of factors such as: 

1. A method to determine capacitance, C;, of the circuit (Fig. 15.2). 

2. A generator ‘‘X”’ to produce the driving voltage e; having variability 
in frequency and negligible internal impedance. 

A current indicator that introduces a minimum of reactance and 

resistance. 

4. A circuit or method to indicate P/ directly, or with a minimum of 
calculations. 

5. A number of other factors associated with the above which will be 
mentioned at the logical times. 

To construct a measuring circuit to determine the anti-resonant imped- 
ance by means of a series circuit so as to avoid any unnecessary measure- 
ments and computations involves the following basic principle. Excluding 
p, Fig. 15.2 is essentially equivalent to the circuit used in Q meters. The 
ratio of voltage e, to the driving voltage e; is the voltage stepup or the Q of 
that part of the circuit containing the resistance. In this case, the resistance 
is in the crystal which at the operating frequency has an effective Q of 


w 


Qi = - (15.7) 
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where L and K represent the effective values of the crystal. Equation (15.4 
will be found to embody the relation 


PI = (;: (15.8 


where X, is the reactance of C,, the capacitance introduced by the circuit at 
the operating frequency. If e; is kept constant, e, will at all times be 
proportional to Q;. By insertion of an attenuator network, whose attenua- 
tion varies with frequency in the same manner as does the reactance of (,, 
between terminals 46 and the voltmeter, the meter indication will be 
proportional to the product of these quantities or proportional to PI. With 
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Fig. 15.3—Measuring circuits of the performance index meter. 








suitable calibrations, therefore, it should be possible to get indications of PJ 
as readily as is now done for Q. 

The circuit shown in Fig. 15.2 is now best redrawn as in Fig. 15.3A. The 
crystal embodying elements Z;, C;, Ri, and Cy of Fig. 15.2 is now repre- 
sented in Fig. 15.3A by the dotted.rectangle and as having effective induct- 
ance, L, and effective resistance, R, both of which are functions of frequency. 
Capacitance, C;, is simulated by capacitors C, plus C, and C; in series where 
C, represents the capacitance of the crystal socket. Zero internal impe- 
dance of the generator is simulated by maintaining the driving voltage con- 
stant at all times and at all frequencies. To facilitate explanation, the 
measured voltage e; at the place shown is considered to be the driving volt- 
age from a zero internal impedance generator. 

Instead of using an ammeter to indicate current in the circuit, a voltmeter 
is utilized to measure voltage across an element under such conditions as 
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not to introduce disturbing capacitance. Splitting the series capacitance 
into two parts, C, and C;, the latter fixed and large compared to C,, provides 
the impedance element across which the voltmeter is connected. The input 


capacitance of the voltmeter is incorporated in the magnitude of (;. A 


capacitance attenuator, A, of known or calibrated values interposed on the 
input of the voltmeter enables the voltmeter to be used to indicate voltage 
ratios in terms of the attenuator calibration. 

The measuring voltmeter and a shunting capacitance, C, , are connected in 
the plate circuit of the amplifier tube, V-1. This circuit provides sufficient 
gain to furnish an output voltage of measurable magnitude and also provides 
an output voltage inversely proportional to frequency. The indication of 
the output voltage is proportional to PJ. 

The utilization of a vacuum tube in a circuit leading to a quantitative 
measuring instrument such as the voltmeter across C, involves determination 
of tube constants or calibration. The determination of these constants is 
best evaluated experimentally. A calibrating circuit for that purpose is 
shown in Fig. 15.3B. A capacitance, Cy, of high impedance in series with 
comparatively negligible resistance, R., is connected across the driving volt- 
age terminals of e; with a voltmeter measuring e; giving a reading e,.. The 
second subscript ‘‘c’’ indicates calibration conditions. By connecting the 
input circuit of V-1 across this resistance, the attenuation variation with 
frequency of the Ry — Cy network cancels the attenuation variation with 
frequency in the plate circuit of V-1. The ratio of e;, to e,. will then be in- 
dependent of frequency. In the “calibrate” circuit (Fig. 15.3B), the capac- 
tior attenuator, -1,, interposed in the grid circuit is set at unity (minimum 
insertion loss) for a given deflection of the meter indicating e,. In the oper- 
ate circuit (Fig. 15.3A), the attenuator is readjusted so that voltage ep» 
produces the reading of e, as obtained in the calibrate position. The quan- 
titative action of the amplifier then may be expressed in terms of C,, R,, 
C, and a reading from the attenuator 1, as will be shown later, and it is 
constant and independent of frequency. By placing this resulting constant 
in an equation, which will also be derived later, the value of PJ may be de- 
termined in terms of such constant, of the reading of attenuator .1, and of a 
reading on the scale of C, that has been calibrated in terms of C). 

To facilitate still further the operation of the PJ meter, the voltage e; is 
produced as shown in Fig. 15.4 by arranging for the oscillator to have its 
frequency controlled by the crystal through feedback from capacitor, C;. . 
Automatic volume control is provided such that the amplitude of e; is 
essentially constant at all times and at all frequencies. The circuit is con- 
structed to oscillate at the desired frequency, and adjustment for insuring 
this operation is provided in the form of a phase shifting circuit with variable 
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capacitor, C,. After a crystal has been inserted in its proper place, osci |a- 
tions will begin, but may be slightly above or below the resonant frequency 
of the crystal plus C;. By adjustment of C, the frequency can be shifted 
the slight amount necessary for resonance. This is observed by placing 
switch S in the P/ position and making the adjustment to give maximum 
deflection of e,. 


SWITCH "S" 








ATTENUATOR 
LIMITING LINEAR 
AMP PWR AMP — PI AMP 












































| hig RECT. 


AND AMP 


Fig. 15.4-—Diagram of Performance Index meter. 








15.20 Derivation oF PI Crrcurr EQuATION 


The following circuit relations derived from Figure 15.3 show first, that 
the ratio of e, to e, is a function of the Performance Index of the crystal, 
and second, that the calibration circuit permits an absolute evaluation of 
its magnitude. 

At resonance, the effective circuit (, designated as (2 , is determined from 


Os mn €o + & - a) (1 +. a) (15.9) 


ej e 8 


Since the circuit Q includes the capacitance of C, as a part of the crystal, it is 
necessary to express (2 in terms of the crystal’s properties (see Fig. 15.1). 


Since Q; ( = *) of the crystal is independent of C,, the relationship between 


Q, and Qzis readily obtained by equating the expressions for the anti-resonant 
impedance first, when C, is considered to be in shunt with the series capac- 
itor, C,, and second, when C, is considered as part of the crystal. This re- 
sults in 

mii 
C; 


QO. = Ci 


(15.10) 


where 


enabli 


Expre 


hence 
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where 
. Ci, - 
C, = C. + ———, 15.11) 
il He So ( 
enabling Q; to be expressed as 
ey C.C. - 
), = : = (15.12) 
“= 2G — 
Expressing é in terms of e, , we have 
Me, | = ty Vv "7 oo Xe, (15.13) 
| ? (15.14) 
| Up ie r Je 
Xe, 
hence 
Cy Wy ae 
p= ai Xe, (15.15) 
Ge ly 
and 
3 =e, (15.16) 


With the above equations substituted in (15.12), we may express Q; as 


é CC CrCr Xo, 5.17 
a es — wae 
? ej ” Gn, (C, come Cur / T | lp | : | 


Now 
) 
Pt DS aR, (15.18) 
wl; 
Therefore 
en Cc Cy Pe c 
PI = P A - : . rhe | 15.19) 
ej Gn, (C, — Cz)* "A +|¥| 
If 


5.20) 


7 
+ | 

ere, 

rel 
Il? 
E 


PI = ep { Cre Cy 


4 15.21 
e; Gn, (C, it Gar ( ) 


The simplified PJ expression (15.21) assumes that the reactance of C, 
is small compared to the plate resistance and plate load resistance of V-1. 
The evaluation of PJ from this expression has three obvious difficulties: 
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(1) C,,/Gm is a quantity that is difficult to evaluate numerically, (2) ihe 
magnitude of PI is measured in terms of the ratio of the two voltages ¢ 
and e,, and (3) the measurement is dependent upon the gain of a vacuym 
tube amplifier, V-1. These difficulties may be materially reduced in their 
consequence by an internal calibration circuit. 

The internal calibration circuit (Fig. 15.3B) consists of a capacitor, ( 
and resistor, Ry, in series. If the reactance of C, is very much greater than 
R,, and the plate resistance of V-1 is very much greater than the reactance 
of C, , the calibration is essentially independent of frequency. 

The internal calibration circuit (Figure 15.3B) enables the evaluation of 
CpG, to be carried out. The additional subscript, c, indicates “calibrate” 
conditions, 


Cre 


V Ri + XCq 


Cic wC i Ry 


14 et 


Cove 
— = eg. 
a * 

Equation (15.15) remains the same for both “operate” and “calibrate’ 
conditions with the exception of the second subscript reserved for the 
“calibrate” operation. Therefore, by solving for C,,/G,» we find 

Cy| _ ge 


bate y/ + Ea (15.25) 
'p 


Equation (15.25) may be rewritten as (15.26), if (15.23) and (15.24) are sub- 
stituted in (15.25) 


G. _ Cic Rx C4 1 1 


Gy, 7 wake xin P R 2 15.26) 
¢, 1+[%o] /i+(2] (15.26 
'p Xe, 


If (15.26) is substituted in (15.19), it is found that PJ may be expressed 
as follows: 
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[he above equation involves only the original approximation that maxi- 
mum current indicates resonance. If R, and C, are selected such that 
R, << Ne, and if A, equals unity, then 


PI = jes] oa age, (15.28) 
€5 Cpe (C, — C,)° 
From this expression it can be seen that the ?J measurement is inde- 
pendent of calibration of both the e, and e; vacuum tube voltmeters, pro- 
vided that the same voltmeter scale factors are used for the ‘“‘operate’’ and 
“calibrate” conditions. The absolute calibration then depends on the 
magnitude of 1, Ri ,Ca,C,,C,andC,. The “multiply-by” factor that 
R,C4C, 
Cc, — C,)? 


Accurate evaluation of this quantity by capacitance and resistance meas- 


is to appear on the C, dial is determined by the magnitude of 


urements is a little difficult since the denominator represents the square of 
the difference of two small capacitances. When C; is large, the evaluation of 
this factor is helped considerably. This ‘‘multiply-by” factor may be 
experimentally determined by a voltage measuring means which permits 
an evaluation of this factor to a higher degree of accuracy. Substituting 
(15.11) in (15.28) we have 


PI = | 4 A = (1 + “) (15.29) 
pe k 8 


Now by shorting the crystal socket terminals (Fig. 15.3A) and applying a 
voltage e, at the e, generator terminals of external origin (the crystal oscil- 
lator circuit itself may be used if self-excitation is provided), ‘the current i, 
through the capacitors C;, and C, is given as 


oe €,aC,C; - €; wl A ; 
C.+C (1 + 7 (15.30) 
c, 
Now the voltage, e2 , across the series capacitor, Cy , is 
ty f 
1 (15.31) 


The ratio of e;/e: may be expressed as given in (15.32) when (15.31) is sub- 
stituted in (15.30) 


= (4 + o) (15.32) 
2 C 8 = 
If (15.32) is substituted in (15.29), we find 
lp Sic Ra C4 ey 


A 


PI 


15.3 
C. €o ( , 3) 


Cpe ej 
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- a : : 
lhe quantity — is readily determined by the attenuator, A, when the 
e2 


switch, S, (Fig. 15.4) is operated between ‘““M” and “PI” for the above 
described conditions. The absolute calibration then depends upon A, 2,. 


Cy and C;. All four of these quantities may be determined within a few 
per cent. 


15.30 OscILLATOR CORRELATION 


The equivalent crystal circuit has been discussed in so far as the measure- 
ment of PJ is concerned; however, for correlation with an oscillator, the 
behavior of the crystal in that oscillator must be duplicated. Correlation 
of the PJ meter with an oscillator is a function of both amplitude and 
frequency. It is obviously necessary from the derivation of (15.28) that 
the frequency of operation be duplicated, but the necessity for ampli- 
tude correlation can only be explained from the practical consideration that 
the equivalent circuit components of Fig. 15.1 are parameters that may be 
functions of amplitude. Crystals having nonlinear characteristics of the 
type that necessitate amplitude correlation may in part be attributed to 
either the method of mounting the crystal or couplings to other modes of 
vibration whose coupling coefficients are functions of amplitude. 

In most oscillators the voltage across the terminals of a crystal is a func- 
tion of many parameters such as plate voltage, vacuum tubes, etc. With 
an average set of conditions, however, reasonable correlation is obtained 
with the PJ meter for a single adjustment of the generator voltage, e; , for 
all crystals. The magnitude of e; must, of course, be chosen to produce a 
voltage across the crystal equal to the average value obtained in the oscil- 
lator circuit for which the crystal is intended. 

Frequency correlation with an external oscillator is a function of the 
effective capacitance, C;, in shunt with the crystal. In order to duplicate 
the oscillator frequency with the PJ meter, the capacitance, C,, (Fig. 15.4) 
must be adjusted until the frequency of oscillation in the PJ meter is the 
same as that in the oscillator for a crystal having average activity. In 
Fig. 15.4, the capacitance, C,, is variable, and its dial is calibrated in terms 
of both the total effective capacitance across the crystal, C;, and the resulting 

R, C4 Ci 
(Cy — Cz)? 
means of a capacitance bridge connected across the crystal socket terminals 
with the generator impedance shorted. 

The determination of the dynamic or effective capacitance, C;, across the 
crystal for an oscillator may similarly be obtained by adjusting the magni- 
tude of C, in the PJ meter until the frequencies of oscillation in the P/ 


multiplying factor The magnitude of C; may be measured by 








meter | 
of oscil 
oscillat 
inasn 


The 
be rep 
quency 
across 
This it 
by the 
teristic 
The o: 

1. | 

2. ] 


If t 
is ess 
maint 
mum 
2nn Vv 
result 
gener 
simils 
The | 
an ad 
In th 


conti 
from 
of A 
must 
in tl 
abso 
all ti 

TI 
































On 
nd 


li- 


Le RIA ASME IL IES LOLI SING EO ETRY HN itt tpn oan 





PERFORMANCE INDEX OF QUARTZ PLATES 229 


meter and in the oscillator under test are identical for the same amplitude 
of oscillation. By this means, the P/ meter directly indicates the effective 
oscillator capacitance, C,. The amplitude of oscillation must be duplicated 
in as much as C;, is not independent of the amplitude in most oscillators. 


15.40 DerscriIpTION OF OSCILLATOR GENERATING “e,” 

The generator plays no part in the theory of PJ measurement as it could 
be replaced by a signal generator or any other suitable source of radio fre- 
quency energy. It is convenient, however, to utilize the voltage appearing 
across C; as an input to an amplifier whose output represents the generator. 
This in effect constitutes a feedback oscillator whose frequency is controlled 
by the crystal under test. Initial consideration of the over-all charac- 
teristics of the PJ meter oscillator leads to the following requirements. 
The oscillator must, 

1. Be capable of oscillating all crystals usable in other oscillator circuits. 

2. Be capable of operating the crystal over a wide range of shunting 

capacitances in order to duplicate all the frequencies of oscillators 

now in the field. 

3. Be capable of permitting high degrees of AVC control in order to 
maintain the generator voltage constant while the frequency is ad- 
justed for reasonance. 

If the generator voltage, ¢, , is constant, resonance of the crystal circuit 
is essentially indicated by maximum crystal current, and oscillation is 
maintained at that resonant frequency. The adjustment to obtain maxi- 
mum current is such that the phase shift throughout the oscillator loop is 
2xn where n = 0, 1, 2,3, etc. As previously described the phase shift and 
resulting frequency of oscillation are varied by a tuned circuit. The 
generator voltage, e; , is held constant by an automatic amplitude control 
similar to the automatic volume control which is often applied to amplifiers. 
The manual control of the magnitude of the generator, e; , is provided by 
an adjustment of the bias voltage of the automatic amplitude control circuit. 
In this way the maximum or start gain is independent of the setting of the 
amplitude control. 

Automatic amplitude control (commonly referred to as automatic volume 
control, 4 VC) of an oscillator may be applied by the separation of the limiter 
from the linear amplifier. ‘This means that in order to apply a high degree 
of AVC to the PI oscillator (Fig. 15.4), the input voltage of the limiter 
must be above the threshold of limiting by an amount exceeding the variation 
in the 8 path caused by the AVC control. This enables the limiter to 
absorb the changes in the gain of the linear amplifier such that p38 = 1 at 
all times. 

The time constant of the limiter is fast compared to that of the AVC 
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circuit, a condition which permits damping of transients set up by chanzes 
in gain occurring from AVC action. The input of the linear amplitier 
is held constant by the limiter. Gain changes in the linear amplifier pro- 
duced by the variation of C, (Fig. 15.4) are absorbed by AVC, while the 
variation of activity in the crystal is absorbed by the limiting amplifier. 


15.50 EVALUATING PERFORMANCE INDEX 


From (15.28) it can be seen that the attenuator, A, the effective variable 
capacitance, C;, together with the vacuum tube voltmeters, e; and e,, 
offer a number of possible variations in the method of evaluating the con- 
stants used to determine the PJ of quartz crystals. There are, however, 
two principal methods—the first provides direct reading, while the second is 
more accurate but requires an indirect evaluation. 

The first method utilizes a means of calibration of the meter scales directly 
in terms of P7. The attenuator is adjusted such that its indicator reading 
times the multiplying factor associated with the dial attached to C, is some 
multiple of 10. If in the calibrate position, 4, is set at unity, and ep and 
éic are adjusted by varying the capacitive load to some reference deflection, 
then the expression for PJ becomes 

ep 


where 


= Cic Cy RaCa 
. (C; sian es] 


The Performance Index then is indicated by the two readings of e, and e;. 
The absolute magnitudes of e; and e, need not be known since it is possible 
to use as a reference, the arbitrary calibrating deflections of e;, and e€,.. 
The magnitude of e, indicates the significant figures while e, is a multiplying 
factor. 
The second method of evaluating Performance Index eliminates any 
calibration errors in the two vacuum tube voltmeters, e; and e,. This 
nace 2 Cp Cie . ” 
method utilizes the attenuator to adjust : . ] = 1. In the “calibrate” 


iope 

operation, ej, is set to equal e; and then e,, is adjusted for full scale or a 
convenient deflection. In the ‘‘operate’’ position, the attenuator is varied 
until e, equals e,-. In this manner, the two readings of the attenuator 


. . ep *  * 
are used to determine the ratio of ig and the measurement is independent 
pe 
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of the voltmeter calibration. The factor R«C, is a constant; therefore, the 
P.I. equation (15.28) simplifies to 


P.I, = (AA)Ke (1 


mn 
w 
wa 


where 


AA = change in attenuator insertion loss between the “operate” and 


> 


ti : 
“calibrate” conditions in terms of output voltage ratio ~~, given 





Cpe 
A 
as. 
« A. 
: CyRa Ca es Cae ; Bon 
Ky = (C Cc.) where C; is the effective capacitance in series with the 
’ ane zx 


crystal, C;, is the fixed series capacitance and C, is the crystal 
socket capacitance. (See Fig. 15.3) 


15.60 P.I. Meter AppLicaTIONS 


The application of this instrument can be extended to determine other 
properties of both crystal and oscillator. With the aid of a frequency 
measuring means and a capacitance bridge, the P.I. meter may be used to 
determine all the circuit constants designated in the electrical equivalent 
circuit of Fig. 15.1. If the loss in the holder is negligible then the equations 
are considerably simplified; however, in those instances where holder loss 
must be considered, the approximation that Y,;<< R, which may be allowed 
for most cases enables an evaluation of M and Q, that is readily computed. 

The dial controlling C,, that is calibrated in terms of the total capacitance, 
C,, makes possible the calculation of the magnitude of the input impedance 
to the crystal circuit, R, as well as Q1. 


xe 
R= 
PI. | 
, (15.36) 
PI. 
,= | 
V1 x, | 


The magnitude of Q; may also be measured directly from equation (15.12) 
where (; was given as 
a) CeCe e 
O; Se om . es (15.12) 
ey (C, —_ ee 
As may be seen from Fig. 15.4, eo/e; can be evaluated in terms of the 
attenuator calibration by enabling switch, S, to select the “PJ” and “M” 
positions respectively and adjusting the attenuator such that the same 
output meter indication is obtained in the two cases. 
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The quantity Q, makes possible the calculation of Q. where Q. is defined 
as 


O be wl, — Xi 
R. R, (1 r ) (15.37 
t 


It can be shown that (Q. in terms of Q; is given by the following equation if 
R,> X;. 
QO  @Q 

X, Pd. 15.38) 
1-a 1- ate 

R, R, 

The Figure of Merit, M, defined by (15.3) at the series resonant frequency 

of the crystal, a, becomes 


0. = 


1 Ma - 
M = =? 5.39) 
J a wy Co Ry R, (1 


The measurement of M can be determined from ( provided C; and C 
are known. M may be determined from the following expression 


1 C, Co 4 
=~ = PT. Co dnd C; (15.40) 
R, 


If C; is selected such that C, > Cy, then for most cases R, is large com- 
pared toQ,X,. This means that Q;=Q,and R= R,. With this approxi- 


mation, (15.40) becomes 
Cr =) - 
M=0,—11 — 15.41) 
Oral TC, ( 


The relationship between R, and R, as a function of frequency may be 
expressed directly from the input impedance expression of the equivalent 


circuit. This is given as 
R. = ° > - 
w—-w fl 1 (15.42) 
m, 2| + We 
We. — W ae 


If 


QO! ‘ : 
*The relationship R, = R; ( + 2) was derived in Fair’s paper*!. 
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which is a plausible assumption when C, > Cy, and we assume that 


1) + We 
pe , = w, then, 


R, 


saci k ~ a (15.43) 
WwW W@W) 


The relationship between R, and R, does not involve R, and may be 
expressed in terms of Cy and C; instead of frequency. Neglecting R, to 
determine the relationship between capacitance and frequency in the P.I. 
meter as derived in Section 15.92, we find 


oer [i+ 4/1 +4] 
re : (15.44) 


W2 — Wi Co 


where, 


P, = ’ ~ se 
Co ‘5 
(: % ‘) (15.45) 

C; 


If P; > 2 then the expression between R, and R; may be written 
Co : fac 
R. = Ri\ i+ C (15.46) 
t 


The restriction that C, > Cy may be removed if the error between (15.42) 
and (15.46) is taken into account. This error may be expressed as 


1 2 
Per Cent Error in R. = 100 |p| E akan al (15.47) 


The resonant resistance, R,, together with (15.46) provides a means of 
checking the P.I. meter. The magnitude of R; may be determined by the 
substitution method and from this, the value of R, calculated. Fig. 15.5 
represents the agreement between the P.I. meter and those expected from 
resonant frequency measurements. 

The remaining crystal constants 1; and C, (Fig. 15.1) may be evaluated 
from the measurement of w;, w. and Cy. The resonant frequency, w; , 
is defined as 


. 1 wor 
i= 15.48) 
w] LG ( 


The anti-resonant frequency, w., is defined as 


is 1 1 1 . 
; = 15.49 
— o- (4 + “) cai 
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Fig. 15.5—Typical characteristic of a quartz crystal measured by the Performance Index 
meter. 
Solving these equations simultaneously, it is found that 
1 1 ) 
Lh = 773 2 = >. Cf 
Co(we ™ 1) ~wWe Co(we = 1) 4 
‘ P (15.50) 
Ces Co(w2 — wi) L 2Co(we — w1) | 
1=- ‘ al 


wy Ww} ) 


15.70 EXPERIMENTAL DATA 


The performance of the P.I. meter may best be illustrated by experi- 
mental data. The following data indicate the correlation which may be 
obtained between the P.I. meter and various types of oscillator circuits. 


Experimental considerations are extended to 


(1) Frequency and amplitude correlation with a “Pierce” and “Tuned- 
Plate” oscillator 
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(2) The measurement of the effective capacitance, C;, of an oscillator as a 
function of tuning, and 
(3) The variation of P.I. as a function of voltage across the crystal. 
The results presented are not to be considered as generalized data, but are 
intended only to show a set of measurements obtained for a specific set of 
operating conditions for each type of circuit. 


It has been pointed out that the frequency of oscillation is a function of 
R,, #1, #2, Co and C;. Since R;, w1, w2 and Co are explicit parameters 
of the crystal, the capacitor, C;, becomes the only frequency determining 
element in the P.I. meter. From analytical methods to be described in 
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Fig. 15.6—Frequency variations in oscillators as a function of Performance Index. 


Sections 15.80 and 15.92 the frequency difference between the P.I. meter 
and the generalized oscillator (neglecting R,) is given as 


by 2 (we — 1) (1 4 os 


2 
, ) 
M C; [1+ yf + 4 | 
Pi (15.51) 


(@e — w) Ga 
we (1 t 3) 


Since the magnitude of frequency change in the above equation is small 
compared to the variations caused by changes in operating conditions, the 
P.I. meter may be used as a frequency correlation medium. Fig. 15.6 is 
an example of the correlation between the “Pierce” and ‘“Tuned-Plate” 
oscillator and the P.I. meter. The P.I. meter falls between these two 
oscillators in frequency for any crystal activity. It must be recognized 
that Fig. 15.6 is not conclusive to the extent of generalization; however, it is 
indicative of possible correlation with these two popular oscillator circuits. 
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Figures 15.7 (A) and (B) show the amplitude correlation between the 
Performance Index meter and the grid current of the ‘‘Tuned-Plate”’ and 
“Pierce” oscillators, respectively. The P.I. vs. grid current characteristi: 
was arbitrarily taken at three frequencies—4.5 mcs, 5.66 mcs and 7.81 mes. 
The change in grid current of the oscillator shown in Fig. 15.7 (A) with 
frequency is caused by the varying L-C ratio in the plate circuit. (See 
curves A, B and C for constant crystal activity.) The curves A and D 
represent the effect of changing bands by switching coils, varying the 1-( 
ratio 2 to 1 in the plate circuit for the same crystal frequency. Fig. 15.7 
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Fig. 15.7—-Typical oscillator characteristics. 


(B) represents the correlation between P.I. and grid current of the “Pierce” 
type oscillator. The results of this correlation indicate that the grid cur- 
rent is essentially independent of the operating frequency. 

The measurement of P.I. is independent of the level of crystal vibration, 
provided that the electrical equivalent circuit parameters of Fig. 15.1 
become constants; however, in actual practice these are not constants, 
particularly R,. Variations of this type, as previously discussed, make 
it necessary to duplicate the amplitude of oscillation of the P.I. meter with 
the oscillator. Fig. 15.8 represents the variation of P.I. as a function of 
voltage across the crystal terminals for five crystals arbitrarily selected. 
It is readily observed that P.I. may be a random function of amplitude. 
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Fig. 15.8—Observed variation of Performance Index as a function of voltage across the 
crystal terminals. 
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Fig. 15.9—Observed variation of effective crystal capacitance in a Miller oscillator. 
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Equation (15.44) indicates that the capacitance, C;, determines the oper- 
ating frequency between w; and a» of any given crystal. The capacitance. 
however, may include reflected reactances from associated circuits or pos. 
sibly from circuits unintentionally coupled to the oscillator. Normally, 
crystals are adjusted to frequency for a specified value of C;. This makes 
it of interest to measure the magnitude of C; over the range of the manua 
oscillator tuning adjustments, as well as over the frequency range of t! 
oscillator. Fig. 15.9 shows the circuit capacitance, C,, plotted as a functio 
of tuning of the plate circuit of a Tuned-Plate oscillator. Tuning of the 
plate circuit is expressed in terms of percentage of maximum grid current 


( 


15.80 Crrcurr ANALYsIS INVOLVING THE ACCURACY OF P.I. 
MEASUREMENTS 

The method of P.I. measurement just described involved a number of 
unverified approximations. These approximations under the majority of 
conditions will be proven to be justified and the resulting expressions for 
percentage error will be obtained. 

It is necessary to apply a method of analysis that is most readily adaptable 
to the crystal circuit. The analysis described in this section involves the 
use of Conformal representation as a means of determining (1) the behavior 
of the equivalent crystal circuit, (2) the error resulting in P.I. from assum- 
ing operation at the resonant frequency rather than the frequency for mini- 
mum impedance, and (3) the comparison of frequency of oscillation in the 
P.I. meter with other oscillator circuits. 

Generally, the variations of reactance, Y, and resistance, R, of the equiv- 
alent circuit (Fig. 15.1) are plotted as a function of frequency, and the 
analysis of the impedance, R + jX, between the resonant frequency, «, 
and anti-resonant frequency, w., are handled in precisely the same way as 
any linear passive element. This was essentially the procedure used to 
dezive the equations in section 15.60. The analysis required to evaluate the 
errors leads to rather an elaborate study; however, in Section 15.81, it will 
be shown that it is very helpful analytically if the impedance of the crystal! 
is plotted in the form of a circle diagram, that is, with the ordinate repre- 
senting reactance, X, and the abscissa representing resistance, R. 


15.81 Conformal Representation 


Conformal Representation or Mapping is a convenient tool which for this 
application enables the physical operating condition to be expressed quanti- 
tatively from its graphical counterpart. Physical interpretation also 
makes it possible to draw many other conclusions that by other methods 
prove clumsy and laborious. 
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The basis for this analysis depends upon the ability to utilize the following 
equation to represent any impedance whose frequency of operation is con- 
trolled by a crystal. This equation is known as a linear fractional trans- 
formation 

Wat + 84 

y + 6Z 

The terms a, 3, y, and 6 represent complex constants and Z represents a 
complex variable later to be chosen to represent a linear function of fre- 
quency. Since W and Z represent the dependent and independent variable, 
they may also be considered as representing two separate planes. The 
abscissa and ordinate of these two planes represent their real and imaginary 
components respectively. The planes are linked by (15.52), that is, this 
equation will transform a specific point from one plane to the other. 


The constants a, 8, y and 6 for the equivalent crystal circuit are determined 
by writing the expression for Z, in the form of (15.52). For example 
(neglecting R,), Z, from Fig. 15.1 may be written as 


1 
__ [ei(a-2)] 
Z lemme (15.53) 


” jul j; “ye 
Le sl («1 pa (<. “3 ))I 


4¢ 


By substituting (15.48) and (15.49) in (15.53), this impedance may be writ- 
ten as 
, 1 [wRi + jlilw — wi)(w + «)] 


‘fo = = ; (15.54) 
jwCo [wRi + jlilw — w2)(w + w2)] 


Since the operating frequency, w, represents some frequency between 
and w2, and w > wz — w:, we can make the following approximations 
in this operating range. The symbol w, is defined as the average operating 
radian frequency. 


_ a1 + ow 
2 


wa 


We a OF } 
If (15.55) is substituted in (15.54), factor w, out, add and subtract w from the 
imaginary component in both numeratorand denominator, we may write Z, as 


fern | : Lak. 
< oi o—-wW . (w — we) 
2 a Ee +d iil o» | - E <,| 2 pecadllas 


ll + li = 


(15.56) 


(w — wy) 
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The constants a, 8, y and 6 may be written immedately from (15.56): 
however, for purposes of simplification let 





1 
T= 
w, Co 
2 = 
Z=j 7 (o— wy) (15.57) 
21; . , ine 
M = R (w. — wi) (See equation 15.50) 
1 
then, 
sg = eae (15.58) 


J+ U4 


Now if Z, may be represented by W in (15.52) the remaining constants 
must be 


a=rTt+ 7M y¥=j) 
, (15.59 

fs 5=j) 

Now Z represents the frequency variable and graphically represents a line 

coincident with the J-axis in the Z-plane, and W’, the corresponding im- 

pedance variable, represents a circle in the W-plane. 

the center of the circle, Wo, in the W-plane is given by 


a 
W, == E (15.60)* 
iP y 
i 


when ¥ and 6 represent the conjugate functions ot y and 6 res,«ctively. 
The radius of the circle is given by 


The coordinates of 


ed | 


Finn 

ou (ee (15.61)* 
. v4 
. Ss 


*E. C. Titchmarch, “Theory of Functions,” Oxford 1932, pp. 191-192. Note: These 
equations are not derived in Titchmarch; however, by taking the limit as the diameter of 
the circle in the Z-plane approaches infinity, (15.60) and (15.61) result. 
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Now substituting the values given in (15.59) in (15.60) and (15.61) to deter- 
mine the coordinates of the center of the circle and the radius respectively, 
we find 
Wo = 
(15.62) 


As might be expected, the radius of the circle is equal to the real component 
of the coordinates of the center of the circle. This indicates that Fig. 15.10 
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RESIS TANCE 


— REACTANCE + 


ENLARGED VIEW 
OF AREA 
NEAREST ORIGIN 





Fig. 15.10—Circle diagram for crystal circuit combinations 


aa 


graphically represents the circle diagram where » = 7. The impedance, 


Z., is represented for a given frequency, as a vector from the origin to 


the corresponding point on the perimeter of the circle. From Fig. 15.10 
the following crystal properties may be deduced. 
1. The anti-resonant impedance designated as O — /P is simply ¢ + 
Vo? — y?, and when evaluated equals Ri M? — 1) = R,M’. 
2. The resonant resistance, O — P; given by ¢ — Vo? — 7? isequal toR;. 
3. The maximum positive reactance between w, and w: is represented by 
the distance from P; to Ps. It is given as o — » which equals 


2 
& M* _ i) | 
2 
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4, The condition which must exist when the crystal reactance is zero 
between w; and we occurs when ¢ = n or M = 2. 

It follows that the error of measuring, say the series resistance R, , by 
varying the frequency for maximum transmission and assuming true 
resonance when the crystal is between two low non-inductive resistors is 
associated with the difference between the length of the two vectors O — ?, 
and O — P,. The per cent error caused by the crystal capacitor, Co, by 
this method of measurement of R; is given as 


7 Pp 
Per cent error of R; = 100] 1 — —* | 
O — P; 


™ 100| 1 _ Ve ta =| 


o— Vo? — 9 


(15.63) 


If o and » are substituted in (15.63), we find 


ez 


Per cent error of R; = 100 (15.64) 
Py | 


This difference in amplitude was caused by the difference in frequency 
between resonance and minimum impedance. This frequency difference 
may be determined by transforming the points 7; and P; into the Z plane 
by (15.52) and subtracting them arithmetically. 

In order to express the coordinates of any point in the W-plane by its 
real and imaginary components let 


W = Ro t+ jXo (15.65 


Now the coordinates of ?; may be expressed as 
. ’ 


R= o[1- 4/1-(2 


The coordinates of P; are similarly given by 


T 1 


iii 1 /ie (2) 
_ 7 F al 


5 1 ia 


bo +) 


(15.67) 
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he coordinates of these two points represent values of W; now (15.52) may 
be solved for 2. 
‘ a— yW anti 
Z =~. (15.68) 
bw -— B 


Substituting in values for a, 3, y and 6 given by (15.59), we find 


° 21 , ‘ M(w — Wo) . a(n + Xo) 
Ww — ws) ] = —]2 


as * ae pee +X +R 


Z (15.69) 


The real component of Z must be zero since the function of Z is coincident 
with the Y-axis. By substituting values of Ry and No from (15.66) and 
(15.67), we find 


(15.70) 


M / 2 
[1 +(3) -1] 
o 
Subtracting (15.71) from (15.70) to get Aw we have 


M 


(Wo @)) 


Ae oe 
where 
V1 — (n/c)? — V1 + (n/a)? 
(n/a)? 


) Jo 
Now the lim A =1. When 1 = = then Aw = 2 
o M 


nio—0 
15.82 Circuit Analysis Involving Crystals 

The same procedure could be followed for the impedance, Z;, in Fig. 
15.3; however, the impedance expression conforming to (15.52) may be 
written directly if a more general expression is derived for impedances 
added in parallel or series. 

Paralleling the impedance, W, with an impedance, 7, (Fig. 15.11) modifies 
the constants in (15.52) but not its form providing 7 is essentially constant 
between w; and w:. For parallel impedances the impedance equation 
becomes, 

WT a+ BZ a’ + BZ 
W+T a B\, vy+8Z (15.74) 
+> T + r Z 
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The numerator remains the same; however, the denominator has an addi- 


tional term added to both y and 6. 

































































Fig. 15.11—Various crystal circuit combinations. 


In the same way, a series element, S, may be added as shown in Fig. 15.11. 
The impedance of this combination is similarly modified and the input im- 


pedance may be expressed as 


ven E +5 (, rn “)] +4 E +5 (s + )] Z 
si y 7 1 
W+T a B\, 
= 


_ a’ + BUS 
a >" + aA 
The addition of a series element leaves the denominator unchanged but adds 
a term to the numerator. The form of (15.75) is exactly the same as (15.52) 
except that the magnitudes of the constants, a, 8, y and 6 have been modified. 
As an example, apply (15.74) to obtain an expression for the impedance 
Zan (Fig. 15.2). This expression for impedance may be written directly 
if the following values are substituted in (15.74). 


uw 
~~! 
st 


a = 7+ jMr | 
B=r | 
y=b6=j (15.76) 
ee | 
fs 
- WT [7 + jMz] + [r]Z 
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From this expression we can again see the same form as (15.52) except that 
the coefficients of a, 8, y and 6 are modified by the capacitance, C,, in parallel 
with the equivalent crystal circuit. It has been previously explained that 
the frequency of oscillation in the generalized oscillator is determined by the 


anti-resonant frequency of the impedance, Z4,, and that this impedance 
represents the exact definition of P.I. P.I. therefore is represented in 
magnitude by O — P, in Fig. 15.10. 

Equation (15.75) enables us to write the expression for the input im- 
pedance, Z;, for the P.I. meter directly and furthermore, it enables us to 
compute the error produced by the adjustment of C, to minimum impedance 
rather than unity power factor as assumed in the derivation of (15.28). 
This error obviously will be a function of M, Cy, C, and R,. Writing the 
impedance expression for 4; directly, requires a little modification in that 
the crystal is shunted by two elements, C, (the crystal socket capacitance) 
and R, (the holder loss) as well as having the combination in series with 

CC) 
the capacitor, C3, where C; = ~,, ~.. Shunting the crystal with C, 
CC, + CG; 
modifies (15.77) only in that C, = C,; from this we can use (15.75) directly 
ina form readily adaptable to the determination of our original coefficients, 
a”, B”, y” and 6”. Adding a series capacitor, C3, and a shunt resistor, 
R,, modifies (15.77) in a manner specified by (15.75). Here, 


CN) 
t+ jMr s-i(1+¢) 


| 


gms T «R; : 
m4; 14 & hae it 

Y 4 ¢ 0 J ¢ ‘0 : jul 3 
Substitute these constants given by (15.78) in (15.75) then a”, 8”, 


6” become, 


ha re C:\ . MrX., 
wats(v+%)=[r+ x. (1+ G) 4 


MC.X., _ 
¢ 0 
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Substitute these values in (15.60) and (15.61) to obtain the coordinates of 
the center of the circle, Wo, and the radius, o. 


Co Co + Cz MrC, 
M 2 == 1 Phas, 
Wo — I( Co + Cs bin z)- _ 2i( Md C3 +5 2 + _R, ¢ 4] 


~~ MC, + 
ali+ (Zz) + Me ” a 


: 
j 
; 
; 
} 
(15.80) | 
i 
& 
f 
i 
: 





Co 
M 
ony ieee a f (15.81) 
1 : ey 
+(j i) "Rak 
From the above expressions, o is not exactly equal to the real component 
MC # 
of Wo. If however, 2 ss K - , then the above equations reduce to FF 
R, Co = c. +t 
Co ( Co + Cz 4, ; 
M ———~ — 2j{1+ —— + . : 
Wo = [Me na A, NC = Ri Cs /. (15.82) § 
ae 
P(r + ote) 
and : 
|, MC» | 
= — Gt Ge : (15.83) 
211 0 T 
[: ase eo | 


Now the radius, ¢, equals the real component of Wo. The per cent error 


of P.I. resulting from tuning to minimum impedance rather than unity 
power-factor in terms of Fig. 15.10 is given as, 


Fig. 
Me ee ‘. _O-P 
Per Cent Error of P.I. = 100 E Ox “J Re 
| E fe y/ -() | (15.84) 
= 100} 1 — e—__*____ wd 
L [¥s+(t)-1]] 
oO 
If 


Now since 7 represents the imaginary component of (15.82) then 


2o = M 


; bee ee, ‘) (15.85) 
(+ E\+§ pos ue 
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Fig. 15.12—-Inherent error of the P.I. meter due to tuning for an indication of minimum 
impedance rather than unity power factor. 


Rewriting (15.84) we have 


pie 
Per cent Error of P.I. = 100 — | 


[4/1+ 


M 


a 
. Co X, (15.87) 
1 = M : : 
(1+ O44 %) 


The error given by (15.86) is plotted in Fig. 15.12 as a function of P2. 


(15.86) 


If Cp > C, then 
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Fig. 15.13--Exterior view of crystal Performance Index meter. 
15.83 Frequency Errors 

As suggested in Section 15.81, conformal representation simplifies the 
mathematics required for the determination of frequency errors. In Section 
15.81, the difference between the resonant frequency and the minimum im- 
pedance frequency was computed for the equivalent crystal circuit. This 
same procedure could be used to compute the frequency difference between 
the antiresonant frequency of the generalized oscillator circuit (Fig. 15.2) 
and the minimum impedance frequency of the P.I. meter for the same value 
of C;. Comparison of the frequency of these two oscillators is plotted in 
Fig. 15.6 together with the measured values obtained from a ‘Pierce’ and 
a “Tuned Plate” oscillator. This frequency comparison involves setting 
for Aim (Fig. 15.2) and one for the impedance, 
Z, (Fig. 15.3) similar to Fig. 15.10. The impedance equations for both 74, 
and Z; would be arranged such that they have the same function of 7 in 
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order to have a common Z plane. In this way, transformation of operating 
points such as the “‘antitresonant frequency” operating point (/’s) for the 
Z,» impedance circle, could be subtracted from the ‘*minimum impedance 
frequency”’ operating point (/;) of the P.I. meter impedance circle in the 
common Z-plane. As in section 15.81, the frequency difference represents 
the arithmetic difference between P; and P; in the Z-plane in terms of 
2hi 
R; 


This curve was computed for the case when R, is negligible. The deriva- 


(w — w2). As an example, look at the calculated curve in Fig. 15.6. 


tion of (15.51) given in section 15.92 precisely follows the procedure just 
described. 

It is of interest to note that (15.27) may also be derived by Conformal 
means. It is more laborious than the usual circuit equations of section 
15.2; however, it does provide a check of the methods used. 


15.84 Errors of Other Approximations 


Further consideration of P.I. meter errors leads to the assumptions made 
in the derivation of (15.27). The derivation of (15.27) assumed that the 
resistor, Ry, was non-reactive. While actually it can be made essentially 
noninductive, we have neglected the effect of the input capacitance of the 
attenuator that is shunted across its terminals. The error from neglecting 
this capacitance in (15.28) is given by the following expression 

Per Cent Error = 100] 1 — V I i O' , (15.88) 
wR(CuQ + Ca) 

Where Q, equals the reactance of the shunt capacitance of the attenuator, 
C, , divided by the magnitude of the calibration resistor, Ry . 

It is interesting to note in the derivation of (15.28) that Ry was assumed 
to be very much less than Ve, which introduces an error of, 


R. - 
; a 


(15.89) 


/ 
Per Cent Error = 100] 1 — / soak, 


15.90 Derivation of Circuit Equations 
a 


15.91 Derivation of Equation (15.42) 

Other equations used in this paper may best be developed from Fig. 
15.3. By analysis of the input impedance, Z, , the basis for the development 
of (15.42) is as follows: 
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: 1 
ime ix ; [x (1 7 -c) | 
 jwlr * joo Tis 
[a+ jfon-*(h + YI] 
w \¢ 0 ( 1 


By substituting (15.48) and (15.49) in (15.90), we find 
171 , 1 oR + jLile” — oF : 

Z; aes = = 15.9 

jw E ™ C5 E + jLy(w? — wi) J] (15.91 


This may be expressed in the form, 


wR 4 jl, [= - Ws) re (w” es wi) | 
t 0 


From Fig. 15.2 


Ae 


(15.90 





Z; = H (15.92 
oad —wly(w? — ws) + jw? Ri 
where 
CoC; 
H = — ‘ 
Co + C; 


" . : 2 W2 —~ W) 
Now adding and subtracting w: to the -.——_ term we have 
0 


wk, : (w" — w) (w” — w>) (w> = a 
qt a ‘aac: nao 


—wly(w” — ws) + jw’ Ri 


Rationalizing (15.93) and equating it to R, and substituting in 1, = 


‘ 15.93) 
Z —_ ( : 


>, (obtained from (15.50)), we find 


Co(w? oP w}) 
oe =e w sf (15.94) 
at 
Ws — Wj w, M 


15.92 Derivation of Equations (15.51) and (15.44) 


Equation (15.51) makes possible the theoretical computation of the fre- 
quency difference between the generalized oscillator and the minimum im 
pedance frequency adjustment of the PJ meter. The derivation assumes 
that R, is negligible and that the total capacitance across the crystal ter- 
minals is lumped in series with the crystal. 

The impedance, Zs, in the generalized oscillator, Fig. 15.2, was given by 
(15.77). This equation may be expressed as follows: 


+jM)+Z 


MC, +i] +7 (15.95) 
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From this expression, as previously explained, (see Section 15.81), the 
following values for o and n may be determined. 
M Co } 
¢= ; - | 
204+C 4 (15.96) 
n =1 


The anti-resonant impedance of Z, is represented by O — P, in Fig. 15.10. 
The left hand term of (15.95) for the Ps operating point becomes 


Zapwl(Co+ Cy) =o + Ve _ (15.97) 


Substituting this value in (15.95) and solving for Z, we find 
Cot C 


f= 
: 


If Ks is expanded and all except the first two terms are neglected, 7 
may be expressed as 


Z= —j]M E Co _ | (15.98) 


where 


(1 + 4] 
—@, t+ i= (15.99) 
a 


The next step is to obtain a similar expression to (15.99) only for the 
minimum frequency impedance of Z; in Fig. 15.3 with p disconnected. For 


1 
this application S ae C and7 = x, Substituting these values, as well as 
js 
those in (15.59), in (15.75), we have 
; 1 | 1], 
Heist jer + — | +ir tte 
' w( ce t ‘ 


wl (15.100) 


J+ 94 
From this equation, values for o and 7 may be determined as described in 


Section 15.81. 


(15.101) 
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By the same procedure just described for evaluating Z from the jn- 
pedance expression Z,,, the value of Z corresponding to the minimum 
impedance operating point, P;, (Fig. 15.10) must be determined. For 
this operating condition 4; may be expressed by (15.65). The coefficients 
of this operating point are given by (15.67) with the above values of ¢ and 
9 (Equation 15.101). Utilizing (15.68) to solve for 4, we get 


se wee [ bi /' + | 
a: Co ) 
(: + 4 ? 


(minimum impedance) (15.102 


now 


: ot aan 
Z4-Z,;= (Aw) 
Ri 
where Aw = the difference in radian frequency between the frequency of 
oscillation in the generalized oscillator (anti-resonant frequency of the im- 


pedance, 4,4) and the frequency of oscillation in the PJ meter (minimum 
° ° ° . ry! 
impedance frequency of the impedance, Z;) 


a inf ‘0 2 mn af, oe 
oe TC ‘| Mw TG, 


+ y' + es 


(15.51) 


It is of interest to note that (15.102) becomes (15.44) when the value of 
Z (15.69) is introduced. 
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Lightning Protection of Buried Toll Cable 


By E. D. SUNDE 


A theoretical study of lightning voltages in buried telephone cable, of the 

liability of such cable to damage by lightning and of remedial measures, together 

with the results of simulative surge tests, oscillographic observations of light 

ning voltages and lightning trouble experience 

INTRODUCTION 

RACTICALLY all of the toll cable installed since 1939 has been of the 

carrier type and most of it has been buried in order to secure greater 
immunity from mechanical damage. It was realized, however, that bury- 
ing the cable would not prevent damage due to lightning and that, on ac- 
count of their smaller size, more damage was to be expected on the new car- 
rier cables than on the much larger voice-frequency underground cables 
then in use. Moreover, when damage by lightning does occur, such as 
fusing of cable pairs or holes in the sheath, it is not so easy to locate and 
repair as on aerial cables, since excavations may have to be made at a num- 
ber of points. Studies were therefore made of the factors affecting damage 
of buried cables by lightning and remedial measures were devised and put 
into effect in cases where a high rate of lightning failures was anticipated on 
new installations, or was experienced with cable already installed. Most 
of the cable installed was thus provided with extra core insulation, and shield 
wires were plowed in on many of the new routes. 

It was recognized early in these studies that more effective lightning pro- 
tection might be secured by providing the lead sheath with a thermoplastic 
coating of adequate dielectric strength and an outside copper shield, and 
that such cable might be required in territory where the earth resistivity is 
very high. This type of cable has recently been installed on a route in high- 
resistivity territory where experience has indicated that other types of con- 
struction would probably be inadequate and, since it has advantages also 
from the standpoint of corrosion and mechanical protection, it may be used 
also where lightning is not of such decisive importance. 

When lightning strikes, the current spreads in all directions from the 
point where it enters the ground. If there are cables in the vicinity they 
will provide low resistance paths, so that much of the current will flow tothe 
cables near the lightning stroke and in both directions along the sheath to 
remote points. The flow of current in the ground between the lightning 
channel and the cables may give rise to such a large voltage drop that the 
breakdown voltage of the soil is exceeded, particularly when the earth 
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resistivity is high. The lightning stroke will then are directly to the cables 
from the point where it enters the ground, often at the base of a tree. Fur- 
rows longer than 100 feet have been found in the ground along the path of 
such arcs. 

The current entering the sheath near the stroke point is attenuated as it 
flows towards remote points. Since a high earth resistivity is accompanied 
by a small leakage conductance between sheath and ground, the current 
will travel farther the larger the earth resistivity. The current along the 
sheath produces a voltage between the sheath and the core conductors, 
which is largest at the stroke point. This voltage is equal to the resistance 
drop in the sheath, between the stroke point and a point which is sufficiently 
remote so that the current in the sheath is negligible. Since the current 
travels farther along the sheath the higher the earth resistivity, this resist- 
ance drop will also increase with the earth resistivity. The maximum 
voltage between sheath and core is thus proportional to the sheath resistance 
and, as it turns out, to the square root of the earth resistivity. Carrier 
cables now being used are of smaller size and have a higher sheath resistance 
than full-size voice-frequency cables, and for this reason they are liable to 
have more lightning damage, particularly when the earth resistivity is high, 

To secure experimental verification of certain points of the theory pre- 
sented here, staged surge tests were made on the Stevens Point-Minneapolis 
cable, one of the first small-size buried toll cables to be installed. The 
results of these tests, which have already been published,’ are here compared 
with those obtained theoretically, on the basis of the earth resistivity meas- 
ured at the test location. Lightning voltages on this cable route were also 
recorded by automatic oscillographs and the results of these observations 
are also briefly discussed together with the rate of lightning failures experi- 
enced on this and other routes. 

The first part of the paper deals with voltages between the cable con- 
ductors and the sheath due to sinusoidal currents and surge currents. The 
second part deals with the liability of damage due to excessive lightning 
voltages and with certain characteristics of lightning discharges of impor- 
tance in connection with the present problem, such as the impedance 
encountered by the lightning channel in the ground, the rate of lightning 
strokes to ground and to buried structures and the crest current distribution 
for such strokes. In the third part remedial measures are discussed, to- 
gether with lightning-resistant cable. 


I. VoLTAGES BETWEEN CABLE CONDUCTORS AND SHEATH 


1.1 General 


Cable installed in the ground is designated “underground” when placed 
in duct, and “buried” when not in duct. Buried cable is sometimes pro- 
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vided with steel tape armor for protection against mechanical damage. 
While such armor may also reduce voltages due to low-frequency induction, 
mainly because of the high permeability of the steel, this is not true in the 
case of lightning voltages. The magnetic field in the armor due to lightning 
current in the cable is rather high, and the corresponding permeability fairly 
low. The armor resistance is, furthermore, quite high compared to that of 
the sheath, so that the effect of the armor may be neglected in considering 
lightning voltages. The tape or armor is usually separated from the sheath 
by paper and asphalt, but is bonded to the sheath at every splice point. 
Strokes to ground, or to the cable, may give rise to large currents in the 
armor and thus to excessive voltages between the armor and the sheath 
some distance from bonding points. The resulting arc may fuse a hole in 
the sheath or dent it, due to the explosive effect of the confined arc, and 
insulation failures may be experienced on this account. Such failures are 
not considered here since they are usually confined to a single point and are 
thus of less importance than insulation failures due to excessive voltages 
between the core conductors and the sheath, which may be spread for a 
considerable distance along the cable. 

For protection against corrosion, buried cables are usually jute-covered 
(asphalt, paper and jute) and in some cases have thermoplastic or rubber 
coating. The leakance of jute-covered sheaths is usually large enough so 
that the cable may be assumed to be in direct contact with the earth and 
the leakance is, furthermore, increased at the time of lightning strokes by 
numerous punctures due to excessive voltage between sheath and ground. 
This effect is large enough so that even rubber-covered cable may be re- 
garded as in direct contact with the soil in the case of direct strokes and 
sometimes also for strokes to ground in the vicinity of the cable, as discussed 
later. 

In order to calculate the voltage between the sheath and the core con- 
ductors of a buried cable, due to a surge current entering the sheath or the 
ground in the vicinity of the cable, it is convenient to consider at first a 
sinusoidal current. ‘The voltage due to a unit step current may then be 
obtained by operational solution and, in turn the voltage for a current J(¢) 
of arbitrary wave shape, by means of either one of the integrals: 


t 
VG) = | JI'(t — r)S(r) dr 
Jo 


() 
= [ J(t — r)S'(r) dr 


“0 
where S(¢) is the voltage due to unit step current and S’(¢) the time deriva- 
tive of this voltage. The second of the above integrals is more convenient 
in the present case. 
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Photographic observations indicate that a lightning discharge is usually 
initiated in the cloud by a so-called “stepped leader,” except in the case of 
ee : ‘ ~ : . 3 ‘ . a : 
discharges to sufficiently tall structures’ where this leader is initiated a: 


the ground end. After the leader reaches ‘the ground, or the cloud in the 
case of a tall structure, a heavy current “return stroke” proceeds from the 
ground toward the cloud at about ;' the velocity of light. The main surge 


of current in the return stroke, which usually lasts for less than 100 micro- 
seconds, may be followed by a low current lasting for '5 second or so. 
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Percentage of Lightning Strokes in Which Current Exceeds Ordinate 


Figure 1—Distribution of crest currents in lightning strokes. 
Curve 1: Currents in strokes to transmission line ground structures, based on 4410 
measurements, 2721 in U. S. and 1689 in Europe. 
Curve 2: Currents in strokes to buried structures, derived from curve 1. 


There may then be a second leader, which does not exhibit the stepped 
character of the first leader and always proceeds from the cloud, and a 
second return stroke. This may be followed by a third leader and so on, 
the average number of strokes in multiple discharges being about 4 and the 
average time interval between strokes about ;'5 second. Single-stroke dis- 
charges are, however, most common, discharges having more than 6 strokes 
being quite rare although discharges with as many as 40 strokes have been 
observed. 


The crest value of currents in lightning discharges varies over wide limits. 
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Measurements of current in the ground structure of transmission lines* 
indicate that a relationship as shown in Fig. 1 exists between the crest cur- 
rents and the percentage of discharges in which they occur. In the same 
figure is shown the crest current distribution for strokes to buried structures, 
which is derived in Part II from the curve for strokes to transmission line 
ground structures. Although measurements of wave shape are not exten- 
sive, they indicate that the current reaches its crest value in 5 to 10 micro- 
seconds and that it decays to half its maximum in 25 to 100 microseconds, 
the average being about 50 microseconds. An average wave shape is 
assumed in this investigation. 

In some 80 per cent of all lightning discharges the cloud is negative, so 
that the flow of current is from the earth toward the cloud. Further details 
about lightning discharges are summarized in recent surveys * which also 
contain an extensive list of references. 


1.2 Direct Strokes—Current Propagation Along Sheath 


As mentioned in the introduction and discussed further in Part II, a 
lightning stroke to ground may arc to a buried cable in the vicinity, in which 
case virtually all of the current will enter the sheath near the stroke point. 

When a sinusoidal current J enters the sheath at « = O and the sheath is 
assumed to extend indetinitely in opposite directions from this point, the 
sheath current at the distance x is given by the following approximate ex- 
pression 

J -r: 
I(x) => e¢ (2) 


where I’ is the propagation constant of the sheath-earth circuit and is given 
by the following expression, derived in a previous paper.” 


Wie. nas 1} 
r = s liw(iw + 1 px)|° (3) 


where: v = Velocity of propagation along sheath 
= (2/vx)’ meters per second 


vy = Inductivity of earth = 1.256-10 °° hy/meter 
« = Capacitivity of earth = €-8.858-10 ” fd/meter 
¢ = Dielectric constant of earth 


p = Earth resistivity, meter-ohms 
In deriving the above formula the resistance of the sheath is neglected in 
comparison with its external reactance, which is permissible for frequencies 
in the range of importance, and the sheath is assumed to be half buried, 
that is, with its axis in the plane of the earth’s surface. The latter assump- 
tion gives rise to a comparatively small error when the formula is applied 
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to cables buried at depths up to one meter or so, the propagation constant 
for a cable buried at infinite depth being larger than that given above by 
a factor of 1/2. 

It is assumed that the current is propagated from the cable up the light- 
ning channel with infinite velocity. The voltage between the cable con- 
ductors and the sheath obtained in this manner reaches a crest value after 
some 50 to 100 microseconds, or after the current in an actual lightning 
channel has traveled from the ground to the cloud. The error due to this 
assumption is thus probably quite small as regards the crest voltage, al- 
though the wave front will be somewhat slower when the actual velocity 
of propagation is considered. 


1.3 Direct Strokes—Voltage for Sinusoidal Current 


The current along the sheath gives rise to an electric force along the 
latter. The electric force along the inner surface of the sheath is given by: 


: Tee ae 
E(x) = I(x) = > se ™ (4) 
where s is the mutual-impedance of the sheath-earth and core-sheath cir- 
cuits. 

The latter mutual impedance is equal to the ratio of electric force along 
the inner surface of the sheath at any point, to the total current along the 
sheath at the same point, and for low frequencies equals the direct-current 
resistance of the sheath. It is given by the following slightly approximate 
formula :” 

. 1, . . \} » 
Z = R(iwy)’*/sinh (iwy)’ (5) 
where: w = 2nf and 
R = Unit length d-c resistance of sheath, ohms/meter 
y = vb/2raR 


6 = Thickness of sheath, meter 
a = Radius of sheath, meter 
v = Intrinsic Inductivity of sheath 


1.256 X 10° henrys/meter 

The current in the core-sheath circuit and the voltage between core and 
sheath due to an impressed field E(x) along the core (inner surface of 
sheath) are obtained from the following equations, which are the general 
solutions of the transmission line equation for the core-sheath circuit. 


J(x) = [A + P(wle** — [B + O(w)e™™ (6) 


U(x) 


ll 


KoA + P(x)le'* + KolB+ O(x)]Je™™ (7) 
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where yp and Ay are the propagation constant and the characteristic im- 
pedance of the core-sheath circuit and 
1 


P(x) IK E(x)e'" dx = 
_ “0 


z | 
2Ko 
z 1 


' a at 2 
Q(x) = xt E(x)e'" dx = — 


Since the current must be zero at x = 0, it is necessary that 1 = B. 
To make the curent vanish when x becomes infinity, B must equal 


° 
~ 


eo” “er +h 


between core and sheath becomes: 


With these boundary conditions the voltage 


2) a Sg " 
U(x) = > 


ror (Te = (10) 

The propagation constant Ty is much smaller than I’, and may be taken as: 
Py = [(Ro + iwls)ioC] 

(11) 


= : liw(iw + Ry Lo]! 


Ry = Unit length resistance of core-sheath circuit, ohms/meter 
Ly) = Unit length inductance of core-sheath circuit, hy/meter 
C = Unit length capacitance of core-sheath circuit, fd/meter 

to = (1/LoC)? 


1.4. Direct Strokes—Lighining Voltage at Stroke Point 


The largest voltage between sheath and core conductors is obtained for 
x = 0, and for this case (10) becomes: 


i J Z 
Qo) = — 
ow) 2T+To 


Ry'/sinh (iwy)! 
| aa \ Je 

(iw a 1 px)” a (tw fe Ry Li )" 
; a 


aT) 


For sufficiently high frequencies, so that (iwy)’ > 1,iw > 1 px, iw > Rolo, 


. . 4 . | ° 
and sinh (iwy)’ = 3 exp (iwy)’ expression (12) becomes 


U(O) = JRy° 
yee 
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For sufficiently low frequencies, so that (iwy): <1, iw < i/px, iw < Ry/L 
and sinh (twy)’ 2 (iwy)’'(1 + iwy/6), expression (12) becomes 
: FR 1 
U(0) = (14 


2a + 8 (iw)*(t + wy 6)’ 
where 
a= 2p, B - RoC 
For small values of time, corresponding to large values of iw, the functio; 
S’(t) defined before, as obtained by operational solution of (13) is “” 


1 VU% ) 4 
S’(t) = Ry’ . e (15 
( 4 v + % (4) 


For large values of time, corresponding to small values of iw, the functio 
is obtained by operational solution of (14) and equals: (Reference 10, pair 
542) 


S’ t) calf R 6 “hl /6t ) 1¢ 


where, with (6¢, ‘y)' = uU: 
h(u) = —ie” erf (iu) = —=e“ [ e dr (17 
V 0 


erf being the error function. 

Values of the function /(#) are given in Table I. 

In Fig. 2, curve 1 shows the function S’(t) calculated from (15), and 
curve 2 that calculated from (16), for a cable of 1.4” diameter, using con- 
stants as indicated in figure. The constants apply to a cable on which 
measurements have been made of the voltage between sheath and core 
conductors, at a location where the measured earth resistivity was 40) 
meter-ohms. The function S’(¢) corresponding to equation (12) is obtained 
with sufficient accuracy by drawing a transition curve, 3, between curves 
1 and 2. 

If the impedance ¢ is taken equal to the direct-current resistance R of the 
sheath and if the velocity of propagation along the sheath and along the 
core are assumed to be infinite, so that T = (iwa)’ and Ty) = (iw6)', the 
following expression is obtained 


nee 1\) (18 
SA 2(aé +B) Ant ¥ 


In the following it will be shown that (18) is accurate enough for practica! 
purposes. 

The wave shape of the current in lightning strokes may be approximated 
by an expression of the form: 


JI® = Ie" — &”. (19) 
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With a = 013-10", 6 = .5-10°, a current of the wave shape used in the 
measurements referred to above is obtained. This current reaches its 
crest in 10 microseconds and decays to its half-value in 65 microseconds, 
and is fairly representative of the average wave shape of lightning stroke 
currents. In the following, the voltages are for convenience referred to a 
crest value of 1000 amperes, which is obtained when J = 1150 amperes, 
the latter current being the initial value of each of the two exponential 
component currents included in (19). 


TABLE I 


, a , 
/ 2 2 
V hiu) = e™ [ e dr = “WV 
? j 


< “0 


= u when u < 11 


~ 


1 
an when u > 10 
= 28 


4283 
3014 
2232 
1782 
1496 
1293 
1141 
1021 
O845 
0630 
0503 


~ 


0 
0993 
.1948 
.2826 
.3599 
4244 
4748 
.5105 
5321 
5407 
.5381 


of 
2 
ia 
4 
a 
.6 
a 
8 
9 


—_ 
SKOU ee Wwiee 


~ 
> 


A more complete table for the range between « = 0 and u = 4 is published in Bericht- 
everhandlungen Akademie der Wissenschaften, Leipzig, Math-Phys. Klasse, Vol. 80, 1928, 
pages 217 to 223. 


In Fig. 3, the dashed curve shows the measured voltage and curves 1 and 
2 that calculated for the above surge current for two conditions. In cal- 
culating curve 1, S’(¢) was taken as curve 3 of Fig. 2, the voltage being 
obtained by numerical integration in accordance with (1); in calculating 
curve 2, z is taken as the direct-current resistance of the sheath and the 
velocities of propagation are assumed to be infinite, so that S’(f) is given 
by (18). In the latter case the following expression is obtained for the 
voltage by solution of (1): 


V(t) = la ‘h(v/at) — b *h(+/o0)| (20) 


2(a' + 8B’) 


where the function h(u), « = +/at or +/ dt, is detined as before. 
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Comparison of curves 1 and 2 of Fig. 3 shows that (20) is accurate enoug! 
for practical purposes, so that the voltage may be taken proportional to the 
direct-current resistance of the sheath. Since 8° is only 2.5 per cent oj 
a’, propagation in the core-sheath circuit may be neglected in comparison 
with propagation along the sheath-earth circuit, so that it is permissible 
to take the voltage proportional to the square root of the earth resistivity. 





3000 
~ 2000 
S 
ro) 
8 ’ 
gS s(t) 
i 
ro) 
1000 
0 
1 2 5 10 20 50 100 
t - Microseconds 
Figure 2—Approximate solution for S’(t). 
1: Calculated from formula for small times. 
) 


: Calculated from formula for large times. 

3: Transition curve giving approximate solution for S(t). 
Earth resistivity, p = 400 meter-ohms. 

Radius of cable, a = 1.75 cm. 

Sheath thickness, 6 = 2.4 mm. 

Sheath resistance, R = .92-10-3 ohms /meter. 

Core-sheath cap. Co = .96-10~° fd/ meter. 

Core-sheath resist. Ro = R  .92-10°3 ohm, meter. 
Velocity % = 2-108 meter ‘sec. 

Velocity 7 = 1-10* meter, sec. 


Furthermore, from (20) it is seen that when a and 6 are divided by the same 
factor k, so that the wave shape of the current remains the same but the 
duration of the current is increased k times, the voltage is increased +/ 
times. Thus, if the surge current had reached its crest value in 20 micro- 
seconds and its half-value in 130 microseconds, the voltage would be in- 


creased by 1/2, and the crest voltage would have been reached after 120 
rather than 60 microseconds. 
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If the breakdown voltage of the core insulation is assumed to be 2000 


volts, the above cable would be able to withstand a stroke current of about 
30,000 amperes before the insulation is punctured. From Fig. 1 it is seen 


80 


amperes 


oe 
peas 


120 


Microseconds 


Figure 3—Comparison of measured, shown by dashed curve, and calculated voltage 
between sheath and core conductors, shown by curves 1 and 2, for surge current as shown 
and cable constants as given in Fig. 2. 

1: Calculated from formula including skin-effect in sheath and finite velocity of propa 


gation. 
2: Calculated from formula based on d-c resistance of sheath and assuming infinite 


velocity of propagation. 


that in about 50 per cent of all strokes the crest current exceeds 30,000 


amperes. 
When there are two cables, each will provide shielding for the other, and 


> 


the shielding effect may be calculated as for shield wires (Sec. 3.3). Fre- 


quently the cables are of equal or nearly equal size and are close together. 

It is then accurate enough to use the parallel resistance of the two sheaths 
5 

in calculating the voltage, which will be practically the same in both cables. 








264 BELL SYSTEM TECHNICAL JOURNAL 


1.5 Direct Strokes—Lightning Voltages Along Cable 


It was shown above that with only a minor error the impedance s may be 
taken equal to the direct-current resistance of the sheath and that the 
propagation constants may be taken as 


I = (iwa)', To = (iwp)! 


With this modification expression (10) becomes: 


, J R j (iwa) dx } iwB) drys } 
B = - “ea _ 2 Ga) . 1) 
U(x) ‘8 Bre | (iw) (24 
The corresponding function S’ is:"” 
7 R 1 ; } —ar?4t } —Br2/4t 
S'(x,t) = —- ae —- pe ) (22) 
\ 2(a — B) (1) \ B 


The voltage due to a surge current J(¢), as obtained from (1), may be 
expressed as: 


, R ioe ee oe ere ; 
V(x, t) = eo la’ g(a’ x, 1) — 8’ g(B'x, 2)] (23) 


: i \ 
where, with o = a’y or B’x, 


t 1 } P 
g(o,t) = [ J(t - »( ) e° “dr (24) 
Jo 1T 


For a current as given by (19), the latter integral may be expressed in 
terms of error functions of complex arguments, for which, however, no tables 
are available at present. Curves for the function g, as obtained by numer- 
ical integration are shown in Fig. 4. 

When the core conductors are connected to the sheath at « = 0, the con- 
stants 4 and B of (6) and (7) are obtained from the following boundary 
conditions: At x = 0, V(0) = Osothat A = —B. As before, B = —Q(~). 
The voltage between core and sheath is then given by: 


, J RT y tad 
l (Xx) 2 > ° (e , — @ To ) 
2 I” ney r 
0 - 
(25) 
JR 7 (iwa) ir —(iwB) dr 
os - a , le 4 ee al (iwB | 
:- — I, 
In this case the derivative of the voltage due to unit step current is: 
f R ' 1 : —ar? r? 
Sox, ) = = meen (26) 
2(a — B) \at 
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Figure 4—Function g(o, t) = J(t — r){ — ] e-? 
0 aT 


J(t) = I(e*t — e*) 
a = 1.3-10',b = 5-108 J = 1150 amperes 


Ra 
4 . = 7 2» — g(Rix >7 
Vo(a, 2) a —p) Be 9 g(3° x, t)| (27) 
In Fig. 5 the crest values of V(x, ¢) and Vo(x, /), calculated for the cable 
considered before, are plotted against x, together with those observed in 
the tests. When the voltage (0, ¢) at the point where current enters the 


sheath is great enough to break down the insulation of a core conductor, the 
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latter will be in contact with the sheath by virtue of arcing. Under this 
condition, the voltage V(x, 4) between this conductor and the sheath wil! 
increase with distance along the cable as shown in Fig. 5. A maximum is 
reached a fairly short distance from the original fault, and beyond this 
point the voltage slowly decreases. After a puncture of the insulation where 
current enters the sheath, other failures may therefore occur, not neces- 
sarily at the point where V(x, ¢) is largest, but sometimes at points nearer 
or much farther away where the insulation may be weaker. A. single 
lightning stroke may thus cause insulation failures over a considerable dis- 
tance along the cable. 
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Figure 5—Comparison of measured variation of voltage between sheath and core, along 
cable, as shown by points and dashed curve, with calculated variation shown by curves 
1 and 2. 

1: Conductor not connected to sheath. 

2: Conductor connected to sheath at point where surge current enters sheath. 


1.6 Direct Strokes—Vollage Due to Long Duration Current 


As mentioned before, a current of low value and long duration may exist 
on the lightning channel after the main discharge. This current is usually 
of such long duration that the resistance of the sheath must be considered 
in calculating the current propagation along the cable. The propagation 
constant in that case becomes 


r = ((R+ iwc) = ? + r) (*) (28) 
, 2p 


R, L and G being the unit length resistance, inductance and leakance of the 
sheath-earth circuit. Neglecting propagation in the core-sheath circuit, 
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the voltage between core and sheath at the stroke point due to a sinusoidal 
current J; becomes, 
’ JiR 
U0) => =4,(2) Rp + RD 
2v 


2] 


The corresponding voltage for a unit step current J; is: 


V3(0; 2) = ne (2e) erf (Rt L)’ 
2v 


where erf is the error function. 
For large values of time, when Rt/L > 1, (30) becomes 


: L 
Vi(0, 2) =J7,R( 2 (31) 
= (25) 


The latter expression is valid when ¢ exceeds about 2 milliseconds and 
thus applies for the long duration current of a lightning stroke, since the 
latter usually lasts for about 100 milliseconds. For a current of 1000 
amperes, the core-sheath voltage for a cable of 1.4’ diameter is about 700 
volts. In many strokes the long duration current may be several hundred 
amperes, and a substantial voltage may then exist between core and sheath 
for .1 second or so. Thus, while this current component does not increase 
the crest voltage, it substantially increases the likelihood of permanent 
failure when the insulation is punctured by prolonging the current through 
the puncture. 


1.7. Strokes to Ground Not Arcing to Cable 


Let it be assumed that the current enters the ground at the distance y 
from a buried cable and that conditions are such that it does not arc to the 
latter. The flow of current in the ground gives rise to an electric force in 
the ground along the cable, and thus to currents in the sheath and to voltages 
between core and sheath. When the earth is assumed to have uniform con- 
ductivity, the earth potential at the distance r from the point where current 
enters the ground is given by: 


Ve = JQ¢(r) = Jp/2rr 
where 
Earth resistivity in meter-ohms 
r = (x + y°)’ = Distance in meters 
The sheath current and the voltage between sheath and core may in this 
e e . P ll 4 - 
case be obtained from published formulas,” provided propagation along 








268 BELL SYSTEM TECHNICAL JOURNAL 


the core-sheath circuit is neglected in comparison with propagation along 
the sheath-earth circuit, which is permissible. The voltage between sheath 
and core conductor differs by the factor s/Z from the voltage between 
sheath and ground as given in Table II, case 3 of the paper referred to, 
being defined as before and Z being the unit length self-impedance of the 
sheath-earth circuit. At a point opposite the lightning stroke, x = 0, 
; the voltage between core and sheath is in this case given by: 


Rr f° : R p 
UO, y) = J— (rie “dx = J - (Ty) (33) 
. Ld Ont Z2r , 
where Z = I°/G and G is the unit length leakance of the sheath-earth cir- 
cuit. The leakance is given by the approximate expression: 


1 
G =(2° log ) (34) 
T Ta 


a being the radius of the sheath and log = log,. 
The function ®(T'y) is given by the approximate formula: 


i+ Ty 


(Ty) = log Ty 


(35) 


Inserting (34) and (35) in (33), the latter expression may be written: 
UO, vy) = UO, aly) (36) 


where ('(0, a) = U(O) is the voltages when the current enters the sheath 
directly (vy = a) and: 


Ay) = (Ios me ”) / ve 1 Ta (37) 


where = (twv 2p)' = (iwer)’, 

The rigorous solution of the time function corresponding to (36) would 
be rather complicated. Since, however, A is the ratio of two functions, 
each of which varies logarithmically with I’, and thus varies only slightly 
with iw, an approximate solution is obtained by replacing iw with 1/¢ in 
(37). For instance, the solution of an operational expression p~" is ¢"/n! 
while the solution of p~" log p is [W(1 + m) + log 1/t\t"/n!, y being the 
logarithmic derivative of the gamma function. For representative values 
of nm and t(n < 1,¢ < 10“), p is less than 5% of log 1/t, so that a good 
approximation is obtained by replacing p by 1/¢ in log p, which in this 
illustration simulates the factor A(['y). With this approximation: 


V0, y, t) = VO, a, dAly(a/t)'] (38) 
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In Fig. 6 is shown the variation in the voltage calculated from (38), 
together with that observed in the tests referred to before. That the meas- 
ured decrease in the voltage is smaller than calculated is due to the fact 
that the earth resistivity at the test location increases with depth. Earth 
resistivity measurements made by the four-electrode method show that the 
resistivity is about 400 meter-ohms for electrode spacings up to about 20 
feet and then gradually increases, reaching about 700 meter-ohms at 300 
feet, 1200 meter-ohms at 1000 feet and approaching 1500 meter-ohms for 
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Figure 6—Reduction in voltage between sheath and core with increasing distance from 
cable to point where current enters the ground. 

1: Measured when remote ground representing cloud is at a distance of 1000 ft 

2: Calculated for uniformly conducting earth with remote ground at distance of 1000 ft. 

3: Calculated for uniformly conducting earth with remote ground at infinity. 


large electrode spacings. The measured variation in voltage with separa- 
tion is in substantial agreement with that calculated for an earth structure 
of this type in the manner outlined in Section 1.9. 
1.8 Discharges Between Clouds 

In considering voltages due to discharges between clouds, the lightning 
channel is assumed to parallel the cable. Due to magnetic induction, the 
lightning current will give rise to an impressed electric force along the 
cable sheath. Without much error it may be assumed that there is no 
impressed force outside the exposed section of the sheath and that the 


electric force in the exposed section due to a sinusoidal current J is E (x) = 
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JM, where M is the unit length mutual impedance of the lightning channel 
and the sheath. The resulting sheath current J(x) is obtained from (6 
and (7), when E is replaced by E’, ly by T and Ko by K, the characterist i 
impedance of the sheath-earth circuit. The constants .1 and B are found 
by observing the voltage between sheath and ground is zero at x = 5 2. 
The electric force along the core is given by E(x) = RI(x), and the voltage 
between the core conductors and the sheath is obtained by a second ap- 
plication of (6) and (7), the constants 4 and B being determined from the 
condition that the latter voltage must equal zero at x = s/2. The voltaye 
between the core conductors and the sheath at the distance x along the cable 
beyond one end or the other of the lightning channel projection on the 
cable is then: 


x2. 2 lor : Pez ; 
U(x) = fee s : . tk r ll el "| (30) 
2Z(T" — To)L Tos I's 


the sign of the voltage beyond one end of the channel being opposite to 
that beyond the other end. 

Since T > Ty, the last bracket term may be neglected. It was shown 
previously, that attenuation along the core-sheath circuit within a distance 
of one mile, which is representative of the length s, is quite small, so that 
1—e '" =Tys. With these modifications: 


JRMI°?s 
zr? — 7) ° 


rh. To: roe 
U(x) = si (40) 
The earth-return impedances M and Z are’ given by the following 
. . 12 
approximate expressions: 


Sh OE EE (41) 
2m (h” + y’)(iwa)’ 
P fi 

Ea log 2 a (42) 
2r altwa)* 


where a is defined as before, log = log, and: 


h = height of lightning channel above ground 
y = horizontal separation of lightning channel from cable 


II 


The expression for M holds when a(h? + y*)' > 5, a condition which is 
satistied in the important part of the frequency range. 
Inserting (41) and (42) in (40): 


ue TRe! |. ee 
x) = "Thee 3 
lade eae (t) ‘ ~ 
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where 


V 2hs 
” ” j ? 
(h + y) log \ = 


a(iwa)? 


Comparison with (21) and (23) shows that in this case: 


V(x, ) =~ Rat 


(8° x, t) (45) 
2a — p) 8" 


In the above solution, ~ was assumed constant. Actually it changes 
slightly with frequency and, for reasons mentioned before, it is accurate 
enough for practical purposes to replace tw with 1/¢ when calculating u. 

The maximum voltage is obtained at « = 0, i.e., at a point opposite one 
end or the other of the lightning channel, and comparison with (22) shows 
that this voltage differs from that obtained in the case of a direct stroke 
by the factor u, since B° <a’ so that the second term may be neglected in 
(23). The above factor has the following approximate value: 


sh 
4 
w+» 


p=.l (46) 

Since each cloud has an equal and oppositely charged image at the dis- 
tance / below the surface of the ground, the electric field between cloud and 
ground is substantially equal to that between the clouds when s = 2h. 
For a discharge to take place between clouds, rather than to the earth, the 
length s would, therefore, have to be less than 2h; so that with y = 0 the 
factor would not be expected to exceed u = .28. Thus, for the cable previ- 
ously considered, failures due to discharges between clouds would not be 
expected except for currents in excess of 100,000 amperes in g lightning 
channel approximately above and parallel to the cable. 

Maximum voltages of opposite signs are obtained at the two ends of the 
lightning channel, the voltage at the mid-point being zero. As the distance 
x from one end of the lightning channel increases, the voltage diminishes 
rather slowly in the same manner as shown in Fig. 5 for Vo(x, ¢). 


1.9 Stratified Earth Structures 


In the foregoing, the earth was assumed to have a uniform resistivity p. 
In many cases the average resistivity near the surface along a route may be 
substantially greater or smaller than the resistivity at greater depths, and 
this condition may affect the nature of lightning troubles, as will be shown 
below. 

The function Qo(r) appearing in equation (33) represents the earth 
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potential at a point due to unit current entering the earth at a distance ; 
from the point; i.e. Qo(7) is the mutual resistance between two points on thy 


earth’s surface separated by the distance r. 


the following expression : 


1 r . 
Qo(r) = 5 . lp2 + (p1 — pre?” | (47 
<7 
where 
pi = Resistivity of upper layer, meter-ohms 
p2 = Resistivity of lower layer, meter-ohms 
Yo = k/2d 
d = Depth of upper layer, meters 
k = Constant depending on the ratio p;/ps 


When r is sufficiently small compared to d, the above expression ap- 
proaches the limit Qo(7) = p:/2mr and when r is sufficiently large compared 
to d the expression becomes Qy(r) = p2/2mr. Expression (47) gives a fair 
approximation to the function Qo(r) as given by curves calculated from 
rather complicated integrals.” Earth resistivity measurements by the 
so-called “‘four-electrode measurements” are based on measurements oi 
Qo(r) and the results of such measurements may usually be approximated 
to the same degree of accuracy by (47) as by the curves applying accurately 
for two-layer earth, for the reason that the earth structure usually departs 
considerably from an ideal two-layer earth. 
constant & is about as follows: 


For various ratios p;/p. the 


10 1 a 
1.84 1.16 At 


pi/p2 = 100 
> k= 2 


Inserting (47) in (33) and proceeding as before, the voltage between core 
and sheath due to a stroke at the distance vy may be written: 


= V(0,a,1) ” A(v) + (pr — p2)u(y) 


V(O, y, ¢ | 
0, 9, p2 A(a) + (pi — po)u(a) 


(48) 


where I'(0, a, /) is the voltage for a direct stroke calculated for an equivalent 
earth-resistivity : 


pe = pda) + (p1 — podu(a) (49) 


and where 


Ay) = log [1 + Iy)/Ty] 


(50) 


In the case of a two-layer f 
horizontally stratified earth the function Qo(r) may be approximated by f 
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u(y) = | : eve dx 
» ¢ 


1 + y(yo + T) 
Mi y(yo + T) 


2 6 ™* X(y) when yoy>1 (53) 


=] 


when yoy <1 (52) 


In applying the above expressions, a rough value of p, is first assumed in 
calculating A(@) and u(a) and a more accurate value next obtained from 
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Figure 7—Reduction in voltage between sheath and core with increasing distance from 
cable to point where current enters ground. 
1: Upper layer of 400 meter-ohms and 30 ft. depth. Lower layer of 4000 meter-ohms 
and infinite depth. 
2: Uniformly conducting earth. 
3: Upper layer of 1500 meter-ohms and 30 ft. depth. 
Lower layer of 150 meter-ohms and infinite depth. 


(49). If the value of p, thus obtained differs materially from the assumed 
value, a second calculation may be required. In the expressions for \ 
and u the resistivity p, is to be used in calculating I’, the latter being taken 
as (v/2p.t)' where / is the time to crest value of the voltage as before. 

In Fig. 7 is shown the manner in which the voltage decreases with in- 
creasing separation for three assumed earth structures. The resistivities 
and the depth of the upper layer were selected such that the equivalent 
earth-resistivity, and thus the voltage in the case of a direct stroke, is the 
same in all cases and equal to 1000 meter-ohms. It will be noticed that 
in the case where the resistivity of the lower layer is high, the voltage due 
to a stroke at a distance of 200 ft. is 50°] and at a distance of 1000 ft., 25°; 
of the voltage due to a direct stroke. When the cable is small, insulation 
failures may thus be occasioned by strokes to ground at considerable 
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distances from the cable, although the resistivity near the surface up to 
depths of say 50 ft. is only moderately high. 

It is seen, however, that when the earth resistivity of the lower layer js 
low, failures due to strokes to ground not arcing to the cable are rather 
unlikely, even when the earth resistivity near the surface is rather high. 
On account of the higher surface resistivity, however, a greater number of 
strokes would be expected to arc to the cable for a given equivalent re- 
sistivity, than when the conductivity is uniformly distributed. On the 
other hand, many strokes which would arc to the cable if the earth were 
uniformly conducting may channel through the surface layer to the good 
conducting lower layer, so that the incidence of direct strokes is reduced 
on this account. Experience indicates that the latter factor tends to 
predominate, so that lightning damage is not ordinarily severe when the 
resistivity is low at depths beyond 20 ft. or so. 

In the case of discharges between clouds the coupling between the light- 
aing channel and the cable depends, in the frequency range of importance, 
to a great extent on the resistance of the lower layer. Thus, when the 
resistivity of the lower layer is very high the voltages may possibly give 
rise to insulation failures in the case of small cables, while this is not likely 
to occur when the resistivity of the lower layer is small or when the earth 
structure is uniform and of moderately high resistivity. 


1.10 Cables with Insulated Sheaths 


Assume that a short length Ax of insulated sheath is placed on the ground 
and that a voltage is applied between the sheath and a remote ground. 
When the applied voltage is greater than the breakdown voltage of the 
insulation, arcing to ground will take place at numerous equidistant points, 
provided the insulation and the earth are assumed to be uniform. The 
voltage between the sheath and adjacent ground increases from zero at a 
point where arcing takes place to a maximum value midway between two 
points at which arcing occurs, the maximum value being equal to the 
breakdown voltage of the insulation. Midway between two arcing points 
the potential in the earth (referred to infinity) may with negligible error 
be calculated as though the leakage current through the numerous arcs 
were uniformly distributed along the sheath. This potential in the ground 
would then be AJ/GAx, where AJ is the total leakage current and 1/GAX 
the resistance to ground of the sheath without insulation, G being the unit 
length leakage conductance. Midway between the arcing points the 
potential of the sheath to a remote ground is then: 


V = Vo + Al/AxG (54) 
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Let d/y/dx be the leakage current through the arcs and d/,/dx the leakage 
current due to capacity C between the sheath and the adjacent ground. 
For sinusoidal currents the following equations then hold, when Z is the 
unit length impedance and G the unit length leakance for a sheath in direct 
contact with the earth: 


dIy dl, 1 , , ~e 
a o = J 55) 
(2 r a) _ 
—(, + I,)Z = “ (56) 
dx 
es th = Vy (57) 
iw’ dx 


In the last equation it is assumed that the voltage between sheath and 
ground is equal to the breakdown voltage of the insulation, although this 
is not true in the immediate vicinity of the arcs. 

Eliminating V’ the following equation is obtained: 


d’Iy1 : dy 1 j “ 
i in 4) — b = 0 58 
é a" z) 7 (4 eo ) (99) 


— iwCG 
—_ . ; Y — 
ie Sal” ed G+ iat 


where: 


Equation (58) is satisfied when: 
Ip => Ape sd + Ba * 
ae " (59) 
I, = Aye “3 + Bye “i 


where I and I’; are the propagation constants for a sheath in direct contact 
with the ground and for an insulated sheath without breakdown, respec- 
tively. 


r = (GZ)', r, = (vz)! 

For a sheath of infinite length the By and B, terms vanish, so that: 
I(x) = Io + 1, = Age? + Ave *™ (60) 
The constants A) and A, are obtained from the foilowing boundary 

conditions: 
Atx = 0 I(x) = (0) = Ao+ Ai (61) 
Asx— « I(x) > Aye '* = . et (62) 
1 
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where A, = (Z/Y)! is the characteristic impedance of the insulated she.th 
without breakdown. 
From (61) and (62) 


V Vi 
A, = K, Ay = T(0) x x (63 
So that: 
I(x) = Toye Tt? + U8 (he — ety (64 
Ky 


For a rubber insulated cable the breakdown voltage Vo would be in the 
order of 30,000 volts and the characteristic impedance A, would be in the 
order of 100 ohms. The maximum current which could flow on the sheath 


sataaraptepcsaciees ot 


rue OAS IOS FS 


without breakdown, V9/A,, is then about 300 amperes, while in the case of F 
an average lightning stroke the current /(0) would be about 15,000 amperes F 


(i.e. 30,000 amperes total). The first term in (64) gives the attenuation 


along a sheath in direct contact with the ground. The current given by 
this term would diminish from 15,000 amperes to about 2,000 amperes F 
within a distance of 3 mile or so, for a typical lightning stroke wave shape 
and an earth resistivity of 1000 meter-ohms. At distances of several miles J 


from the stroke point the first term will vanish and the current will be de- 
termined by the second term, since exp (—I,x) will vanish much more 
slowly than exp (—TI'x). 

In the case of a stroke to ground the impressed electric force in the ground 
along the sheath is 


Ej(x) = —dV.(x)/dx (65) 
where J’, is the earth potential due to the lightning stroke current and is 
given by 
2n(x? + y?)! 


Instead of equation (58), the following equation is obtained for the cur- 
rents in the sheath 


d*Iy1 dy 1 5 ; 
(3 C In z) + (Z Y 12) = —£)(x) (67 


V(x) = (66 


Writing £o(x) = cE (x) + (1 — c) F(x), the solution of the latter equation 
may be written as the sum of two solutions of the form given by (6) and (7). 
After the constants Ao , By applying to the current J) and the constants -/, 
and B, applying to the current 7; have been determined from the boundary 
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conditions, in the same manner as before, the total sheath current may be 
written in the form: 


I(x) = cI,(x) + (A — o)I;(x) (68) 


where J, is the current for a grounded sheath, 7; the current entering a 
perfectly insulated sheath by virtue of its capacity to ground and c is given 
by: 


= V/V (0) (69) 


where V° is the potential difference between sheath and adjacent ground 
without breakdown at x = 0, which is substantially equal! to the earth po- 
tential. The above relationship for the constant ¢ is obtained by applying 
equation (57) at x = 0 to the general solution for J; , V" being given by: 


x.) 


G ./ ‘2 tm 
: | E,(x)e ' dx (70) 


V° 


G + 1w0e “0 


The voltage between core and sheath of an insulated cable may be written 
in a similar manner when /, is replaced by V,, the voltage for a cable in 
direct contact with the ground, and J; replaced by V; , the voltage for a cable 
insulated from ground. 

From (68) and (69) it will be seen that when V’.(Q) is much greater than 
the breakdown voltage of the insulation, the current entering the sheath is 
nearly the same as for a sheath in direct contact with the ground. Thus, 
when the earth resistivity is 1000 meter-ohms, and the stroke current 30,000 
amperes, the earth-potential at a distance of 30 meters (100 feet) is 160,000 
volts. The impulse breakdown of the insulation may be in the order of 
30,000 volts, so that the current entering the sheath will be substantially the 
same as for a cable in direct contact with the soil. When the earth-potential! 
at x = Ois only slightly larger than the breakdown voltage of the insulation, 
however, the current entering the sheath through punctures in the insula- 
tion is fairly small. 

In the above derivation the voltages and currents were assumed to vary 
sinusoidally which, of course, is a rather rough approximation in a phenome- 
non where breakdown occurs after the voltage reaches a certain instan- 
taneous value. While the derivation is not accurate, it does indicate under 
what conditions an insulated cable behaves like a cable in direct contact 
with the ground. 


1.11 Oscillographic Observations of Lightning Voltages 


To obtain data on the characteristics of lightning voltages in buried 
cable, five magnetic string oscillographs were installed for one lightning 
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season along a 50-mile section of the Stevens Point-Minneapolis route. 
The oscillographs, which were arranged to trip where the voltage exceeded 
100 volts, recorded lightning voltages due to some 600 strokes on 38 days, 
the Weather Bureau average being 33 thunderstorm days for this region in 
the same months. The character of the voltages varied widely from sharp 
transients of a few millisecond duration to slowly changing voltages lasting 
.2 seconds from one zero value to the next, voltages due to multiple dis- 
charges being quite common, the interval between voltage peaks in such 
cases being in the order of .1 second. Of the disturbances, 90°% lasted for 
more than .1 second, 50° for more than .4 and 10° for more than 1.25 
second, the maximum duration being 2.3 seconds. By way of comparison, 
the observed duration of discharges to a tall structure (3) were, in respec- 
tively 90, 50 and 10 of all cases, in excess of .08, .3 and .6 second, the maxi- 
mum duration being 1.5 second. The maximum voltage recorded was 940 
volts and was probably due to a stroke to ground near the cable. About 2°, 
of the voltages were in excess of 500 volts, most of these and the lower volt- 
ages being due to discharges between clouds, as indicated by the opposite 
polarity of the voltages at the two ends of the test section. The wave shape 
of the voltages at the ends of the section were much the same as at inter- 
mediate points, even for the sharpest surges recorded, the attenuation along 
the core-sheath circuit being quite small. It is possible that substantially 
higher voltages than observed may obtain in the case of severe discharges 
along a path parallel to and directly above the cable, and that such voltages 
may produce cable failure if the core insulation is below normal. 

While the oscillographs were arranged to trip on 100 volts, a smaller 
voltage was recorded in 40°% of all cases, as the peaks were too fast to be 
recorded by the type of oscillograph used. It is also possible that for this 
reason fast voltage peaks in excess of the maximum given above may have 
escaped measurement. 


II. LiGHtTNING TROUBLE EXPECTANCY 
2.1 General 


In estimating the liability of a cable to lightning damage, it is assumed 
below that once the core insulation is punctured, as it is likely to be at several 
points, at least one permanent failure will occur. The lightning trouble 
expectancy curves presented here thus give the number of times lightning 
damage is likely to occur, without consideration of the extent of the damage 
on each occasion. Each case of lightning damage usually involves several 
pairs and, based on experience, repair of each such case would require about 
four sheath openings. Damage due both to direct strokes and strokes to 
ground is included. Discharges between clouds have been neglected as a 
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source of lightning damage, however, as the voltages are likely to be in- 
sufficient unless the insulation is below normal. 

The curves of lightning trouble expectancy calculated here are significant 
only if troubles on a long cable route are considered over a period of several! 
years, So that the mile-years covered are in the order of 1000 or more. 


2.2 Incidence of Strokes to Ground 


To estimate the lightning trouble expectancy it is necessary to consider 
the incidence of strokes to ground in the vicinity of the cable, the number 














Figure 8—Map showing the average number of thunderstorm days per year. 


of such strokes that will arc to the cable and cause damage in this manner 
and the number that will give rise to failure without arcing to the cable. 
Magnetic link measurements (14) indicate that high tension transmission 
lines will be struck by lightning about 113 times per 100 miles per year, on 
the average, the minimum incidence in one year being about one half and 
the maximum about 1.6 times the average value. The above average inci- 
dence is based on observations covering about 1600 mile-years and applies 
for lines traversing areas where some 35 thunderstorm days are expected 
per year, as indicated by data issued by the U. S. Weather Bureau and 


collected over a period of 30 years,” and shown in Fig. 8. The above data 
on strokes to transmission lines may be used to estimate the rate of lightning 
strokes to ground, provided the width of the zone within which a transmis- 
sion line will attract lightning can be determined. 
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Based on laboratory observations on small-scale models’ "* *™) a | 
above ground will attract lightning strokes within an average distance on 
each side of the line which is about 3.5 times the height of the line when (he 
cloud is positive and about 5.5 times the height when the cloud is negative. 
When the average height of a transmission line ground structure above 
ground is taken as 70 feet, a 100-mile line will thus attract positive strokes 
(i.e., strokes originating from a positive cloud) within an area of about 93 
square miles and negative strokes within an area of approximately 14.5 
square miles. 

About 15° 7 of the strokes to transmission line ground structures have 
positive polarity,” ” so that the average rate of positive strokes to ground 
would be about 1.8 and that of negative strokes about 6.6 per square mile 
per year. The rate of positive and negative strokes to ground would thus 
be about 8.4 per square mile per year in areas where the yearly number of 
thunderstorm days is about 35, corresponding to about 2.4 strokes per square 
mile per 10 thunderstorm days. 

Based on the above date, the ratio of negative to positive strokes to 
ground in open country would be 3.6. The ratio of negative to positive 
strokes has been determined by various investigators in different ways.’ 
The ratio derived from measurement of field changes during thunderstorms 
varies between 2.1 and 6.5, that obtained from voltages observed in antennas 
is about 2.9, while point discharge recorder measurements indicate a ratio 
of 3.5 and the magnetization of basalt rocks struck by lightning indicates a 
ratio of 2.25. The above data were obtained in the temperate zone; in the 
tropics nearly all strokes have negative polarity. 


2.3 Arcing to Cable of Strokes to Ground 


As the stepped leader of a lightning discharge approaches the earth, 
charges accumulate in the ground under the leader and the resulting flow of 
current in the ground will give rise to a potential difference between points 
in the ground under the leader and remote points. Thus, when the tip of 
the leader has approached within 10 meters off the ground and the leader 
current is assumed to be as high as 500 amperes and the earth resistivity to 
be 1000 meter-ohms, a point directly under the leader will have a potential! 
of 8000 volts with respect to a remote ground. The potential gradient along 
the surface of the ground would, of course, be affected by the presence of a 
buried cable. The total potential involved is, however, so small that the 
effect of a buried cable on the path of the leader would be entirely negligible 
compared to the effect of irregularities in the surface of the earth. As the tip 
of the leader contacts the ground, the potential may be large enough so that 
the leader may arc to a cable located within 2 ft. or so of the leader. Only 
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aiter the return stroke is initiated, however, is the potential large enough so 
that arcing will occur for appreciable distances. 

When a stroke current J enters the ground, the electric force in the ground 
at the distance r from the stroke point is given by 


E(r) = Jp/2er (71) 


o being the earth resistivity. If the breakdown voltage gradient of the soil 
is ¢) , breakdown of the soil will take place until E(r) = eo, or for a distance: 


ry = (Jp/2ne)' (72) 


For a distance ro around the stroke point the soil may then be regarded as a 
conductor of negligible resistivity. 
The resistance encountered by the lightning channel in the ground is then 


p pe) \ - 
Ro = — = . 73) 
2rro QnJ . 


Measurements of the surge characteristics of grounds of fairly small 
dimensions” indicate that the breakdown voltage of the soil may vary 
from roughly 1000 volts per cm to some 5000 volts/cm. The data are, 
however, quite limited and there is no assurance that the above values 
represent the limits. In the first of the papers referred to, measurement 
was made of the resistance encountered when a current of 10,000 amperes 
crest value was discharged into the ground from an electrode suspended 
above the ground. The measured resistance is in satisfactory agreement 
with that obtained from (73) using the earth resistivity and breakdown 
voltage determined from other measurements in the paper referred to. In 
connection with the present study some small scale measurements were 
made of the breakdown voltage of sand between plane electrodes, and of the 
resistance encountered in discharges into the sand over point electrodes. 
These measurements indicated a breakdown voltage of 5000 volts per cm 
for dry sand having a resistivity of 3700 meter-ohms and 2400 volts/cm 
when sufficient water was added to reduce the resistivity to 100 meter- 
ohms. The measured resistance of point electrodes was a satisfactory 
agreement with that calculated from (73) on the basis of the measured re- 
sistivities and breakdown voltages. These experiments were made with 
currents having crest values from 1 to 50 amperes. 

Measurements of the breakdown voltage of various types of soil between 
completely buried spherical electrodes indicate substantially higher voltages 
than those given above, from about 11 to 23 kv/cm.** It is possible that 
for surface electrodes, breakdown at the lower voltage gradients occurs 
along the surface of the ground, so as to form a conducting plane of radius 
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ry. This circumstance would change the preceding formulas only to the 
extent that 24 would be replaced by 4. For a given eo, the resistance 
would then be about 25% higher. 

The radius ro is not necessarily the same as the distance across which 4 
lightning stroke would arc to a cable in the vicinity. Streamers will extend 
in various directions beyond 7) so that ionization of the soil increases the 
conductivity for a greater distance, ry being the effective radius of an equiva- 
lent hemisphere of infinite conductivity. 

The potential difference between the conducting sphere of radius rp and a 
point at the distance 7, is 


- Jp fi 1 Jpr, — % ri 
Vn = =2(- =) =o? (74 
2r \r ry 2r ryY% 


Let it be assumed that streamers extend beyond 7 until the average 
voltage gradient between ry and 7; equals e,. The potential difference is then: 


Von = ex(t1 — 10) (75) 
From (74) and (75): 
J e 5 
y= Pay, (76) 
2rro ey ey 


The latter expression applies when the field is assumed to have a radia! 
symmetry about the lightning channel. When a buried cable is present, 
however, this symmetry is disturbed. Calculations indicate that the po- 
tential of the cable at the point nearest to the lightning stroke will be less 
than 20% of the earth potential at the same point if the cable were absent. 
As a first approximation the cable may, therefore, be considered to have 
zero potential. The potential difference between the sphere considered 
above and the cable is then: 


_ Jel 


Voe = — 
Qr ro 


(77) 


If rz is the distance from the channel to the cable, the potential difference 
is also given by: 


Vow = e1(12 = ro) (78) 


ro = (i a “) (79) 
ey 


From (77) and (78): 


The arcing distance calculated in this manner will be the maximum, while 
that calculated from (76) will be the minimum distance. 

From measurements made of the effective corona radius of a conductor in 
air,” it is found that the latter may be determined on the assumption that 
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breakdown occurs until e¢) = 14000 volts/cem. On the other hand it is 
known that arcing between two conductors a considerable distance apart 
will occur when the average gradient is about 10,000 volts/cm. For air the 
ratio é/¢: is thus about 1.4. For breakdown in the air, the ionization need 
not be very dense in order that the corona envelope may be regarded as 
conducting as regards displacement currents. For breakdown in the earth, 
however, the ionization must be much more complete in order to provide 
substantially better conductivity than the soil. The ratio e9/e; is, therefore, 
likely to be greater, and has been taken as 2 in the following. 

If the latter ratio is assumed for soil, it is seen by comparison of (76) and 
(79) that the presence of a cable may increase the arcing distance as much 
as 50%. 

The arcing distance may be written in the form 


n= (Jp)'g:, m= (Jp)'gz (80) 


1 ; €o 1 ; ey 
os o = * t bay 
a (*) ej : q (1) (1 + *) ( ; 


With e/e,; = 2 the following values are obtained 
€9 = 250,000 500,000 volts/meter 
gm =1.6:10* 1.13-10°° 
q = 2.4:10° 1.7-10° 


Values of gq, and gq about 25% yreater than those givea above are ob- 
tained by assuming that breakdown does not take place in the soil, but that 
a conducting plane of radius 7 is formed by ionization of the air near the 
ground. Expression (71) is in this case replaced by E(r) = Jp/4r° and the 
expressions for g become 


= 1 ; €9 
ce iy dey ey 
1\) eo 
2 1 
q (i) ( = =) 


In the following g = 2.5-10 * has been taken as representative for low- 
resistivity soil and g = 1.5-10 ° for high-resistivity soil. When the current 
is expressed in kiloamperes, the corresponding values are .08 and .047. 

The distance to which a stroke may arc is accordingly taken as: 


where: 


(82) 


p < 100 meter-ohms p > 1000 meter-ohms 
r= 08 (Jp)' meter r= 047 (Jp)! meter 
= .26 (Jp)’ feet = 15 (Jp)' feet 


where J is in kiloamperes. 
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These values are, of course, of an approximate nature, and are only 
indicative of what may be expected under average conditions. In some 
cases the breakdown voltage of high-resistivity soil may be substantially 
lower than assumed, while that of low-resistivity soil may be noticea}ly 
higher. 


2.4 Crest Current Distribution for Strokes to Ground 


When the earth resistivity is taken as high as 5000 meter-ohms and the 
breakdown voltage of the soil is taken as high as 5000 volts/cm, the re- 
sistance encountered by the channel on the ground for a current of 25,000 
amperes is about 250 ohms. If the lightning channel were a long conductor 
already in existence at the initiation of the return stroke and capable of 
carrying the stroke current without being fused, the current would be 
propagated upward with the velocity of light and the surge impedance of 
the channel would be in the order of 500 ohms. Due to the resistance in 
the ground the current would then be some 30°% smaller than for a stroke 
to an object of zero resistance to ground. However, the lightning channel 
may not be regarded in the above manner, but as a conductor which is 
gradually prolonged at about 1/10 the velocity of light, and the impedance 
of the channel is then much larger, perhaps 5000 ohms. The surge im- 
pedance of a long insulated conductor having unit length capacitance 
is (1/Cv), » being the velocity of propagation. When energy is required 
to create the conductor, so that the velocity of propagation is reduced, the 
impedance is increased. Because of the high impedance of the channel, 
the resistance encountered in the ground may, therefore, be neglected as 
regards the effect on the crest current. The crest current distribution 
curve for strokes to transmission line ground structures may thus be used 
also in the case of strokes to ground, although a different distribution curve 
is obtained for those of the strokes to ground which arc to buried cable 
(Section 2.7). 


2.5 Failures Due to Direct Strokes and Strokes to Ground 


In calculating the number of failures due to direct strokes and strokes 
to ground, the earth is assumed to be a plane surface. A tree placed at 
random may attract a lightning stroke toward a cable or it may divert it 
from the cable and the net effect of a large number of trees along a route 
of substantial length is likely to be small. This is also true for variations 
in the terrain. 

When XN is the number of lightning strokes to ground per unit of area, 
and s the length of the cable, the number of lightning strokes on both sides 
of the cable within y and y + dy is: 


dN = 2Ns dy (83) 
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A lightning stroke at the distance y will cause cable failure when the crest 
current i exceeds a certain value which depends on the distance: 

i = f(y) (84) 


The fraction of all lightning strokes which has a crest current larger 
than i will be designated P)(7). The fraction of the lightning strokes dV 


which will cause cable failures is then: 


dn = dNP,(i) = 2NsPy(i) dy (85) 


© The number of cable failures along the length s due to all lightning strokes 
}to ground up to the maximum distance Y that need to be considered is: 


. 
n= ws [ P,(i) dy (86) 
0 


For the purpose of computation it is convenient to change the variable 
in the latter integral from y to i. With y = f ‘(i) = y(i), di = dy-f"(y), 
iy = f(0) and T = f(Y), the following integral is obtained: 


T 
n = 2Ns | PW — | y(i) Pi(i) ai| (87) 
19 

In (87), J is the maximum stroke current that needs to be considered and 
may actually be replaced by infinity, as will be evident later on. The 
current ip , which is the minimum current that will cause insulation puncture 
in the case of a direct stroke, may readily be determined from the breakdown 
voltage of the insulation and the calculated voltage between core and sheath 
per kiloampere, in the manner illustrated in Sectien 2.4. 

In order to evaluate the integral of (87), it is divided as follows: 


iy I 
“= 2vs| YPN _ i yi(i) Poi) di — | vo(i) Pi(i) ‘| (88) 


io a) 


When i < i;, failures of the cables will be due to arcing of the stroke 
to the cable and when 7 > i;, failures will occur before arcing takes place, 
due to the leakage current entering the sheath. Within each of the above 
two ranges the relationship of y to 7 is different and is designated (7) and 
yo(i), respectively. 

As already shown, failures due to arcing will take place when: 


yi) < q(pi)' (89) 


where q is defined as before and p is the earth resistivity in meter-ohms. 
In Section 1.7 it was shown that failures due to leakage current will 
occur when 
i > io/A(y) 
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where 


Ay) = log (' ee °) / 68 (1/Ta) 


I being the propagation constant of the sheath-earth circuit and a the radiy: 
of the sheath. The solution of the latter equation for y is: 


1 2ot\i 1 
) = Vo = , mig/t a S mio/t vin | 
y Vo\t re =| ee 1) = emote 4 (9 


where 


m = log (1/aT) 


I = (v/2pt)’ per meter 
vy = 1.256-10 ° henries per meter 
t 210 ‘sec. = time to crest of core-sheath voltage. 








When (89) and (90) are equated, the following expression is obtained: 


i ide ee (3) 1 
vq- ; 


The value of 7 which satisfies the latter equation is the value 7, detined 
above, and is shown in Fig. 9 as a function of mip for various values o/ 
(2t/vq")’. 

When (89) and (90) are inserted in (88), the latter integral may be ex- 
pressed as follows: 


2#\' 
n = 2Nsp’ (arin — H(i,)| + (*" G(i,, mi:)) (92 


where the current is in kiloamperes and: 
N = Number of strokes to ground per square meter 
s = Length of cable in meters 
q S08 when p < 100, g = .047 when p > 1000. 


I 
H(i) = -| i* Pi(1) di (93) 
Si Bos PT) 1 Pili) di . 
G(i, Miy) = emio ‘ae 1 ag i emo rr 1 (94) 


The first term in (94) equals VPo(J) = ve(Z)Po(T). Since Psi) is nega- 
tive, the above integrals will have positive values. 

The term g{H(io) — H(is)| of (92) gives the portion of failures due to 
direct strokes while the term involving the function G gives the portion of 
failures due to ground strokes not necessarily arcing to the cable although 


eet op eat eendnarrerce seme a A 


Ee cnntieaanress 





{ hey 
( able 


Kiloamperes 


; 


of 


to 


th 








LIGHTNING PROTECTION OF BURIED TOLL CABLE 287 


| they may do so (i.e. many of the currents in excess of i; may arc to the 
cable, although this is not essential in order to produce cable failure). 
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miofé 1/23 
ie olr—1)= (*) 
(93 iil 
If strokes to ground not arcing to the cable were*neglected as a source 
(o4) F of failures, the number of failures would equal: 
Na = 2Nsp'gH (io) (95) 
hega- If, on the other hand, the dielectric strength of the earth were assumed 
to be infinite, so that none of the strokes to ground would arc to the cable, 
soci the number of failures would equal 
on ot 
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Figure 10—Function H(i) when Po(z) is approximated by: 
Po(i) = exp(—ki), k = .038 per kiloampere. 
H(i) = kM (kiybe** + Bae erfc(Ri)4] 
From Fig. 1, it is seen that Po(7), as represented by curve 1, is nearly a 
straight line on semi-log paper and may, therefore, be approximated by: 


(i) Se" (97) 


With & = .038 per kiloampere, a straight line is obtained which coincides 
with curve 1 at i = Oand i = 100 kiloamperes, and this value of & has been 
used in the following. The functions H and G obtained with this approxi- 
mation are shown in Figs. 10 and 11. In obtaining these integrals, the up- 
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per limit 7 may be replaced by infinity, and it is also seen that the firs' 
term in (94) then vanishes. 

The lightning trouble expectancy as calculated from (92) is shown in 
Fig. 12 as a function of the earth resistivity for various sheath resistances. 
The curves are based on 2.4 strokes to ground per square mile, which is 
approximately the number of strokes per square mile during 10 thunder- 
storm days. The number of thunderstorm days per year along a given route 
is obtained from Fig. 8 and thus the number of times lightning failures would 
be expected during one year. 


2.6 Expectancy of Direct Strokes 


The incidence of direct strokes to the cable may be obtained from (96) 
with io = 0 kiloamperes. 
The number of strokes arcing to the cable is thus 


Na = 2Nsp'gH (0) (98) 
The cable will thus attract strokes within an effective distance. 
y = p'gH(0) (99) 
p < 100 meter-ohms p > 1000 meter-ohms 
y = .365 p meters y = .22 p meters 


a4 2 p feet 


7 p feet 


2.7 Crest Current Distribution for Direct Strokes 


The fraction of the strokes to the cable having crest values in excess 
of 7 is given by: 


P(i) = H(i)/H(O) 


ae (100) 
= 2 (*) e* + erfc (ik)! 
v 


and is shown by curve 2 in Fig. 1. It will be noticed that a buried cable 
attracts a greater proportion of heavy currents than a transmission line, 


because of the circumstance that heavy currents to ground arc for greater 
distances. 


2.8 Lightning Trouble Experience 


As mentioned before, lightning damage may be due to denting or to 
fusing of holes in the sheath, or to excessive voltages between the sheath 
and the cable conductors. Only the latter form of lightning failures have 
been considered here, since they predominate for cable of the size now being 
used, particularly in high-resistivity areas, and are likely to extend for a 
considerable distance to both sides of the point struck by lightning and are 
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thus more difficult to repair. For full-size cable in low-resistivity areas, 
however, insulation failures are more likely to occur as a result of sheath 
denting. For instance, along the 300-mile Kansas City-Dallas full-size 
cable route, where the earth resistivity is in the order of 100 meter-ohms, 
and where there are some 50 thunderstorm days per year, failures over a 


period of about 15 years have occurred about .5 times per 100 miles per 


year. Of these troubles 85°) were due to sheath denting as a result of 
arcing between the tape armor and the sheath. Based on (99), 100 miles 
of cable would attract lightning strokes within an area of .5 square mile, 
so that the cable would be struck about six times per 100 miles per year, 
when the number of strokes per square mile per year is 2.4 per 10 thunder- 
storm days (Section 2.2). The rate of lightning failures experienced on this 
route may thus be accounted for by assuming that about 7% of the strokes, 
i.e. currents in excess of 90 kiloamperes as obtained from curve 2, Fig. 1, 
will produce sheath denting severe enough to cause insulation failure, while 
about 1%, ie. currents in excess of 140 kiloamperes, will cause insulation 
failure due to excessive voltage. The latter value is in substantial agree- 
ment with that calculated for a full-size cable when the earth resistivity is 
assumed to be 100 meter-ohms. It is evident from the above examples 
that in low-resistivity areas lightning troubles will not be a problem, and 
this is also borne out by experience on other routes installed in such terri- 
tory during the last few years. 

All cable installed in high-resistivity territory since 1942 has been provided 
with doubled core insulation and with shield wires, in spite of which con- 
siderable damage has been experienced on some routes, as between Atlanta 
and Macon. This appears to have been due partly to the circumstance 
that in many cases the insulation in splices and accessories has not been 
equal to that obtained in the cable through the use of extra core wrap, and 
that in some cases damage has been due to holes fused in the sheath due to 
arcing between the sheath and the shield wires. As an example, along the 
Atlanta-Macon route there are some 70 thunderstorm days per year and 
the average effective earth resistivity is about 1300 meter-ohms. The 
corresponding estimated rate of direct strokes to the cable is about 20 per 
100 miles per year. It is estimated, in the manner outlined in section 3.0, 
that only stroke currents in excess of 80 kiloamperes are likely to damage the 
cable, so that on the basis of curve 2, Fig. 1, cable failures would be ex- 
pected to occur about 2 times per 100 miles per year. The actual rate of 
trouble experienced on this route over two years has been about five times 
higher, so that some 50°; of the strokes to the cable; i.e., currents in excess 
of 30 kiloamperes or so, appear to have caused cable failures, most of which 
occurred in splices and accessories. 
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It is evident from the above examples that careful examinations o/ 
trouble records are required before the observed rate of lightning failures 
can be adequately compared with that obtained from theoretical expectancy 
curves. If the cable as well as splices and accessories actually have a di- 
electric strength as assumed in the calculations, it is likely that the average 
rate of failures due to excessive voltages experienced over a long period wil! 
not be any greater than estimated from these curves. 

Based on experience, an average of 4 sheath openings is required to repair 
damage caused by excessive voltage between the sheath and the cable con- 
ductors, as compared to about 2 sheath openings when the damage is due 
mainly to denting and fusing of the sheath, as in the case of full-size, tape- 
armored cable in low-resistivity territory. Although the damage may be 
confined to one point, it cannot usually be located by a single sheath opening. 


III. REMEDIAL MEASURES 
3.1. General 


From Fig. 12 it is evident that the rate of cable failures to be expected, 
and hence the need for remedial measures, depend greatly on the earth re- 
sistivity. Experience has indicated that lightning damage is likely to be 
encountered even when the surface resistivity is fairly low, provided the 
resistivity beyond depths of 10 or 20 ft. or so is very high. Considerably 
less trouble has been experienced where the resistivity below this depth 
is low, even where the surface resistivity has been high. The lightning 
stroke may then channel through the surface layer to the good conducting 
lower layer, so that direct strokes are not experienced as frequently in 
spite of the high surface resistivity. As a guide in applying protective 
measures, earth resistivity measurements are usually made along new 
cable routes.” 

The curves given in Fig. 12 may also be used to find the lightning trouble 
expectancy when extra core insulation, shield wires or both are used. Thus 
when the insulation strength is doubled the effect is the same as if the 
sheath resistance is halved. If the shield wires reduce the voltage by a 
shield factor n, the effect is the same as if the sheath resistance is multiplied 
by n. Considering direct strokes only, curve 2 in Fig. 1 may be used to 
find the percentage reduction in lightning strokes that will damage the 
cable, when the stroke current which the cable is able to withstand is in- 
creased by extra insulation or shield wires. 


3.2. Extra Core Insulation 


One method of reducing failures caused by lightning strokes to buried 
cables is to increase the insulation between the cable conductors and the 
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sheath, no extra insulation being required between individual cable con- 
ductors. This has already been done for most new installations. The 
cable itself, cable stubs, loading cases, and gas alarm contactor terminals 
are all provided with sufficient extra insulation to double the dielectric 
strength between cable conductors and sheath. For a cable like that on 
which the measurements referred to before were made, such a measure 
would increase the stroke current which would damage the cable from 
30,000 to 60,0000 amperes and would reduce the number of direct lightning 
strokes that could cause failure by direct arcing to the sheath to about 20 
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Figure 12—Theoretical lightning trouble expectancy curves showing number of times 
insulation failures due to excessive voltages would be expected per 100 miles for 10 thunder- 
storm days, for cables having sheath resistances as indicated on curves. Dashed line 
represents full-size cable. 
per cent of the total instead of 50 per cent (see Fig. 1, curve 2). The num- 
ber of failures due to direct strokes would thus have been reduced 2.5 times. 


3.3. Shield Wires 


Another method, employed in addition to the extra insulation where 
excessive lightning damage would otherwise be expected, is to bury shield 
wires over the cable. These conduct away part of the lightning current 
and thus reduce the amount that flows along the sheath. These wires 
may be plowed in with the cable, as has been done on several new routes, 
or may be installed afterward. When the wires and cables are plowed in 
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together, the arrangement shown in Fig. 13 has been used. The percentage 
of the current carried by the wires depends to a greater extent on their 
inductance relative to that of the sheath, than on their resistance. Two 
wires are employed, rather than a single wire of smaller resistance, in 
order to obtain a lower inductance than with a single wire. 

On the route between Stevens Point and Minneapolis, where the shield 
wires were installed after the cable was in place, two 165-mil wires about 
twelve inches apart were plowed in some ten inches above the cable for a 
distance of eighty miles. Surge measurements made after these wires were 
installed indicated that the wires reduced the voltage between sheath and 
core conductors about 60 per cent, in substantial agreement with theoreti- 
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Figure 13—Position of shield wires, S, when they are plowed in with cable, C. 
cal expectations. The cable would then withstand 75,000 amperes rather 
than 30,000 without shield wires. If a cable of normal construction were 
used on the above route, the shield wires would thus reduce the number 
of direct lightning strokes that would be expected to cause insulation failure 
about 3.3 times, from 50 per cent of the total without shield wires to about 
15 per cent with shield wires (see Fig. 1, curve 2). If the insulation 
strength is only 1000 volts, however, about 80 per cent of all strokes would 
be expected to cause failure without shield wires and 40 per cent with 
shield wires, so that the reduction would be substantially smaller, as ac- 
tually appears to be the case on the above route. 

Aside from reducing core insulation failures, shield wires also minimize 
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damage to the sheath covering in the case of strokes to ground near the 
cable, particularly in the case of thermoplastic or rubber-covered cables. 
As mentioned in section 1.10, considerable current may enter the sheath of 
an insulated cable in the case of strokes to ground near the cable, due to 
numerous punctures in the insulation of the sheath. If conditions along 
the cable were uniform, the current through each puncture would be so 
small that the insulation would not be damaged. Due to variations in the 
resistivity of the soil and the presence of buried metallic structures, how- 
ever, concentrated arcing may occur, and the insulation may then be 
damaged in spots, so that corrosion may be initiated. Shield wires reduce 
the voltage between the sheath and ground and thus the likelihood of 
damage to the sheath insulation. 

It can be demonstrated theoretically, and has been proved by measure- 
ments, that when current enters shield wires next to a buried cable in good 
contact with the earth, the current in the sheath and the voltage between 
sheath and cable conductors is negligible. The reason for this is that 
current induced in the sheath by the shield wire current is equal and oppo- 
site to the current entering the sheath by virtue of leakage through the 
ground. Negligible voltages between sheath and core conductors would 
thus be obtained if shield wires were installed at such a distance that light- 
ning strokes would be intercepted and direct strokes to the cable prevented. 
To prevent arcing to the cable of strokes to the shield wires, the separation 
between cable and shield wires would have to be at least 6 feet when the 
earth resistivity is 1000 meter-ohms, and greater for higher resistivities. 
Such wires cannot, therefore, be as easily installed in one plowing operation 
with the cable as shield wires at a smaller spacing. They are, therefore, 
not considered here, although they have been installed in one instance in a 
fairly short section where repeated lightning damage had been experienced. 

When the sheath, as well as the shield wires, is in intimate contact with the 
earth, the propagation constant for the shield wires is the same as that for 
the sheath. In the case of a direct stroke to the cable or the shield wires, 
the voltage between sheath and shield wires will be so large that they will 
be in contact with each other at the stroke point by virtue of arcing. The 
current in the sheath is then: 


I(x) = (101) 


J Zn — Zy et 
2 Zn + Zen — 222 


where J is the total current at « = 0, and 
Zu = Ri + iwl, = Unit length impedance of sheath 
Ze = R,+ iwlx» = Unit length impedance of shield wires 
Zp = twlipo = Unit length mutual impedance of sheath and shield 
wires 
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\ 
The voltage between sheath and core conductors at x = 0 then becomes: aa 
—— 4 [ an 
V(0) = J =. ae a 2a = i) : : 
2 at + BZ + Ze — 2Z2 \iw (102 Ban 
+) #e ist 
2a Bt \iw) p + Bo ‘ “ 
where mob 
n = (Le — Ly)/(Ln + Le — 212) (103) t a 
fm ab 
a = R2/(L2 — Li) (104) FE th 
Bo = (Ri + Ro)/(Ln + Le — 2112) (105) 


ret 
The function S’ is in this case 


, = fe iY . = =-foas ] 
S’(0, t) = 2 at + Bl n |(4) + ai h(B) t’) (106) 


where the function / is defined as before. or 
When the shield wires are at a sufficient distance from the sheath, so that s 6th 

proximity effects may be neglected, the self and mutual inductances are 

as follows: 


vO 
wi 


dt 
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In — Ly vs log a (107) ‘ 
la ~ Ln = log (108) . 
Bs eet "a. reo Pp 
where . ze. 
log = log, and: , titi 
v = 1.256-10°° henries per meter ee. 

ru = Radius of sheath 

tx = Radius of shield wire 

ry. = Distance between sheath and shield wire. 

With more than one shield wire, rx. is the geometric mean radius and ry» 
their geometric mean separation from the sheath. When there are two M 
cables, as is frequently the case, ry, is the geometric mean radius of the cables t 
and R is their combined sheath resistance. c 

The surge voltage obtained by solution of (1), for a current as given by 
(19), is in this case: 

G 

IRn a — © 477448 
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When ap = (§o, the voltage with shield wires differs from that obtained 
without shield wires by the factor 7. With two 165-mil wires 12 inches apart 
and 10 inches above the cable considered before, a9 = 1.3-10°, 8) = 1.5-10° 
and» = .47. In this case ap differs only slightly from 8» so that the voltage 
is reduced by the factor n. Measurements made before and after the shield 
wires were installed indicated a reduction factor of .40 (i.e. the voltage with 
was .4 times the voltage without shield wires). The reasons for the smaller 
observed factor is partly that the shield wires are in more intimate contact 
with the earth than the sheath, and partly that the resistivity of the soil 
above the cable, where the shield wires are located, is somewhat smaller 
than the resistivity at the depth of the cable. 

With 104-mil wires, a9 = 3.2-10°, 8) = 2.3-10° and » = .49. When the 
reduction factor is determined more accurately by calculating the crest 
voltage with shield wires from (109) and comparing it with the crest voltage 
without shield wires, a value of .52 is obtained as compared with .47 for 
165-mil wires. It is thus seen that, within certain limits, the voltage re- 
duction provided by shield wires depends to a comparatively small extent 
on the size of the wires, the resistance of 104-mil wires being about 2.5 times 
that of 165-mil wires. 


3.4 Lightning-Resistant Cable 


As mentioned before, buried cable may be covered by jute, thermoplastic, 
or rubber for protection against corrosion. The coating may be damaged 
by gophers or by lightning, and severe corrosion may be experienced at 
points where the coating is ruptured, particularly when thermoplastic, or 
rubber coating is used. Even when the earth resistivity is low and protec- 
tion against core insulation failures due to excessive voltage would not be 
required, the sheath coating may be damaged rather frequently. 

Shield wires may effect a substantial reduction in core insulation failures 
and may also prevent damage to the coating in the case of strokes to ground 
at some distance from the cable. In the case of direct strokes, however, 
arcing between the shield wires and the sheath will damage the sheath coat- 
ing and may also fuse a hole in the sheath, although there may be no insula- 
tion failures due to excessive voltage between the sheath and the cable 
conductors. 

Reduction of damage to the coating and to the sheath occasioned by light- 
ning, rodents, or corrosion and protection against core insulation failures 
occasioned by excessive voltage or crushing of the sheath may be secured 
by providing the sheath with a thermoplastic or rubber coating and an out- 
side copper shield. If various auxiliary equipment connected to the sheath, 
such as load coils, gas pressure contactors and terminals are also properly 











298 BELL SYSTEM TECHNICAL JOURNAL 


insulated from ground, currents will not then enter the sheath, except 
through the capacitance to the outside shield. ‘The voltage across the core 
insulation will then be so small that core insulation failures will not occur, 
unless the voltage between the outside shield and the sheath is large enough 
to puncture the coating. Thus, if the resistance of the outside shield were 


1.1 ohms per mile (10-mil copper shield around 1.2” cable) and if the impulse 


breakdown voltage of the coating were 20,000 volts, breakdown of the in- 
sulation would not be expected except for currents in excess of 130 kilo- 
amperes when the earth resistivity is as high as 4000 meter-ohms. For a 
cable of 2” diameter with a 10-mil copper shield, breakdown of the thermo- 
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Figure 14—Thermoplastic covered, copper jacketed cable. 


plastic insulation would not be expected except for currents in excess of 200 
kiloamperes when the breakdown voltage of the coating is 20 kv and the 
earth resistivity is 4000 meter-ohms, or when the breakdown voltage is 
10 kv and the earth resistivity 1000 meter-ohms. The above type of cable 
is also advantageous in that low-frequency induced voltages and noise due 
to static are reduced. 

A cable of the above type is now being installed for a distance of about 
180 miles along the Atlanta-Meridian route, where the effective earth re- 
sistivity varies between 1000 and 4000 meter-ohms and lightning storms 
occur frequently. The diameter of the lead sheath of this cable is about 
2”, and the sheath is covered with two layers of thermoplastic each 45 mils 
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thick, with diametrically opposite seams having an overlap of about }” 
The outside 10-mil copper shield is corrugated to facilitate bending and has 
a }” overlap at the seam. The thermoplastic coating is flooded with 
thermoplastic cement to reduce moisture absorption. A photograph of this 
cable is shown in Fig. 14. 


SUMMARY 


Current in the sheath of buried cables, due to direct lightning strokes, 
strokes to ground near the cables or discharges between clouds, gives rise to 
voltages between the cable conductors and the sheath. The voltages are 
practically proportional to the square root of the earth resistivity and to 
the direct-current sheath resistance. For the latter reason they are sub- 
stantially larger for carrier cables of the size now used than for the much 
larger voice-frequency cables. While direct strokes are usually most im- 
portant, strokes to ground must be considered when the cables are of small 
size, even when the surface resistivity is low, provided the resistivity at 
greater depths is high. Under the latter conditions it is possible that for 
small cables, discharges bet ween clouds over the cable may also cause failure. 

For cables with thermoplastic or rubber coating the voltages between 
the sheath and the core conductors are much the same as for jute-covered 
cables. The coating of such cable is likely to be damaged by direct strokes 
and strokes to ground near the cable, in which case corrosion of the sheath 
may occur at such points. 

Based on theoretical lightning expectancy curves, the incidence of light- 
ning troubles increases faster than the sheath resistance or the earth re- 
sistivity. When the breakdown voltage of the core insulation is doubled by 
use of extra core wrap, or when shield wires are installed in situations where 
lightning damage is anticipated or has been experienced, a substantial re- 
duction in lightning failures is to be expected. Shield wires will, however, 
not prevent damage to the sheath and the sheath coating. 

Where the earth resistivity is very high and lightning storms occur fre- 
quently, doubled core insulation together with shield wires may not provide 
sufficient protection, even for cable of substantial size. Protection against 


various forms of lightning damage may then be secured by use of thermo- 
plastic sheath coating of adequate dielectric strength together with an out- 
side concentric copper shield. 
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Abstracts of Technical Articles by Bell System Authors 


Ultra-Short-Wave Receiver for the Cape Charles-Norfolk Multiplex Radio- 
telephone Circuit... D. M. Biack, G. Ropwin and W. T. WIntTRINGHAM. 
The requirements for an ultra-short-wave receiver for use in a multiplex 
radiotelephone link circuit are outlined. The technical details of a receiver 
designed to meet such requirements in the circuit between Cape Charles 
and Norfolk, Virginia, are described. 

Ultra-Short-Wave Multiplex? Cuartes R. Burrows and ALFRED 
Dectno. The technical requirements of a twelve-channel ultra-short-wave 
multiplex system are discussed and the means of meeting them are described. 
The intermodulation between channels in equipment based on this design 
has been reduced to the point where it is possible to use twelve-channel 
radio systems in the toll plant. By employing a sufficient amount of en- 
velope feedback, the transmitter can be operated with a high modulation 
factor without the use of spread sidebands. 

Airplane Vibration Reproducer’ G. R. CRANE. This paper describes 
a reproducer set designed for use in the reproduction for analysis of multiple 
track film recordings. It is capable of reproducing simultaneously 13 
variable-area tracks recorded side by side on standard 35-mm. film. Re- 
corded signals between 5 and 3000 cps are accurately reproduced and may 
be analyzed for frequency components, amplitude, and phase relation. 

Airplane Vibration Recorder J. C. Davipson and G. R. Crane. This 
paper describes a portable film recorder capable of simultaneously recording 
13 variable-area tracks on 35-mm. film. It is intended for use in the analy- 
sis of airplane vibration or similar studies in which it is desirable to record 
disturbances (mechanical, acoustical, or electrical) from a number of sources 
in such a manner that the resultant record can be analyzed for frequency, 
amplitude, and phase relation. Film speeds of 12, 6, or 3 in. per sec. are 
available. 

Application of Sound Recording Techniques to Airplane Vibration Anal- 
ysis® J. G. Frayne and J. C. Davipson. This paper describes methods 
which have been developed for analysis of the various vibration components 
present in airplane structures. The complex wave forms are recorded on 
standard motion picture sound negatives during flight. These films later, 
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after proper development, are analyzed electrically, making possible a 
complete analysis on the ground and thereby reducing materially the time 
devoted to flight test, and also simplifying the process of analysis of complex 
wave forms. 

Ultra-Short-Wave Transmitter for the Cape Charles-Norfolk Multiplex 
System.® R. J. Krrcuer and R. W. Frus. Design features of an unat- 
tended ultra-short-wave double-sideband multiplex transmitter are de- 
scribed. Forty decibels of envelope feedback is utilized over the 12- to 
60-kilocycle band of the twelve type-K carrier-signal channels which 
modulate the last stage of the transmitter. Accessibility of apparatus 
and ease in maintenance contribute toward obtaining maximum reliability 
of the equipment in commercial service. 

Paper Capacitors Containing Chlorinated Impregnants. Stabilization 
by Anthraquinone.” D. A. McLean and L. Ecrerron. This paper shows 
anthraquinone to be an effective stabilizer for capcitors having paper di- 
electrics containing chlorinated impregnants when aluminum electrodes are 
used and d-c. potentials are applied. One half per cent of anthraquinone 
prevents formation of the usual carbonized brown spots in the paper, and 
diminishes corrosion of electrodes and instability of leakage current. It 
increases the life under accelerated testing conditions by factors of four to 
one hundred fold, depending upon materials used and conditions of test. 
This development has added appreciably to the reliability of paper capacitors 
containing chlorinated impregnants, particularly for military equipment 
where high temperatures and high voltages are often encountered simulta- 
neously. Solubility of anthraquinone in the usual chlorinated impregnants 
is limited. Where greater solubility is desired, the more soluble chloro 
and methyl derivatives can be used. 

Reflex Oscillators.® J. R. Pierce. This paper discusses qualitatively 
the behavior of reflex oscillators. Power production, electronic tuning, 
variation of frequency with resonator voltage, effect of modulation coeffi- 
cient, and influence of load are considered. Two brief mathematical appen- 
dixes are included. 

Cape Charles-Norfolk Ultra-Short-Wave Multiplex System.2 N. F. 
ScHLaack and A. C. Dickreson. This paper describes the general features 
of a radio multiplex system which has been installed between Cape Charles 
and Norfolk, Virginia. The radio-frequency equipment operates in the 
vicinity of 160 megacycles. The system employs the 12 telephone channels 
of the type K cable carrier system which are in the frequency range 12 to 
60 kilocycles. 


6 Proc. I. R. E., February 1945. 
7 Indus. & Engg. Chem., January 1945. 
8 Proc. I. R. E., February 1945. 
® Proc. I. R. E., February 1945. 
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